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^Five  different  designs  of  broadband  piezoelectric  ultrasonic  transducers 
ire  evaluated  and  compared^to  determine  the  effects  of  various  transducer  desigr 
parameters  on  the  transducer  response.  The  transducers  tested  include  the  I 
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FC-300  except  that  they  J>a*e  backing  layers  of  0.0  cm  (no  backing),  1.270  cm 
(0.500  in),  and  1.905  cm  (0.750  in),  ^•espaptively.  Pairs  of  each  of  the  trans¬ 
ducers  are  clamped  face-to-face  and  tested  with  narrow-band  continuous  sinusoid 
and  tone  bursts  (0.4  to  3.0  MHz) ,  and  with  a  broadband  input  signal  (j.O  isen 
pulse)./-.  >4 

^jwdel  is  developed  which  includes  the  effects  of  the  generativn  sr.d 
detection"of  stress  waves  by  the  piezoelement,  the  propagation  of  wa^es  through 
attenuating  transducer  elements,  and  the  interaction  of  stress  waves  the 
interfaces  between  those  elements.  The  model  is  used  to  predict  the  output 
signals  for  each  of  the  two-transducers-together  systems,  and  the  predicts-: 
results  are  compared  to  the  measured  results  for  each  system. 

There  are  some  discrepancies  between  the  measured  and  predicted  c-esul::, 
These  discrepancies  are  eliminated,  however,  by  the  introduction  of  time  delays 
in  the  model  as  waves  interact  at  the  surfaces  of  the  piezoelement  ir  cur.  out¬ 
put  transducer.  For  the  transducers  with  titanium  carbide  wearplates ,  waves 
which  are  transmitted  across  the  wearplate-piezoelement  interface  are  delayed 
by  0.08-0.09  Msec,  and  reflected  waves  are  delayed  by  0.20-0.24  ysec.  (By 
comparison,  the  transit  times  of  the  wcarplate  and  piezoelement  are  0.04  and 
0.25  ysec,  respectively.)  These  time  delays  at  the  wearplate-piezoeler.ent 
interface  are  reduced  by  about  40%  when  the  titanium  carbide  wearplate  is 
replaced  by  an  aluminum  wearplate.  A  time  delay  is  also  introduced  for  waves 
traveling  through  the  backing  layer,  in  addition  to  the  transit  time  of  the 
layer  itself;  this  delay  ranges  from  0.75  ysec  for  the  AE-101  to  1.50  and  1.60 
ysec  for  the  AE-101  and  FC-500,  respectively. 

Two  additional  factors  are  shown  to  influence  the  correlation  between 
measured  and  predicted  output  signals.  First,  the  voltage  across  the  output 
piezoelement  decays  exponentially  with  time  due  to  the  flow  of  current  through 
the  oscilloscope.  The  time  constant  associated  with  this  decay  is  found  to  be 
1.0  ysec  for  the  FC-500,  and  about  5  ysec  for  the  other  transducers.  Secondly, 
there  is  a  multiplicative  amplitude  factor  which  depends  on  the  piezoelectric 
constant  and  capacitance  of  the  piezoelement,  which  determines  the  amplitude 
of  the  output  signal.  Since  the  piezoelectric  constant  and  capacitance  are 
difficult  to  measure,  the  amplitude  factor  is  evaluated  for  each  system  simply 
by  dividing  the  amplitudes  of  tb*'  measured  and  predicted  output  signals.  The 
amplitude  factor  is  found  to  be  fairly  consistent  from  one  transducer  to 
another;  over  the  frequency  range  0. 4-2.0  MHz,  the  amplitude  factors  for  all  of 
the  transducers  fall  between  3.5  and  7.0  x  1016  N/m3. 

Several  other  transducer  parameters  are  studied  to  determine  how  they 
influence  the  output  signal  for  a  particular  ultrasonic  test  system.  The 
effects  of  changes  in  backing  impedance,  backing  thickness,  wearplate  impedance 
and  piezoelement  thickness  are  evaluated,  and  a  set  of  guidelines  for  transduce 
design  is  developed.  In  addition,  the  role  of  couplant  layers  and  adhesive¬ 
bonding  layers  in  the  system  is  Investigated,  and  it  is  found  that  couplant 
layers  have  no  effects  on  wave  propagation  in  the  system  at  the  frequencies 
used  (0. 4-3.0  MHz). 
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CHAPTER  1 

INTRODUCTION  AND  LITERATURE  REVIEW 

1.1  INTRODUCTION 

Ultrasonic  testing  is  a  neans  of  assessing  the  mechanical 
integrity  of  a  structure  by  studying  the  propagation  of  high- 
frequency  stress  waves  through  that  structure.  An  ultrasonic 
test  is  performed  by  applying  a  known  stress  waveform  at  some 
point  on  the  structure,  and  then  measuring  the  resulting  stress 
either  at  that  same  location  or  at  some  other  location.  The 
relationship  between  the  input  and  output  stress  waveforms  can 
then  be  used  to  determine  the  mechanical  integrity  of  the 
structure. 

A  familiar  example  of  ultrasonic  testing  is  the  testing  of  a 
ping  pong  ball  for  cracks.  Bouncing  the  ball  on  the  table  intro¬ 
duces  a  stress  impulse  to  the  ball,  and  the  resulting  stress  waves 
in  the  ball  produce  sound  waves  which  propagate  to  the  ears  of 
the  observer.  With  experience,  one  can  detect  even  a  small  crack 
in  the  ball  by  listening  to  the  sound  it  makes. 

Of  course,  it  is  easy  to  detect  the  presence  of  overt  flaws 
in  a  specimen  by  such  a  method,  but  a  quantitative  determination 
of  strength  in  unflawed  and  marginal  specimens  is  much  more 
difficult.  Very  accurate  knowledge  of  the  input  and  output  stress 
waveforms  is  required  to  make  such  a  determination,  and  the 
expertise  to  interpret  the  results  of  ultrasonic  tests  in  terms 
of  the  strength  properties  of  the  specimen  is  needed  as  well. 
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Typically,  the  correlation  of  strength  properties  and 
ultrasonic  test  results  is  performed  by  fabricating  and  testing  a 
large  number  of  specimens  with  varying  types  and  degrees  of 
flaws.  The  results  of  the  ultrasonic  tests  are  compared  to  the 
results  of  destructive  mechanical  tests,  and  a  safe  limit  for 
certain  ultrasonic  testing  parameters  is  determined,  unless  a 
thorough  analysis  accompanies  the  experiments,  however,  these 
empirically-determined  "safe"  limits  can  be  misleading,  incom¬ 
plete,  or  even  incorrect.  Methods  must  be  developed  which  enable 
the  determination  of  strength  properties  directly  from  the  results 
of  ultrasonic  tests. 

There  is  clearly  a  great  need  for  improved  quantitative 
understanding  of  the  behavior  of  ultrasonic  test  systems  and  of 
the  significan^j^f  ultrasonic  test  results.  An  accurate  charac¬ 
terization  of  broadband  piezoelectric  ultrasonic  transducers 
is  a  crucial  step  toward  that  goal. 

1.2  REVIEW  OF  THE  ANALYSIS  AND  CHARACTERIZATION  OF  PIEZOELECTRIC 

ULTRASONIC  TRANSDUCERS 

Any  attenpt  to  characterize  an  ultrasonic  transducer  must 
begin  with  an  accurate  understanding  of  the  behavior  of  its 
piezoelement.  The  fundamental  equations  describing  the 
piezoelectric  effect  were  first  summarized  by  Voigt  in  his  early 
volume,  "Lehrbuch  der  Krystallphysik" ,  and  appear  in  a  slightly 
revised  form  in  Cady  [1],  A  significant  contribution  to  the  field 
was  made  by  Mason  [2]  with  the  discovery  that  the 
electromechanical  behavior  of  piezoelectric  crystals  could  be 


described  by  an  equivalent  electric  circuit,  with  one  electrical 
connection  and  two  mechanical  connections  corresponding  to  the 
pair  of  electrodes  and  the  two  faces  of  the  crystal,  respectively. 
Mason's  equivalent  circuit  has  since  been  employed  to  investigate 
many  transducer  phenomena:  the  coupling  between  different  natural 
modes  of  vibration  of  piezoelectric  crystals  [3] ;  the  effects  of 
this  so-called  mode  conversion  along  with  impedance-matched 
backing  materials  and  frequency- tuned  wear  plates  on  the  frequency 
response  and  bandwidth  of  transducers  [4,5,6];  and  the  nature 
of  the  overall  transmission  coefficient  of  layered  transducers 
[7],  including  the  effect  of  adhesive  bonding  between  layers  [9]. 
Hie  strength  of  Mason's  model  lies  in  the  application  of 
straightforward,  well-understood  methods  of  circuit  analysis  to 
an  otherwise  very  difficult  electromechanical  problem. 
Unfortunately,  many  of  the  equivalent-circuit  parameters  are  not 
easily  measured,  and  it  can  be  difficult  even  to  estimate  values 
for  these  quantities  [9,10].  In  addition.  Mason's  model  provides 
little  if  any  intuitive  insight  into  the  mechanical  aspects  of 
transducer  operation. 

Other  electctrical  models  have  been  developed  to  facilitate 
transducer  analysis.  A  Thevenin-equivalent  circuit  approach  has 
been  suggested  [11],  and  a  transmission-line  model  [12]  has  been 
employed  by  Desilets,  et.  al.  [13]  to  evaluate  broadband  transducer 
efficiency,  and  by  Dotti  [14]  to  predict  impulse  response. 

Despite  such  widespread  use,  however,  the  value  of  equivalent 
circuits  in  understanding  mechanical  phenomena  is  quite  limited. 

Sane  work  has  indeed  been  done  toward  predicting  transducer 
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behavior  from  Che  principles  of  applied  mechanics.  Linear 
elasticity  [15,16]  and  computational  finite-element  methods  [17] 
have  been  applied  to  analyze  the  natural  modes  of  piezoelectric 
plate  vibration.  The  two-dimensional  wave  equation  has  been  used 
as  a  starting  point  for  at  least  one  characterization  [18].  These 
analyses  are  certainly  of  value  in  promoting  physical 
understanding,  but  provide  little  aid  in  describing  the 
interaction  of  the  transducer  components  with  one  another. 

Among  the  most  significant  work  in  predicting  the  actual 
electromechanical  behavior  of  a  transducer  system  is  that  of 
Redwood  [19,20].  Concerning  earlier  work  based  upon  the 
simplified  Mason  equivalent  circuit,  Redwood  wrote: 


"This  equivalent  circuit  is,  however,  an  approximation, 
being  valid  only  at  frequencies  near  the  mechanical  resonance 
of  the  transducer,  and  can  therefore  be  applied  only  to 
problems  involving  continuous  waves,  and  even  then  only  over 
a  very  limited  band  of  frequencies.  This  circuit  is 
certainly  of  no  value  in  the  calculation  of  transient 
performance,  and  here  it  is  necessary  to  use  the  exact 
equivalent  circuit.  The  form  of  the  exact  circuit  is  knowi, 
and  its  derivation  can  be  found  in  the  literature,  but  it  has 
rarely  been  discussed  or  made  use  of.  It  incorporates  a 
section  of  a  transmission  line,  and  the  basic  purpose  of  this 
is  to  represent  the  time  delay  which  is  necessary  for 

mechanical  signals  to  travel  from  one  face  of  the  transducer  j 

to  the  other  face."  [19] 


Redwood  observed  that  the  shortcomings  of  much  of  the  earlier  work 
with  resonant  transducers  were  due  largely  to  the  simplifying 
assumptions  that  had  restricted  both  the  bandwidth  and  the 
usefulness  of  Mason's  equivalent  circuit.  In  fact,  the  above- 
referenced  transmission-line  models  were  no  more  than  redevelopnents 
of  that  portion  of  Mason's  model  which  had  been  neglected  in  the 
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succeeding  resonant-transducer  analyses,  revived  because  of  the 
growing  interest  in  broadband  transducers  in  recent  years. 

Redwood  presents  an  excellent  hybrid  electrical-mechanical 
analysis  of  transducer  response,  including  both  the  resonant 
effects  from  the  simplified  Mason  model  and  the  time  delays  from 
the  exact  circuit.  The  technique  developed  in  [19]  is  then 
utilized  in  [20]  to  predict  transducer  response  in  the 
transmission  and  reception  of  pulses,  and  in  particular  that 
portion  of  the  response  caused  by  the  superposition  of  reflected 
waves  produced  at  the  faces  of  the  piezoelement.  This  approach 
serves  as  the  basis  for  the  present  work. 

Some  additional  physical  significance  may  be  given  to  the 
transmission-line  model  by  considering  the  actual  phenomenon  which 
it  describes — the  behavior  of  stress  waves  at  impedance-mismatched 
interfaces.  Altman  and  Beyer  [21]  consider  the  reflections 
generated  at  the  faces  of  the  piezoelement  in  a  typical 
transducer,  and  Hill  and  El-Dardiry  [22]  eaploy  a  chain-matrix 
computation  to  determine  the  frequency-dependent  transmission 
coefficient  for  a  five-layer  transducer.  Folds  and  Loggins  [23] 
treat  the  oblique  incidence  of  waves  on  transducer  surfaces,  and 
Doll  [24]  utilizes  these  concepts  to  predict  the  time  delays  and 
relative  electrical  signal  amplitudes  in  an  ultrasonic 
through-transmission  test  system.  Lee  and  Williams  [25]  consider 
such  a  system  and,  assuming  a  thin  specimen,  obtain  a  closed-form 
solution  for  the  steady-state  output  signal  amplitude  resulting 
from  an  infinite  number  of  reflections  of  stress  waves  at  the 
transducer-specimen  interfaces.  This  work  provides  the  tools  for 


an  empirical  determination  of  the  maximum  ratio  of  output  voltage 
to  input  voltage  for  a  given  through-transmission  system.  This  is 
an  important  step  toward  the  determination  of  the  transduction 
ratio  (stress  produced  in  the  specimen  per  unit  voltage  introduced 
to  the  transducer,  and  vice-versa)  for  a  particular  transducer- 
specimen  system. 

In  addition  to  the  work  in  transducer  analysis,  much  effort 
has  also  been  devoted  to  better  understanding  the  character  of 
the  ultrasonic  stress  field  produced  by  the  transducer.  Topics 
of  interest  in  this  area  include  the  effect  of  the  presence  of  the 
specimen  on  the  frequency  response  of  the  transducer  [26];  the 
effect  of  piezoelement  boundary  conditions  [27]  and  the  wearplate 
geometry  [28]  on  the  directionality  of  the  stress  field;  the 
limitations  of  the  assumption  of  plane-wave,  piston-source  wave 
propagation  within  the  transducer  [29] ;  and  the  usefulness  of 
finite-element  methods  in  predicting  the  resulting  pressure 
distribution  in  the  specimen  (30). 

As  it  becomes  increasingly  desirable  to  attach  quantitative 
significance  to  ultrasonic  measurements,  many  techniques  are  being 
proposed  for  the  evaluation,  calibration,  and  standardization  of 
ultrasonic  transducers.  Among  these  are:  the  use  of  reciprocity 
methods  to  determine  the  transducer  sensitivity  in  a  particular 
application  [31,32];  the  theoretical  [33]  and  experimental  [34,35] 
determination  of  loaded-  and  unloaded- transducer  frequency 
response;  and  the  experimental  measurement  of  the  distribution  of 
ultrasonic  stress  in  the  near  field  of  the  transducer  [36] . 
Methods  for  determining  the  integrity  of  the  transducer-specimen 


bond  have  also  been  developed  [37].  A  good  treatment  of  current 
transducer  evaluation  procedures  may  be  found  in  [38]. 

1.3  SCOPE  OF  THE  PRESENT  WORK 

The  previous  section  shows  that,  although  some  work  has  been 
done  in  the  area  of  transducer  analysis  and  characterization, 
much  of  the  work  uses  simplified  electrical  equivalent  circuits 
and  is  of  little  value  in  gaining  quantitative  insight  into  the 
mechanical  aspects  of  the  behavior  of  ultrasonic  test  systems. 
Although  some  have  studied  particular  mechanical  topics  such  as 
the  propagation  of  stress  waves  through  layers  of  dissimilar 
materials,  there  is  still  the  need  to  develop  a  better  quantita¬ 
tive  understanding  of  the  effect  of  transducer  geometry  and 
material  properties  on  the  behavior  of  the  transducer  as  a  unit 
and  as  part  of  a  system. 

The  present  work  addresses  this  need.  A  theoretical  model 
of  a  multilayer  broadband  piezoelectric  ultrasonic  transducer 
will  be  developed,  based  on  the  physics  of  the  piezoelectric 
effect  and  on  the  theory  of  wave  propagation  in  elastic  materials. 
Experiments  will  be  performed  on  transducers  with  different 
geometries  and  material  properties  to  determine  the  effects  of 
these  factors  on  the  input-output  characteristics  of  the  trans¬ 
ducers.  The  results  of  the  experiments  will  be  compared  to  the 
predicted  results  from  the  model  in  order  to  gain  additional 
insight  into  the  operation  of  piezoelectric  transducers  and 
ultrasonic  test  systems. 
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CHAPTER  2 

ANALYSIS  OF  ULTRASONIC  TEST  SYSTEM  AND  DEVELOPMENT  OF  MODEL 

2.1  INTRODUCTION 

A  schematic  representation  of  a  typical  ultrasonic  through- 
transmission  test  system  is  shown  in  Fig.  1.  Notice  that  each 
transducer  consists  of  three  components:  (1)  the  piezoelement, 
which  converts  stress  to  voltage  and  vice-versa?  (2)  the 
wearplate,  which  protects  the  piezoelement  from  direct  contact 
with  the  specimen;  and  (3)  the  backing  layer,  which  absorbs 
stress  waves  from  the  piezoelement.  In  addition  to  these  three 
components,  each  transducer  is  potted  with  a  urethane  plastic. 

The  function  of  the  potting  is  to  protect  the  wires  leading  from 
the  piezoelement  to  the  connector,  but  it  is  of  interest  here 
because  it  has  finite  impedance,  so  that  the  "free"  surface  of 
the  backing  layer  is  not  entirely  free.  There  is  also  a  thin 
layer  of  couplant  between  the  transducers  and  the  specimen,  which 
improves  the  transmission  of  stress  waves  across  these  interfaces. 

An  ultrasonic  test  is  initiated  by  applying  a  known  time- 
varying  voltage  signal  to  the  input  transducer.  This  voltage 
produces  a  stress  waveform,  vfiich  propagates  through  the  specimen 
and  into  the  receiving  transducer.  The  receiving  transducer  then 
converts  the  stress  waveform  back  into  an  output  voltage  signal, 
which  is  measured  and  compared  to  the  input  voltage.  These  three 
processes— generation,  propagation,  and  detection  of  stress  waves — 
are  analyzed  separately  in  the  following  sections,  and  the 
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results  are  then  combined  to  produce  a  complete  system  model. 


2.2  NOMENCLATURE 

In  order  to  keep  track  of  the  enormous  number  of  stress  waves 
which  are  produced  within  the  system,  names  have  been  assigned  to 
the  waves  according  to  the  following  rules  (see  Pig.  2): 


(1)  The  input  transducer  is  designated  transducer  1,  and 
the  output  transducer  is  called  transducer  2. 

(2)  The  backing,  piezoelement,  and  wearplate  of  each 
transducer  are  designated  B,  p,  and  w,  respectively. 
The  specimen  is  denoted  by  the  letter  S. 

(3)  Waves  traveling  to  the  right  are  indicated  by  R;  left¬ 
traveling  waves  are  named  L. 

(4)  Waves  arriving  at  the  interfaces  between  transducer 
components  or  at  the  wearplate-specimen  interfaces 
have  names  ending  in  I;  waves  which  are  generated  at 
or  reflected  from  these  Interfaces  and  propagate  away 
from  them  are  indicated  by  the  letter  0. 


These  names  represent  the  instantaneous  amplitudes  of  the 
particular  waves  to  which  they  correspond.  For  example,  P^LO 
refers  to  one  of  the  waves  at  the  piezoelement-backing  interface 
in  the  receiving  transducer,  as  indicated  by  the  first  part  of 
the  name,  ?2 B2.  The  rest  of  the  name,  LO,  says  that  the  wave 
travels  away  from  the  interface  and  to  the  left.  Similarly,  wave 
W^R  would  be  the  right- traveling  wave  in  the  wearplate  of  the 
input  transducer. 

Some  of  the  other  variables  that  appear  in  the  analysis  are 
listed  below,  along  with  the  SI  units  of  each: 
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Lg  *  thickness  of  backing  layer  (cm) 

Lp  *  thickness  of  piezoe lenient  (an) 

Ly  *  thickness  of  wearplate  (cm) 

L  *  thickness  of  specimen  (an) 

S 

2 

Z  *  acoustic  impedance  of  backing  (gm/an  /sec) 

B 

2 

2^  *  acoustic  impedance  of  piezoelement  (gm/an  /sec) 

2 

1  *  acoustic  impedance  of  wearplate  (gm/an  /sec) 

w 

2 

Z  »  acoustic  impedance  of  specimen  (gm/an  /sec) 

s 

c  a  wavespeed  in  backing  material  (an/sec) 

B 

cp  *  wavespeed  in  piezoelement  (an/sec) 

c ,  *  wavespeed  in  wearplate  (an/sec) 

w 

cg  «  wavespeed  in  specimen  (cm/sec) 

C0  *  capacitance  of  piezoelement  (farads) 
h  »  piezoelectric  constant  (V/m) 

2 

S  »  cross-sectional  area  of  piezoelement  (an  ) 

R  «  input  resistance  of  oscilloscope  (ohms) 
a  (u>)  »  attenuation  of  backing  material  (nepers/an) 
ag(^)  *  attenuation  of  specimen  (nepers/cm) 

The  values  of  some  of  these  parameters  vary  among  the  prototype 
transducers  tested  in  this  study.  Some  of  the  transducers  have 
titanium  carbide  wearplates,  for  example,  and  some  have  alininin 
wearplates.  It  is  possible,  therefore,  for  the  input  and  output 
transducers  in  the  same  system  to  have  different  wearplate 
impedances,  in  which  case  these  impedances  are  denoted  by  Z 

W1 

and  Z,  ,  according  to  the  first  rule  given  earlier.  No  such 
distinction  is  made  for  variables  which  have  the  same  value  in 
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every  transducer,  such  as  piezoelement  thickness  and  backing 
impedance. 

2.3  GENERATION  OF  STRESS  WAVES  FROM  INPUT  VOLTAGE 

Applying  a  step  function  of  voltage  across  the  piezoelement 
of  an  ultrasonic  transducer  produces  four  stress  waves  (see 
Appendix  A,  section  A.7) .  Using  the  above  nomenclature,  the 
amplitudes  of  these  four  waves  are  given  by  eqns.  (A82)  in 
Appendix  A  as: 


P.V^ro 


P1W1L0 


B1P1R0 


B1P1L0 


(2-1) 


where  h  is  the  piezoelectric  constant  in  V/m,  C0  is  the  electrical 
capacitance  in  farads,  S  is  the  surface  area  of  the  piezoelement 
in  an  ,  and  V0  is  the  amplitude  of  the  step  input  of  voltage.  The 
value  of  the  time  constant,  x,  is  given  by  eqn.  (A79)  in  Appendix  A: 


,  h2C2 

I  „  0 

Ih + 

* 

[VLi"dl 

T  S 

l1 

(ZB  *  Zp 

w 

p  +  Z*L 

i 

Eqns.  (2-1)  give  the  amplitudes  of  the  four  stress  waves 
produced  in  the  input  transducer  by  a  step  input  of  voltage. 
These  may  be  represented  in  a  general  form  as: 


I  -  KVo  exp  (i) 


(2-3) 


where  T  is  stress,  and  K  is  a  constant.  In  order  to  extend  this 
step  response  to  a  continuous,  time-varying  input  voltage  v(t) , 
a  stepwise  approximation  must  be  used,  as  shown  in  Fig.  3.  This 
translates  V(t)  into  a  series  of  n  discrete  steps,  each  of 
amplitude  AV(i),  which  are  separated  in  time  by  At,  known  as  the 
time  interval  of  the  discretization.  If  each  of  these  steps  AV(i) 
produces  a  stress  response  such  as  that  given  in  eqn.  (2-3),  the 
stress  at  any  time  t*kAt  is  simply  the  superposition  of  the  stress 
responses  produced  by  all  of  the  preceding  voltage  steps: 


T(kAt)  -  k|AV(1)  exp |(k  -  1)  ^r} 

+  AV (2)  exp  |  (k  -  2)  ^}+  •  •  • 

+  AV(k  -  1)  expj^p-j  +  Av(k)l 

*  (2-4) 

It  is  convenient  to  write  the  series  as  a  summation: 


k 

T(kA7)  -  K  l  AV(i)  exp 
i*l 

(2-5) 

This  sunnation  gives  the  time-dependent  amplitude  of  a  stress 
wave  generated  at  the  piezoelement  of  the  input  transducer  for  an 
arbitrary  time-varying  input  voltage,  V(t).  Evaluating  the 
constant  K  in  eqn.  (2-3)  for  each  of  the  four  stress  waves 


given  in  eqns.  (2-1) ,  and  substituting  into  eqn.  (2-5)  gives: 


PjWjROGcAt) 


P1W1L0(kAt) 


~?W)  ^  l.  u>  “P  I 

>  wl/  s 

(  ZP  \  hC  r 

z  +z"  J— 

\p  wi/  s  i-i 


(k  -  i)  ~- 


V li)  exp  { \k  -  i) 


B1P1R0(kAt)  -  L-^-y  )  S,  £  AV(i)  exp  j  (k  -  i)  ^  [ 
\  B  p/  S  i"l  ' 


B1P1L0(kAt)  -  .z 


T~P~^  hCo  i  exp  |(k  -  i)  rr 

>  P/  S  »-l  '•  1 


where  the  stepwise  function  V(i)  is  obtained  from  V(t)  by: 


(2-6) 


AV(i)  -  V(iAt)  -  V((i  -  l)At) 


(2-7) 


Note  that  the  anplitudes  given  in  eqns.  (2-6)  are  solely  those 
resulting  from  the  application  of  a  voltage  to  the  piezoelement. 
If  other  waves  are  already  propagating  in  the  system,  their 
amplitudes  must  be  added  to  those  in  eqns.  (2-6)  to  obtain  the 
actual,  tine-varying  amplitudes  of  waves  P  W  RO,  P  K  LO,  B  p  ac, 

11  1  u.  11 

and  B^LO. 


2.4  PROPAGATION  OF  STRESS  WAVES  WITHIN  THE  SYSTEM 

As  soon  as  the  stress  waves  are  produced  at  the  piezoelement 
of  the  input  transducer,  they  begin  to  propagate  through  the  system. 
As  they  propagate,  they  are  affected  in  several  ways:  (1)  time 
delays  are  introduced  due  to  the  wave  transit  times  across  the 
transducer  and  specimen  layers;  (2)  a  decrease  in  amplitude  shows  up 
due  to  attenuation  in  certain  layers;  (3)  reflections  are  generated 
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at  impedance-mismatched  interfaces  between  layers.  These  effects 
are  treated  separately  below. 


2.4.1  Wave  Propagation  in  Piezoelement  and  Wearplate 

The  velocity  of  propagation  of  plane  compressional  waves  in 
both  piezoelectric  and  nonpiezoelectric  materials  is  given  by  eqn. 
(A5)  in  Appendix  A  as: 

2  E 
c  -  — 

p  (2-8) 
where  c  is  the  wavespeed  in  m/s,  E  is  the  elastic  modulus  of  the 
material  in  N/m  ,  and  p  is  the  mass  density  in  kg/mJ.  If  a  layer 
is  nonattenuating  and  norvdispersive,  it  propagates  waves  at  a 
constant  velocity  without  changing  their  shape.  Ttie  only  effect 
on  the  waves  is  the  introduction  of  a  time  delay,  \,  where: 


Here,  X  is  the  so-called  "transit  time"  of  the  layer  in  seconds, 

L  is  the  thickness  of  the  layer  in  meters,  and  c  is  the  wavespeed 
in  m/sec. 

Attempts  have  been  made  [39J  to  determine  the  attenuation  of 

the  lead  metaniobate  piezoelement  and  the  titaniun  carbide 

wearpl  ite  used  in  the  prototype  transducers.  Comparatively  no 

attenuation  has  been  observed  for  the  wearplate,  and  negligible 

attenuation  has  been  reported  for  the  piezoelement  (see  Fig.  C6) . 

The  piezoelement  and  wearplate  may  therefore  be  considered  to  be 

nonattenuating  for  the  purposes  of  this  study,  so  that  if  X  and 

1 

Xy  represent  the  transit  times  of  wearplates  1  and  2  as  found 
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from  eqn.  (2-9) ,  and  if  Xp  represents  the  piezoelement  transit 
time,  the  following  equations  hold  for  the  piezoelements: 


P  V  Rl(t  +  Xp)  -  B1P1RO(t) 
B1P1LI(t  +  Xp)  -  P^LOCt) 
P2B2RI(t  +  Xp)  »  W^P2R0(t) 
W2P2LI(c  +  Xp)  =  P2B2LO(t) 

and  for  the  wearplates: 


(2-10) 


W  SRI(t  +  x,t  )  -  P ,M,RO(t) 

1  1  1 

p1w1T,i(t  +  \  )  m  v1-0^) 

W2P2RI(t  +  )  -  SW2RO(t) 

SW2LI(t  +  X^,  )  -  W2P2LO(t) 

2  (2-11) 

Eqns.  (2-10)  and  (2-11)  account  for  the  time  delays  introduced  by 

the  nonattenuating  piezoelements  and  wearplates,  respectively. 

For  example,  the  first  of  eqns.  (2-11)  shows  that  the  wave  leaving 

the  piezoelement-wearplate  interface  in  the  input  transducer  at 

time  t,  P^W^RO(t) ,  arrives  at  the  interface  between  wearplate  and 

specimen  (W.SRI)  at  time  t+X  ,  that  is,  after  the  transit  time 

l  w: 

required  to  pass  through  the  wearplate. 


2.4.2  Wave  Propagation  in  Backing  Material  and  Specimen 

In  addition  to  introducing  a  time  delay,  some  materials  also 
attenuate  the  waves  which  pass  through  them,  thereby  reducing  their 
amplitude.  Fig.  C6  shows  that  the  backing  material  used  in  the 


-29- 


prototype  transducers  exhibits  significant  attenuation  at  all 
frequencies,  with  a  very  strong  frequency-dependence  above  1.0  MHz. 

The  attenuation  of  stress  waves  in  material  layers  is 
described  by  the  following  equation: 

T 

”  exp{-a( 

(2-12) 

where  is  the  amplitude  of  the  stress  wave  that  enters  the  layer, 
T0  is  the  amplitude  of  the  exiting  wave,  a(v)  is  the  frequency- 
dependent  attenuation  parameter  in  nepers/m  (where  w  is  the  radian 
frequency) ,  and  L  is  the  thickness  of  the  layer  in  meters. 

Using  eqns.  (2-9)  and  (2-12),  the  equations  describing  wave 
propagation  through  the  attenuating  backing  layers  may  be  written: 

B  P  RI ( t  +  X_  )  -  B.RO  (t)  exp  (-a_(u)L.  ) 

11  1  o  d 

B..LI  (t  +  X  )  *  B  P.LO(t)  exp(-=t  (u.')L  ) 

1  O  11  D 

B.RI  (t  +  X  )  *  P  B  RO(t)  exp(-a  (w)L  ) 

2  o  Oj  2  2  B  is ^ 

P  B  LI(t  +  XR  )  -  B  LO  (t)  exp  (-a  (w)L  ) 

2  2  b2  2  o  B  B2  (2-13) 

Here,  (w)  is  the  attenuation  of  the  backing  material  in  nepers/m, 

B 

X  and  X  .  represent  the  transit  times,  in  seconds,  of  the  backing 

W1  W2 

layers,  and  and  L„  are  the  thicknesses  of  the  layers  in  meters. 
The  wave  names  subscripted  0  denote  those  waves  at  the  rear  surface 
of  the  backing  layer,  that  is,  tne  surface  farthest  away  from  the 
piezoelement. 

In  general,  the  specimen  will  also  exhibit  frequency-dependent 
attenuation,  with  attenuation  parameter  as(w).  If  the  impedance. 


wavespeed,  and  thickness  of  the  specimen  are  known,  eqn.  (2-9)  can 
be  used  to  obtain  the  transit  time  >  for  the  specimen  and  the 
following  propagation  equations  may  be  written: 

SW2Rl(t  +  \  )  -  W  SR3(t)  exp  t-ua(t.)Ls) 

W,  SLI  t  +  X.)  *  SVLLO(t)  exp  (-a  (u)L  ) 

1  S  S  a  (2—14) 

Here,  as  in  eqns.  (2-i7>,  ;c  <  i«  tna  attenuation  of  the 

specimen  in  nepers/m,  X  is  the  transit  time  of  the  specimen  layer 

O 

in  seconds,  and  Lg  is  the  thickness  of  the  specimen  in  meters. 


2.4.3  Reflection  at  Impedance-Mismatched  Interfaces 

Whenever  a  wave  strikes  an  interface  between  two  materials 
of  different  acoustic  impedance,  a  portion  of  its  amplitude  is 
transnitted  across  the  interface,  and  a  portion  is  reflected  back 
(see  Appendix  A,  section  A.4).  If  a  stress  wave  with  amplitude  Tj 
traveling  in  a  mediun  with  impedance  Z^  encounters  a  medium  with 
impedance  Z2»  eqns.  (A31)  in  Appendix  A  give  the  amplitudes  of  the 
resulting  reflected  and  transnitted  waves: 


T 

tr 


2  Z. 


Z,  +  z„ 


(2-15) 

where  Ttr  represents  the  wave  transmitted  across  the  interface, 
and  Tr  represents  the  reflected  wave. 

Applying  eqns.  (2-15)  to  the  impedance-mismatched  interfaces 


between  wearplates  and  piezoelements  gives  the  following: 
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P^RO 


P1W1LO  - 


w2p2ro 


W-P.LO 
2  2 


2  2^  (P^RI)  +•  (2p  -Z^XP^LI) 
(ZP  +  Zv^ 

(ZM  -  ZpKP^RI)  +  2  Zp  (P^LI) 

“  <ZP 

2  ZP  (W2P2RI)  +  (2W  ‘  ZP)(W2P2LI) 

(zw2  +  ZP) 

(Zp  -  Zw  )(W2P2RI)  +  2  Zw  (W2P2RI) 

1  (V  v  2 


(2-16) 


At  the  wearplate-specimen  interfaces,  the  equations  are: 


I^SRO  « 

VT^SLO  = 
SK,RO  - 

s:;2lo  = 


2  Z  (V  SRI)  +  (ZW,  "  Z  )(W_SLI) 

SI  1  S  i-  _ 


(ZW,  +  z  ) 
"1  s 

(Zs  - 

•  zw^a^sRi)  + 

2  Zv^SLI) 

<ZKj  *  Zs> 

2  ZV;, 

(SV0RI)  +  (Z  - 
z  s 

Zw^(SW2LI) 

(ZS  +  z^2) 

(ZWt 

-  z  )  (SW -,RI)  + 

S  a. 

2  Z  (SW.LI) 
s  ^ 

(zs  +  Zw.) 

5  L. 


And  at  the  rear  face  of  the  backing  layer: 


(2-17) 


B.RO  *  (  U-  B)  B  LI 
10  \Zu  +  ZB/  1  ° 


B2L0o  '(z“  *  z“)B2RIo 


(2-ie) 


-V  Sf* ’ 


Note  that  each  of  the  waves  produced  at  an  interface  depends  on 
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bo  th  of  the  waves  arriving  at  that  interface.  For  example,  the 
first  of  eqns.  (2-16)  says  that  the  right-traveling  wave  leaving 
the  interface  between  piezoelement  and  wearplate  in  the  input 
transducer,  P^RO,  is  the  sun  of  the  transmitted  part  of  the 
incident  right-traveling  wave  at  that  interface,  P  W  RI,  and  the 
reflected  part  of  the  left-traveling  wave,  P^WjLI.  For  simplicity, 
the  time-dependence  of  these  quantities  is  not  explicitly  stated 
here,  although  they  do  represent  instantaneous  values  of  continuous 
time-dependent  stresses. 


2.5  DETECTION  OF  STRESS  WAVES  AND  PRODUCTION  OF  OUTPUT  VOLTAGE 
The  output  voltage  Vout  produced  by  a  step  input  of  stress 
in  the  wearplate  of  a  transducer  is  obtained  from  eqn.  (A57)  in 
Appendix  A: 


V 

out 


(t) 


-2  h  -ft) 

(zW2  +  y  I  ~ — l 

l  l2  1 


(2-19) 


where  W2P2RI°  is  the  amplitude  in  N/m2  of  the  step  input  of 
stress,  and  R  is  the  resistance  in  ohms  connected  across  the 
piezoelement  of  the  transducer.  The  values  of  the  time  constants 
t1  and  t2  are  given  in  eqn.  (A56) : 


_1_ 

T1 


(2-20) 


where  CQ  is  the  capacitance  of  the  piezoelement  in  farads,  h  is 

the  piezoelectric  constant  in  V/m,  and  Zu  and  Z.  are  the 

2  V 


3 


characteristic  impedances  of  the  vearplate  and  piezoelement , 
respectively. 

It  will  later  prove  useful  to  write  the  output  voltage  in 
terns  of  the  stress  waves  inside  the  piezoelement.  If  M^p^Rl 
is  the  only  stress  wave  entering  the  transducer,  eqn.  (2-16) (b) 
gives 


2  w2p?ri  y^RG 


(2-21) 

Making  this  substitution  for  V^P.RI  in  eqn.  (2-19)  gives  the 
output  voltage  in  terms  of  stress  wave  W^RO,  inside  the 
piezoelement: 


(2-22) 

The  term  in  brackets  in  this  equation  is  a  complicated 
function  of  time,  especially  for  values  of  t  and  x2  which  are 
close  to  one  another.  A  simplifying  assumption  may  be  made, 
however.  Eqn.  (2-20)  shows  that  the  difference  between  the  time 

i 

1  ““ 

constants  depends  on  [hV[RS(Z  +Z  )]}2,  which  is  estimated  to  be 

W2  P 

an  order  of  magnitude  less  than  1/RC  ,  for  the  transducers  used  in 
this  study.  Neglecting  this  term,  in  eqn.  (2-20)  then  gives 
l/Tj*l/RC0,  and  1/t2*0.  Substituting  these  values  back  into  eqn. 
(A55)  and  reinverting  the  Laplace  transform  [40]  produces  the 


final  result: 
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V 

out 


(t) 


-  h  W,P,RO 
2  2  o 


Z 

P 


(2-23) 


This  equation  gives  the  output  voltage  in  terms  of  the  stress 
inside  the  piezoelement.  Although  considerably  simpler  than  the 
previous  equation,  it  nonetheless  shows  some  interesting  behavior. 

First  observe  that  if  the  time  constant  RC0  is  large,  the 
term  exp(-t/RC0)-*l,  and  the  integral  is  simply  equal  to  t.  This 
corresponds  to  a  ramp  of  voltage  for  a  step  input  of  stress,  so 
that  the  voltage  is  simply  the  time-integral  of  the  stress  inside 
the  piezoelement.  Fig.  4  shows  a  rectangular  pulse  of  stress  as 
it  propagates  through  the  piezoelerrent,  and  the  voltage  which  is 
generated  as  a  result.  Note  from  Fig.  4  that  the  voltage  is 
determined  only  by  the  stress  inside  the  piezoelement — when  the 
stress  wave  leaves  the  piezoelement,  the  voltage  returns  to  zero. 

Secondly,  note  what  happens  if  FCo*0.  This  can  only  be 
realized  if  R=C,  corresponding  to  t  short  circuit  across  the 
piezoelement.  The  exponential  in  eqn.  (2-23)  then  goes  to 
exp(-  »)«3,  so  that  vout  is  always  0.  No  voltage  is  accumulated 
across  the  piezoelement  because  of  the  short  circuit. 

The  third  case  is  the  most  interesting,  where  RCq  is 
somewhere  between  i!  and  «■ .  For  this  case,  the  voltage  is  also  an 
integral,  but  it  is  not  precisely  the  integral  of  the  stress 
waveform.  It  is,  instead,  the  convolution  of  the  stress  waveform 
with  an  exponentially-decaying  function  of  time.  A  closed-form 
solution  of  the  integral  gives: 
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(2-24) 


so  that,  rather  than  a  ramp,  the  output  voltage  shows  an  inverted 
exponential  rise  converging  to  a  maximum  value  (see  Fig.  5).  At 
the  maximum,  a  steady-state  condition  exists,  where  the  decay  of 
voltage  is  exactly  compensated  by  the  increase  in  voltage  due  to 
the  additional  stress  continually  entering  the  piezoelement. 

Thus,  the  voltage  across  the  output  resistance  R  is  a  rather 
complicated  function  of  the  stress  inside  the  piezoelement.  Eqn. 
(2-23)  indicates  that  waves  entering  the  piezoelement  produce  a 
negative  voltage,  so  that  those  waves  leaving  the  piezoelement 
must  conversely  produce  a  positive  voltage.  Including  the  effects 
of  all  of  the  waves  entering  and  leaving  the  piezoelement  then 
gives  the  total  output  voltage: 

r 


-  h  RC 


V  (t) 
out 


W,P,RO 
2  2  o 


VJ  P  LI 
2  2  o 


-  P  B  RI  +  P,B,LO  ► 
2  2  o  2  2  o 


(2-25) 


Here  the  subscript  o  denotes  the  amplitude  of  a  step  function  of 
stress,  arriving  at  the  surface  of  the  piezoelement  at  time  t=0. 

An  equation  analogous  to  eqn.  (2-7)  can  be  used  to  obtain  a 
ste^ise  approximation  of  each  of  the  continuous  time-varying 
stress  waves  which  contribute  to  the  output  voltage.  The  super¬ 
position  metfiod  of  section  2.3  then  gives  the  total,  time-varying 
output  voltage  resulting  from  the  combination  of  all  of  these 


stress  waves: 
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Vout(k*c>  * 


-  h  RC  k 

Z  ~  L' 
P  i-1 


AW2P2RO(i)  -  AW2P2LI(i) 


(2-26) 


This  equation  is  used  in  the  model  to  obtain  the  output  vcitage 
produced  by  a  known  set  of  stress  waves  W'2P2RO(t),  W2P  LI(t), 
P2B2RI<t)#  and  P2B2LO(t). 


2.6  DEVELOPMENT  OF  OVERALL  SYSTEM  MODEL 

The  equations  of  the  previous  three  sections  describe  all 
phases  of  the  conversion  from  input  voltage  to  output  voltage  in 
an  ultrasonic  test  system.  All  that  remains  is  to  combine  these 
equations  into  a  sufficiently  general,  yet  convenient  form  for 
application  to  particular  systems.  This  is  a  straightforward 
task,  except  for  the  bookkeeping  involved  in  tracking  the  many 
waves  generated  within  the  system.  This  section  introduces  the 
algorithm  used  to  solve  this  problem. 

In  Section  2.3  it  was  shown  that  each  specimen  and  transducer 
layer  has  a  characteristic  transit  time  \  ,  given  by  eqn.  (2-9) . 

If  the  transit  times  of  all  of  the  layers  were  the  same,  this 
transit  time  could  be  used  as  the  time  step  At  for  the  discretiz¬ 
ation  of  the  input  signal  (see  eqn.  (2-7)),  so  that  it  would  be  a 
simple  matter  to  account  for  all  of  the  waves  in  the  system  by 
applying  the  interface  equations  at  times  At,  2At,  3At,  and  so  on. 
Each  of  the  individual  transit  times  is  different,  however,  so 
that  a  few  multiple  reflections  can  generate  a  very  large  number 
of  waves  which,  rather  than  superimposing  upon  one  another,  fall 


between  one  another  in  time.  The  number  of  waves  produced 
increases  as  time  passes,  so  that  for  long  simulations  literally 
thousands  of  waves  must  be  considered. 

A  simple  solution  is  to  mathematically  discretize  each  of  the 
transducer  and  specimen  layers  into  snail  segnents  of  equal  transit 
time  X0,  as  shown  in  Fig.  6.  This  transit  time  becomes  the  time 
step  for  the  simulation,  and  may  be  arbitrarily  chosen  as  small  as 
desired.  If,  for  example,  the  piezoelement  has  53  segments,  it  is 
a  simple  matter  for  the  digital  computer  to  keep  track  of  the 
amplitudes  of  the  stress  waves  at  each  of  the  49  locations  between 
the  segments,  and  to  shift  the  amplitudes  from  one  location  to 
another  to  represent  an  elapsed  time  of  one  transit  time  \Q.  For 
instance,  recall  the  first  of  eqns.  (2-10): 


P^RKt  +  \p)  -  B^ROa) 


(2-27) 


If  the  array  p^RCi}  contains  the  amplitudes  of  the  right-t  rave  ling 
stress  waves  between  the  segments  of  the  input  piezoelenent,  and 
A0  is  the  transit  time  of  each  segment,  eqn.  (2-22)  becomes  a 
system  of  equations: 


h:R{i}(t  +  A0)  *  b1P  R0(t) 
P^R;2}(c  +  Xc)  -  piR{i}(t) 
P:R(3}(t  +  A0)  -  PiR{2}(t) 


P1R(50}( t  +  X 0)  -  P1R{49} (t ) 
P1W1Rl(t  +  Ao)  -  PlR{50/(c) 


Although  this  method  increases  the  number  of  equations  involved 
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in  a  short  simulation,  it  offers  the  advantage  that  waves  produced 
at  the  sare  interface  at  the  same  time  are  simply  added  together, 
rather  than  being  tracked  separately.  As  a  consequence,  the  number 
of  calculations  which  nust  be  performed  is  not  affected  by  the 
production  of  multiple  reflections,  and  the  complexity  of  the  problem 
does  not  increase  with  time.  There  is  a  slight  error  e  introduced 
o y  discretizing  the  transit  times  of  the  layers  (see  Fig.  6),  but  the 
error  car:  be  made  arbitrarily  small  by  reducing  A0  and  increasing  the 
nuauer  of  segments  in  each  layer.  Writing  eqns.  (2-10),  (2-11), 

(2-13)  and  (2-14)  in  this  fashion  completes  the  system  model. 

Appendix  B  contains  a  listing,  a  flow  chart,  and  a  sample 
simulation  session  using  the  BASIC  program  that  has  been  written  to 
perform  the  input-output  simulation  of  a  general  ultrasonic  through- 
transmission  test  system. 


CHAPTER  3 


EXPERIMENTS 


Tn  an  ultrasonic  through-transmission  test  system,  a  specimen 
is  clamped  between  the  faces  of  the  transmitting  and  receiving 
transducers.  An  electrical  waveform  known  as  the  input  signal 
excites  the  piezoelement  in  the  transmitting  transducer,  which  in 
turn  generates  ultrasonic  stress  waves  in  the  specimen.  These 
waves  propagate  through  the  specimen  and  eventually  arrive  at  the 
piezoelement  of  the  receiving  transducer,  which  generates  a 
corresponding  electrical  signal  known  as  the  output  signal.  Both 
the  input  and  output  signals  are  monitored  simultaneously  with  a 
digital  oscilloscope,  and  the  data  stored  for  later  analysis. 

3.1  PROTOTYPE  TRANSDUCERS 

Five  different  transducer  designs  are  evaluated  in  this 
study .  The  Acoustic  Emission  Technology  (AET)  PC-500  broadband 
transducer  serves  as  the  standard  after  which  the  other  four 
designs  are  patterned.  As  shown  in  Fig.  Cl  (see  Appendix  C) ,  the 
PC-50C  consists  of  a  lead  metaniobate  piezoelement  sandwiched 
between  a  protective  wearplate  of  titanium  carbide  and  a  highly- 
attenuating  tungsten  epoxy  backing  layer. 

Three  transducers  have  been  fabricated  for  this  work  which 
are  identical  to  the  PC-500  except  for  the  thickness  of  the 
backing.  The  standard  backing  layer  is  0.635  an  (0.250  in) 
thick,  and  the  prototype  transducers  have  backing  thicknesses  of 


0  an  (no  backing),  1.270  an  (0.500  in),  and  1.905  cm  (0.75P  in). 
These  transducers  are  identified  by  the  numbers  AE-100,  AE-102, 
and  AE-103#  respectively,  and  the  manufacturer  (Panametrics, 

Inc.,  221  Crescent  St.,  Waltham,  MA  02154)  claims  that  they  are 
identical  to  the  standard  AET  PC-500  transducer,  except  for  the 
thickness  of  the  backing. 

In  addition,  a  transducer  has  been  fabricated  which  has  a 
different  wearplate  than  the  standard.  In  order  to  minimize  the 
reflections  generated  at  the  piezoelement-wearplate  and  wearplate- 
specimen  interfaces,  aluminum  was  chosen  as  the  material  for  the 
wearplate  of  this  transducer,  designated  AE-101.  This  aluninun 
wearplate  is  the  same  thickness  as  the  standard  titanium  carbide 
wearplate,  so  that  transducer  AE-101  is  again  identical  to  the 
PC-500,  except  for  the  wearplate  material. 

The  dimensions  and  material  properties  of  all  five  of  these 
prototype  transducer  designs  are  sunmarized  in  Fig.  C2. 

3.2  ELECTRONIC  TEST  EQUIPMENT 

A  pulsed  oscillator  (Arenberg  model  PG-650C)  was  used  to 
generate  the  input  signals  for  the  experiments.  This  signal 
generator  is  capable  of  producing  continuous  and  gated  sinusoids 
ranging  in  frequency  from  0  to  4  MHz,  and  broads nd  pulses  as 
short  as  1  usee  in  duration.  The  peak-to-peak  amplitude  of  these 
waveforms  is  continuously  variable  from  0  to  100  volts. 


The  input  and  output  signals  were  monitored  and  stored  on  a 
Nicolet  model  2090  digital  oscilloscipe,  capable  of  sampling  two 
channels  simultaneously  at  20  MHz  and  with  8-bit  resolution.  The 


digitized  waveforms  were  then  transferred  through  an  RS-232  data 
link  to  an  IBM  Personal  Computer  for  permanent  storage,  analysis 
and  plotting. 

3.3  SPECIMENS  AND  EXPERIMENTAL  PROCEDURES 

Transducers  and  specimens  to  De  tested  were  clamped  together 
at  a  constant  pressure  of  0.4  MPa  (twice  the  "saturation  prtessure" 
of  0.2  MPa  reported  in  (41))  on  a  pneumatic  test  starvd.  A  layer 
of  ultrasonic  couplant  (AET  SC-6)  was  applied  at  all  interfaces 
between  transducers  and  specimens  to  inprove  the  transmission  of 
stress  waves  across  those  interfaces.  An  electrical  attenuator 
was  installed  between  the  signal  generator  and  the  oscilloscope  to 
avoid  overloading  the  oscilloscope.  Fig.  C3  is  a  schematic  of 
the  system  used  for  broadband  pulse  experiments;  Fig.  C4  shows 
the  system  set  up  for  tone  burst  experiments. 

Two  categories  of  experiments  were  performed.  First,  all 
transducers  were  clamped  face-to-face  and  tested  with  broadband 
pulses  (1  psec  duration) ,  and  with  gated  and  continuous  sinusoids 
at  frequencies  of  0.4,  0.6,  1.0,  1.5,  2.0,  2.4,  and  3.0  MHz. 

Next,  some  specimens  were  clamped  between  transducers  and  tested 
with  the  same  input  signals,  to  investigate  the  effect  of  couplant 
layers  and  adhesive  bonds  on  the  propagation  of  stress  waves 
within  the  system.  Fig,  C5  shows  a  sketch  of  the  three  specimens 
tested  in  the  experiments. 

Specimen  A  (see  Fig.  C5)  is  simply  a  1.270  cm  (0.500  in) 
thick  aluminum  plate.  It  was  tested  between  both  AE-103  and 
AE-101  transducers,  to  observe  the  difference  in  output  caused  by 
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impedance-matching  the  wearplone  to  me  specimen.  (The  AE-101  has 
an  aluminum  wearplate,  while  the  AE-103  has  the  standard  titanium 
carbide  wearplate.)  Specimen  5  consists  of  two  0.635  cm  (0.250  in) 
thick  aluninun  plates  separated  oy  a  couplant  layer,  and  in 
specimen  C  the  couplant  layer  is  replaced  by  a  thin  (0.005  an) 
epoxy  layer.  Both  of  these  specimens  were  tested  between  AE-101 
transducers,  and  the  results  compared  to  those  for  the  solid  1.270 
cm  plate  to  determine  the  effect  of  the  intermediate  layers  on  the 
output  of  the  system. 


3.4  EXPERIMENTAL  RESULTS 

Appendix  D  gives  the  measured  output  of  each  of  the  systems 
for  broadband  pulse  and  tone  burst  inputs.  No  output  signals 
corresponding  to  tone  burst  inputs  above  1.0  MHz  are  shown, 
because  the  backing  attenuation  is  so  large  above  this  frequency 
that  the  amplitude  cf  the  back.ng  wave  is  negligible.  The  results 
of  the  continuous-wave  experiments  are  also  given,  in  Table  7. 
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CHAPTER  4 


RESULTS  AND  DISCUSSION 


In  this  chapter,  the  results  of  the  evaluations  of  each  of 
the  transducer  designs  are  presented  and  discussed.  First,  the 
experimental  results  for  each  transducer  are  compared  with  the 
predicted  results  obtained  from  the  wave  propagation  model 
developed  in  Chapter  2,  which  makes  possible  the  quantitative 
evaluation  of  time  delays  and  piezoelectric  constants  in  the 
system.  Next,  the  experimental  results  for  each  of  the  trans¬ 
ducers  are  studied  to  determine  the  role  of  each  transducer 
design  parameter  on  the  output  of  the  system.  The  conclusions 
are  confirmed  with  the  wave  propagation  model,  and  the  model  is 
then  used  to  predict  the  effects  of  additional  parameters  on  the 
output.  The  results  of  experiments  to  determine  the  effects  of 
oouplant  layers  and  adhesive-bonding  layers  are  discussed,  and, 
finally,  a  scheme  is  presented  which  enables  the  model  to  predict 
complete  output  signals  for  broadband  input  signals. 


4.1  RESULTS  OF  CHARACTERIZATION  USING  THE  WAVE  PROPAGATION  MODEL 
The  wave  propagation  model  developed  in  Chapter  2  can  be 
used  t  characterize  an  ultrasonic  test  system.  In  order  to 
obtain  good  correlation  between  the  predicted  results  (from  the 
model)  and  the  measured  results,  it  is  found  that  several  time 
delays  must  be  introduced  into  the  model.  These  time  delays, 
along  with  the  time  constant  RC0  for  exponential  decay  of  the 
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cutput  voltage,  and  the  backing  attenuation  parameter  cx^  (w) ,  make 
up  the  list  of  what  are  called  "modeling  parameters"  in  the 
following  discussion.  These  are  parameters  whose  values  may  be 
arbitrarily  varied  to  optimize  the  correlation  between  measured 
and  predicted  results.  The  characterization  of  an  ultrasonic 
test  system,  then,  consists  of  finding  the  values  of  each  of  the 
modeling  parameters  that  give  the  best  agreement  between  measured 
and  predicted  output  signals. 

Appendix  E  shows  how  the  optimum  values  of  each  of  the 
modeling  parameters  may  be  determined  for  a  particular  ultrasonic 
system.  The  procedure  is  an  iterative  one,  in  which  the  value  of 
a  particular  parameter,  say  4>.,  is  repeatedly  altered  in  the 
simulation  until  the  predicted  output  signal  matches  the  measured 
signal  as  closely  as  possible.  This  value  is  taken  to  be  the 
correct,  or  optimum  value  of  4^  for  the  particular  system  under 
consideration.  The  optimum  value  is  then  used  as  an  input  to  the 
model  while  the  next  parameter  is  varied,  and  so  on  until  all  of 
the  parameters  have  been  optimized.  Table  1  lists  the  values  of 
the  modeling  parameters  for  each  of  the  predicted  curves  in 
Appendix  E,  and  Fig.  7  shows  a  flow  chart  of  the  characterization 
process,  where  represents  the  particular  parameter  being 
optimized  at  any  given  time. 

Two  requirements  must  be  met  in  order  for  this  scheme  to 
work.  First,  the  model  itself  must  be  accurate,  so  that  good 
correlation  between  measured  and  predicted  output  signals  can, 
indeed,  be  obtained  if  the  correct  values  of  the  modeling 
parameters  are  used  as  input  data.  Also,  the  effects  of  the 
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parameters  on  the  output  of  the  system  must  be  independent,  so 
that  changing  one  parameter  does  not  change  the  optimim  value  of 
a  parameter  which  has  already  been  optimized. 

It  is  easy  to  demonstrate  that  the  model  satisfies  the  first 
requirement.  Figs.  E8,  E13,  and  E15  (see  Appendix  E)  show  the 
excellent  correlation  between  measured  and  predicted  results 
which  is  obtained  for  carefully-chosen  values  of  the  modeling 
parameters.  These  cases  correspond  to  broadband  pulse,  0.4  MHz 
tone  burst,  and  0.6  MHz  tone  burst  input  signals,  respectively, 
for  the  system  of  two  PC-500  transducers  clamped  face-to-face. 

The  independence  of  the  modeling  parameters  is  not  so  easily 
established,  but  the  discussion  of  Appendix  E  shows  that,  if  the 
parameters  are  optimized  in  the  proper  order  (5,  2,  RCo,  a,  y, 
h~C0/S) ,  the  effect  of  changing  one  parameter  on  other 
parameters  which  have  already  been  optimized  is  minimal. 

Four  systems  are  characterized  using  this  procedure. 

Each  of  the  systems  consists  of  a  pair  of  transducers  clamped 
face-to-face,  the  four  systems  corresponding  to  the  PC-500, 
AE-101,  AE-102,  and  AE-103  transducers.  The  results  of  the 
characterizations  are  shown  along  with  the  experimental  results 
for  each  system  in  Tables  2-5,*  Table  5  contains  a  sumary  of  the 
optimum  values  of  the  modeling  parameters  for  each  system.  The 
effects  of  each  of  the  parameters  on  the  output  of  the  system 
are  discussed  in  the  following  paragraphs. 


4.1.1  Wearplate-Piezoelement  Time  Delay,  6 

As  mentioned  earlier,  several  time  delays  must  be  included 
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in  the  wave  propagation  model  in  order  to  obtain  good  correlation 
between  measured  and  predicted  results  for  a  given  system.  The 
particular  time  delay  discussed  here,  called  £,  is  introduced  as 
waves  pass  through  the  interfaces  between  che  weatplai.es  and 
piezoelements  in  the  system.  The  effect  of  tnis  time  delay  is 
simply  to  increase  the  time  interval  between  the  application  of 
the  input  signal  and  the  first  appearance  of  the  output  signal, 
as  discussed  in  Appendix  E  (subsection  E.1.2).  This  time  delay 
does  not  affect  the  shape  or  amplitude  of  the  output  signal  in 
any  way. 

The  values  of  6  in  Table  6  are  all  within  a.  01  tsec  of  each 
other,  except  for  transducer  AE-101,  which  is  the  only  transducer 
with  an  aluninun  wearplate.  It  appears  that  a  softer  wearplate 
material  like  aluninum  reduces  the  value  of  the  time  delay. 
Otherwise,  the  value  of  £  seems  to  be  fairly  constant  for  trans¬ 
ducers  of  similar  construction  (same  wearplate  and  piezoelement) . 
In  addition,  no  variation  of  £  with  frequency  is  observed  over 
the  range  from  0. 4-3.0  MHz. 

4.1.2  Wearplate-Piezoelement  Time  Delay,  3 

This  time  delay  is  introduced  as  waves  are  reflected  back 
into  the  wearplates  from  the  wearplate-piezoelement  interfaces 
in  the  system.  Fig.  E3  shows  schematically  the  propagation  of 
stress  waves  through  the  system,  vrfiere  the  solid  black  lines 
represent  the  wavefronts.  Observe  that  the  waves  bounce  back  and 
forth  between  the  wearplate-piezoelement  interfaces  several  times, 
and  that  this  results  in  an  echo  in  the  output  signal.  The  effect 
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of  8  is  to  increase  the  time  between  the  echoes,  as  described  in 
Appendix  E  (subsection  E.1.3).  Predicted  output  signals  corres¬ 
ponding  to  several  different  values  of  8  are  shown  in  Fig.  8 
(reproduced  from  Fig.  E6);  these  should  be  ccn pared  with  the 
measured  broadband  pulse  output  signals  for  each  of  the  systems, 
shown  in  Fig.  9. 

The  value  of  S  varies  considerably  from  transducer  to 
transducer;  the  PC-500  shows  a  value  of  0.24  psec,  while  8  for  the 
AE-101  is  only  slightly  greater  than  half  of  that  (0.125  usee). 
Note  that  the  37.5%  reduction  in  8  associated  with  the  aluminum 
wearplate  is  exactly  the  same  as  the  reduction  in  £ .  Clearly, 
both  of  the  time  delays  at  the  wearplate-piezoelement  interface 
depend  strongly  on  the  wearplate  impedance. 

Observe  from  Table  6  that  both  the  AE-102  and  the  AE-103  have 
identical  values  for  8.  This  is  significant,  because  these  are 
the  only  two  transducers  with  matched  backing  layers  and  titaniun 
wearplates  that  were  made  at  the  same  time.  Thus  the  time  delay 
8  seems  to  be  consistent  from  one  transducer  to  another,  if  the 
transducers  have  identical  wearplates  and  piezoelements  and  were 
manufactured  under  the  same  conditions.  The  discrepancy  between 
the  AE-102  and  AE-103  versus  the  PC-500  is  probably  attributable 
to  differences  is  materials  or  manufacturing. 

The  optimization  of  3  was  initially  carried  out  using  only 
the  broadband  pulse  input  signal.  However,  the  optimum  values 
obtained  for  the  broadband  pulse  were  also  found  to  be  the  optimum 
values  for  tone  burst  input  signals  over  the  entire  frequency 
range  of  the  experiments,  indicating  no  dependence  of  8  upon  the 


frequency  of  the  input  signal 


4.1.3  The  Time  Constant,  RC0 

The  parameter  RC0  can  have  a  very  significant  effect  on  the 
output  signal  for  a  particular  ultrasonic  test  system,  as 
described  in  Appendix  E  (subsection  E.1.4).  RC0  is  the  time 
constant  of  the  series  R-C  circuit  made  up  of  the  output  piezo- 
eleinent  and  the  oscilloscope.  The  value  of  RC0  determines  how 
fast  the  output  voltage  decays  from  the  piezoelement;  if  RC0  is 
smaLl  (on  the  order  of  0.1  usee) ,  the  output  signal  is  propor¬ 
tional  to  the  stress  inside  the  piezoelement,  and  if  R C0  is  large 
(5.0  usee) ,  the  signal  is  proportional  to  the  integral  of  stress. 
This  effect  is  discussed  in  more  detail  in  Chapter  2  (section 
2.5);  Fig.  E7  shows  predicted  output  signals  corresponding  to 
several  different  values  of  RCq. 

Transducers  AE-101,  AE-102,  and  AE-103  all  show  a  time 
constant  of  approximately  5.0  usee,  indicating  that  the  output 
voltage  is  virtually  the  exact  integral  of  the  stress  inside  the 
piezoelement.  The  output  for  the  PC-500,  however,  shovs  a 
slight  decay  with  time,  indicating  a  value  for  the  time  constant 
of  about  1.0  usee.  This  difference  is  almost  certainly  a  result 
of  variations  in  transducer  materials;  the  lead  metaniobate  used 
for  the  PC-500  piezoelement  apparently  has  a  lower  relative 
dielectric  constant  than  that  used  in  the  new  prototype 
transducers,  resulting  in  a  lower  value  for  the  capacitance  C;. 
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4.1.4  Backing  Wave  Time  Delay,  y 

The  function  of  the  backing  layer  in  an  ultrasonic  trans¬ 
ducer  is  to  broaden  the  bandwidth  of  the  transducer  by 
suppressing  the  development  of  standing  waves  in  the  piezoelement. 
There  is  an  additional  effect  of  the  backing  layer,  however. 

Fig.  10  shows  the  broadband  pulse  response  for  each  of  the 
prototype  transducers  that  has  a  matched  backing  layer;  observe 
that  waves  travel  through  the  backing  layer  and  eventually  return 
to  the  piezoelement  at  a  reduced  amplitude,  producing  a 
contribution  to  the  output  voltage  at  that  time.  Hie  PC-500  and 
AE-101  transducers  have  0.635  cm  (0.25 0  in)  backing  layers  with  a 
transit  time  of  3.7  sec,  so  that  the  backing  waves  for  these 
transducers  are  delayed  by  at  least  7.4  sec.  The  AE-1R2  has 
twice  the  backing  thickness  of  the  AE-101  and  the  AE-103  three 
times  the  thickness,  giving  arrival  times  for  the  backing  waves 
of  at  least  14.6  and  22  sec,  respectively.  The  measured  arrival 
times  exceed  the  predicted  times  in  every  case,  indicating  the 
introduction  of  a  time  delay,  y,  as  waves  travel  across  the 
backing  layer.  Fig.  Ell  shows  four  predicted  output  signals 
corresponding  to  four  different  values  of  the  time  delay,  y. 

Table  6  shows  values  for  y  ranging  from  0.75  '^sec  for 
transducer  AE-101  to  1.60  Ksec  for  the  FC- 500.  Again  the  PC-500 
shows  substantially  different  behavior  than  the  others,  but  note 
that  the  ratio  of  the  y  values  for  the  AE-101  and  AE-102  is 
identical  to  the  ratio  of  their  backing  thicknesses.  This 
correlation  with  backing  thickness  would  indicate  that,  rather 
than  being  an  interfacial  effect,  the  time  delay  has  something  to 
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do  with  wave  propagation  through  the  backing  layer  itself. 

The  highest  frequency  at  which  a  backing  wave  can  be  detected 
is  1.0  MHz ,  and  this  only  for  the  transducers  with  0.635  an 
(0.250  in)  backings.  Over  this  limited  range  (0. 4-1.0  MHz),  no 
frequency-dependence  is  observed  for  the  time  delay  v. 

4.1.5  Backing  Attenuation,  ot_  (  u) 

D 

The  effect  of  the  attenuation  of  the  backing  material  is  to 
reduce  the  amplitude  of  the  backing  wave  in  the  output  signal  for 
the  system.  Fig.  E12  shows  four  predicted  output  signals 
corresponding  to  four  different  values  of  backing  attenuation,  afi. 

Fig.  11  shows  a  plot  of  the  measured  attenuation  of  the 
backing  material  (taken  from  Fig.  B6) ,  and  the  values  of  a fi(uj) 
obtained  from  the  characterization  of  the  transducers  are  also 
shown.  Note  that  the  predicted  values  of  a_  are  slightly  greater 

D 

than  the  measured  values,  indicating  that  the  measured  backing 
waves  are  actually  smaller  in  amplitude  than  one  would  expect, 
based  on  the  values  of  backing  attenuation  reported  in  (24). 

The  probable  explanation  of  this  slight  discrepancy  is  that 
there  is  some  other  source  of  attenuation  in  the  system  in 
addition  to  the  backing  layer  itself.  This  would  mean  that  the 
stress  waves  entering  the  backing  layer  would  already  have  been 
attenuated  somewhat.  Perhaps  the  source  of  this  attenuation  is 
the  conversion  of  stress  to  voltage  in  the  piezoelement  of  the 
output  transducer,  rt  seems  reasonable  to  expect  that  the 
production  of  electrical  energy  at  the  surfaces  of  the  piezo¬ 
element  might  be  accompanied  by  a  corresponding  reduction  in 
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the  mechanical  strain  energy  inside  the  piezoelement.  An 
experiment  is  reconmended  in  the  next  chapter  to  confirm  this 
hypothesis. 

2 

4.1.6  The  Amplitude  Factor,  h  C0/S 

The  predictions  from  the  wave  superposition  model  do  not 

2 

include  the  amplitude  factor  h  Co/S,  vrtiere  h,  Cq  and  S  represent 
the  piezoelectric  constant,  capacitance,  and  surface  area  of  the 
piezoelement,  respectively.  Both  h  and  CQ  are  difficult  to 
measure,  because  they  tend  to  vary  significantly  even  from  one 

"identical"  piezoelement  to  another  [10].  In  this  analysis, 

2 

therefore,  the  amplitude  factor  h  CQ/S  is  found  by  dividing  the 
measured  output  amplitude  by  the  predicted  output  amplitude, 
resulting  in  the  value  of  the  amplitude  factor  that  gives  a 

predicted  output  signal  exactly  equal  in  amplitude  to  the 

2 

measured  signal.  The  values  for  h  CD/S  found  in  this  manner  are 
given  in  Tables  2  through  5. 

2 

Fig.  12  shows  a  plot  of  h  (^,/S  versus  frequency  for  four  of 
the  transducer  systems  tested  in  this  stxiy.  (The  correlation 
between  the  predicted  and  measured  output  signals  for  the 
unbacked  transducer  (AE-100)  was  not  good  enough  to  allow  the 
calculation  of  the  amplitude  factor.)  These  curves  correspond 
to  the  steady-state  output  amplitude  produced  by  a  continuous 
sinusoidal  input  signal.  All  of  the  curves  fall  fairly  close 
to  one  another  over  the  range  0. 4-2.0  MHz;  the  average  values  of 
the  amplitude  factor  over  this  range  for  each  of  the  transducers 
are  reported  in  Table  6. 


Observe  from  Fig.  12  that  transducers  AE-102  and  AE-103 
again  show  qualitatively  very  similar  behavior,  with  h2Co/S  for 
the  AE-102  exceeding  that  for  the  AE-103  by  a  fairly  constant 
amount  over  the  entire  frequency  range.  This  seems  to  indicate 
that,  all  other  factors  being  equal,  the  sensitivity  of  a 
transducer  (which  is  proportional  to  the  amplitude  factor) 
decreases  for  increasing  backing  thickness.  Referring  back  to 
the  transient  results  for  the  broadband  pulse  input  in  Fig.  10 
shows  that  they  confirm  this  observation,  with  the  AE-102 
(1.270  cm  (0.500  in)  backing)  showing  almost  a  25%  greater  output 
amplitude  than  the  AE-103  (1.905  cm  (0.750  in)  backing). 

This  correlation  of  output  amplitude  with  backing  thickness 
is  surprising.  At  relatively  low  frequencies,  where  the  backing 
wave  makes  a  significant  contribution  to  the  output,  it  is 
understandable  that  a  transducer  with  a  thin  backing  layer  might 
produce  a  larger  steady-state  output  signal,  because  of  the 
contribution  of  the  backing  wave  to  the  output.  This  is  not  the 
case  here,  however,  for  the  difference  in  amplitude  between  the 
AE-102  and  AE-103  transducers  is  approximately  constant  over  the 
entire  frequency  range  of  the  experiments.  In  addition,  the  pulse 
response  shows  the  same  behavior,  and  the  backing  wave  does  not 
enter  into  the  first  peak  of  the  pulse  response.  This  effect 
cannot  be  adequately  explained  without  further  research. 

The  AE-101  transducer  shows  a  sharp  rise  in  amplitude  factor 
above  2.0  MHz.  This  is  probably  due  to  a  change  in  the  transit 
time  of  the  piezoelement  which  is  caused  by  the  presence  of  the 
aluninun  wearplate.  This  effect  is  discussed  in  detail  in 
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Appendix  F  (section  F.2). 

4.1.7  Characterization  of  the  Unbacked  Transducer 

The  differences  between  the  behavior  of  backed  and  unbacked 
transducers  are  so  great  that  a  separate  discussion  is  required 
for  the  unbacked  transducer  (AE-10S)  tested  in  this  study.  The 
difference  is  actually  the  difference  between  a  damped  system  and 
a  resonant  system,  corresponding  to  the  backed  and  unbacked 
transducers,  respectively.  Some  of  the  particular  complexities 
of  a  resonant  system  are  alluded  to  in  the  discussion  of  steady- 
state  output  amplitude  in  Appendix  F  (expecially  section  F.2). 

In  an  unbacked  transducer,  waves  are  reflected  from  the  rear 
surface  of  the  piezoelement  rather  than  being  transmitted  across 
it,  so  that  there  are  two  more  possible  tine  delays  to  evaluate 
than  for  the  backed  transducer.  These  two  time  delays  correspond 
to  the  reflection  of  waves  back  into  the  piezoelement  from  each 
of  its  faces.  Also,  since  no  attenuating  materials  are  present 
in  the  unbacked  system  except  the  slightly  attenuating  piezo¬ 
element,  the  only  loss  of  stress-wave  energy  from  the  system  is  by 
transmission  into  the  potting  material  from  the  rear  surface  of 
the  piezoelement  and  by  conversion  into  electrical  energy  by  the 
piezoelement  itself.  In  a  backed  transducer,  these  effects  are 
insignificant  compared  to  the  attenuation  of  the  backing  layer, 
but  when  no  backing  is  present,  they  become  very  important.  In 
fact,  in  the  AE-100,  some  waves  undergo  as  many  as  30  reflections 
before  their  amplitudes  are  small  enough  to  be  neglected.  This 
means  that  the  time  delays,  transit  times,  impedances,  and  all 
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other  parameters  must  be  evaluated  much  more  precisely  for  an 
unbacked  transducer  than  for  a  transducer  with  a  backing  layer,  or 
errors  will  accumulate  and  eventually  dominate  the  predicted 
results.  As  the  frequency  cf  the  input  signal  nears  the  resonant 
frequency  o I  the  piezoelement,.  this  problem  becomes  increasingly 
more  severe. 

Fig.  13  shows  four  predicted  output  signals  for  the  AE-100 
transducer,  corresponding  to  one  broadband  pulse  and  three  tone 
burst  input  signals.  These  predictions  are  the  best  predictions 
that  can  be  obtained  without  changing  the  model  to  include  the 
time  delays  and  attenuation  factors  which  are  peculiar  to  the 
unbacked  case.  Fig.  13  should  be  compared  with  to  Fig.  14  (taken 
from  Fig.  D3) ,  which  oontains  the  corresponding  measured  output 
signals  for  the  AE-100.  Because  of  the  enormous  number  of 
reflections  in  the  unbacked  transducer,  it  is  difficult  to  explain 
the  differences  between  the  measured  ard  predicted  results  in 
terms  of  particular  transducer  parameters,  but  some  observations 
can  be  made  nonetheless. 

Observe  for  the  pulse  response  (plot  (a)  in  Figs.  13  and  14) 
that  the  measured  output  shows  appreciably  more  damping  than  the 
predicted  output.  The  measured  signal  only  rings  for  about  10 
usee,  compared  to  a  30  usee  duration  for  the  predicted  output 
signal.  Plot  (d)  illustrates  this  highly-resonant  behavior  even 
more  dramatically;  as  the  input  frequency  nears  the  resonant 
frequency  of  the  piezoelement,  the  predicted  output  amplitude 
increases  by  six  orders  of  magnitude  compared  to  the  measured 
output.  Plots  (b)  and  (c)  in  Figs.  13  and  14  show  slight 
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qualitative  similarity  for  measured  and  predicted  results,  but 
correlation  is  rather  poor  without  including  the  two  time  delays 
mentioned  earlier. 

4.2  EFFECTS  OF  OTHER  PARAMETERS  ON  TOE  OUTPUT  OF  TOE  SYSTEM 
During  the  characterization  of  an  ultrasonic  system,  it 
becomes  necessary  to  investigate  the  effects  of  factors  other  than 
merely  the  time  delays  and  time  constants  discussed  in  the 
previous  section.  These  other  factors  are  not  modeling  parameters 
in  the  sense  that  their  values  cannot  be  arbitrarily  varied  to 
improve  the  accuracy  of  the  predicted  results,  but  they  certainly 
affect  the  output  of  the  system  and  it  is  helpful  to  understand 
their  role  in  the  behavior  of  the  system  as  a  whole.  Each  of 
these  factors  is  discussed  in  detail  in  Appendices  F  through  J, 
and  a  summary  of  their  effects  on  the  output  signal  is  presented 
in  the  following  paragraphs. 

4.2.1  Effect  of  Wearplate  Impedance  and  Backing  Impedance 
Appendices  F  and  G  discuss  the  influence  of  these  two 
quantities  on  the  output  of  an  ultrasonic  test  system;  Appendix  F 
draws  conclusions  based  on  the  measured  steady-state  output 
amplitude,  and  Appendix  G  corroborates  the  conclusions  with 
theoretical  predictions.  The  results  show  that  impedance  mismatch 
between  wearplate  and  piezoelement  reduces  the  efficiency  of  a 
transducer;  maximum  transmission  of  stress-wave  energy  across  the 
wearplate-piezoelement  interface  requires  that  the  two  layers  be 
impedance-ma  tched . 
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The  effect  of  impedance  mismatch  at  the  piezoelement-backing 
interface  is  more  dramatic,  if  the  backing  layer  is  not 
impedance-matched  to  the  piezoelement,  standing  waves  can  develop 
in  the  piezoelement  and  cause  it  to  resonate.  A  piezoelement: 
which  is  free  to  resonate  is  much  more  sensitive  to  excitement 
or  near  its  resonant  frequency  than  at  other  frequencies,  so 
that  its  bandwidth  is  relatively  narrow  compared  to  that  of  a 
piezoelement  with  a  matched  backing  layer. 

4.2.3  Effect  of  Backing  Thickness 

Appendix  F  also  discusses  the  effect  of  backing  thickness  on 
steady-state  output  amplitude  of  the  system.  The  result  is  the 
interesting  and  confounding  observation  that  the  steady-state 
output  amplitude,  and  consequently  the  sensitivity  of  a  trans¬ 
ducer,  decreases  with  increasing  backing  thickness.  This 
observation  was  also  made  in  the  previous  section  during  the 
discussion  of  the  amplitude  factor  h  Q,/S,  where  it  was  pointed 
out  that  even  in  the  transient  case,  where  a  broadband  pulse  input 
signal  is  used,  the  amplitude  of  the  output  signal  for  the  AE-103 
transducer  (1.905  cm  (0.750  in)  backing  thickness)  is  25%  less 
than  that  for  the  AE-102  (1.270  cm  (0.500  in)  backing).  This 
effect  cannot  be  explained  except  by  some  peculiar  behavior  of 
stress  waves  at  the  piezoelement-backing  interface.  Further 
investigation  is  required  to  confirm  and  explain  this  result. 

4.2.4  Effect  of  Piezoelement  Thickness 

In  addition  to  considering  the  effects  of  wearplate  impedance 
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and  backing  Impedance  on  the  output  of  an  ultrasonic  test  system. 
Appendix  G  also  uses  the  wave  propagation  model  to  predict  the 
effect  of  changing  the  piezoelement  thickness.  Increasing  the 
thickness  of  the  piezoelement  is  shown  to  increase  the  anplitude 
of  the  output  signal,  but  at  the  cost  of  losing  some  bandwidth. 
Since  the  piezoelement  integrates  the  stress  inside  it  to  produce 
the  output  voltage,  a  thicker  piezoelement  integrates  over  a 
longer  time  interval  and  consequently  filters  out  some  of  the 
higher- frequency  components  of  the  stress  waveform.  The  best 
piezoelement  is  that  which  has  a  short  enough  transit  time  so  as 
not  to  filter  useful  information  out  of  the  output  signal,  yet  is 
thick  enough  to  give  a  reasonable  output  amplitude. 

There  is  another  effect  of  piezoelement  thickness  which  is 
important  in  resonant  transducers.  The  resonant  frequency  of  a 
piezoelement  is  inversely  proportional  to  its  transit  time,  so 
increasing  the  thickness  of  the  piezoelement  will  lower  its 
resonant  frequency,  and  vice-versa. 

4.2.4  Effect  of  Intermediate  Couplant  and  Bonding  Layers 

Appendix  H  describes  the  experiments  used  to  isolate  the 
effects  of  couplant  layers  and  adhesive-bonding  layers  between 
transducer  components  and  at  transducer-specimen  interfaces  in  the 
system.  The  results  show  that  the  behavior  of  stress  waves  at 
impedance-mismatched  interfaces  where  a  couplant  layer  is  present 
is  exactly  as  predicted  by  the  impedance  model  of  Chapter  2  (see 
subsection  2.3.1),  and  that  only  a  very  small  reflection  (less 
than  2%  of  incident  amplitude)  is  produced  at  an  impedance-matched 
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intet face  where  couplant  is  used.  No  time  delays  are  observed  at 
the  couplant  layer  in  either  case,  although  other  researchers  have 
reported  them  at  higher  frequencies  [42] . 

The  presence  of  an  intermediate  epoxy  bonding  layer  is  shown 
to  have  two  significant  and  easily-quantified  effects:  reduction 
of  the  amplitude  of  the  transmitted  wave,  and  production  of  large- 
amplitude  reflections  at  the  layer.  Both  are  due  to  impedance 
nisnatch  between  the  epoxy  and  the  adjacent  layers. 

4.2.5  Effect  of  Frequency-Dependent  Attenuation  on  Output 
for  Broadband  Input  Signal 

Since  the  attenuation  of  the  backing  material  used  in  the 
transducers  is  frequency-dependent,  the  wave  superposition  model 
does  not  give  an  accurate  prediction  of  the  backing  wave  when  a 
broadband  input  signal  is  used.  One  way  to  solve  this  problem  is 
to  use  a  digital  signal  processing  program  such  as  the  Fast 
Fourier  Transform  (FFT)  to  transform  the  input  signal  into  the 
frequency  domain,  and  then  to  selectively  attenuate  different 
frequency  components  as  required  by  the  frequency-dependent 
attenuation  parameter  aD(J  [24], 

A  simpler  procedure  involves  the  use  of  a  digital  low-pass 


filter  to  simulate  the  increase  in  attenuation  at  higher 
frequencies.  Appendix  I  discusses  this  procedure,  which  gives 
satisfactory  results  for  the  broadband  pulse  used  in  this  study. 


CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  CONCLUSIONS 

Five  different  transducer  designs  were  evaluated  in  this 
study.  One  of  the  designs  is  sold  commercially,  the  Acoustic 
Emission  Technology  (AET)  FC-500  transducer,  manufactured  by 
panametrics,  Inc.  of  Waltham,  MA  (model  V105) .  This  transducer 
consists  of  a  lead  metaniobate  piezoelement  which  converts  stress 
to  voltage  and  vice-versa,  a  titanium  carbide  wearplate  which 
protects  the  piezoelement  from  direct  contact  with  the  specimen 
during  testing,  and  a  tungsten  epoxy  backing  layer  which  broadens 
the  bandwidth  of  the  transducer  by  damping  cut  standing  waves  in 
the  piezoelement.  The  other  prototype  transducer  designs  are 
similar  to  the  PC-500,  differing  only  in  backing  thickness  and 
wearplate  impedance. 

The  transducers  were  evaluated  by  comparing  the  results  of 
through- transmission  experiments  to  theoretical  results  predicted 
by  a  mathematical  model  of  a  general  ultrasonic  test  system.  The 
model  was  developed  specifically  for  this  study,  and  is  based  on 
the  theory  of  wave  propagation  in  elastic  materials  and  on  the 
physics  of  the  piezoelectric  effect.  Once  the  accuracy  of  the 
model  was  verified,  it  was  then  used  as  an  aid  in  understanding 
the  effects  of  certain  transducer  parameters  on  the  output  of  an 
ultrasonic  system.  The  conclusions  which  can  be  drawn  from  this 
study  are  summarized  below. 
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Items  (1)  through  (8)  are  the  results  of  the  comparison  of 
measured  and  predicted  results: 

(1)  In  order  to  obtain  good  correlation  between  measured  and 
predicted  results,  it  was  found  necessary  to  account  for 
three  time  delays  in  the  model.  One  delay,  6 ,  is 
introduced  as  waves  are  transmitted  across  the  interface 
between  the  wearplate  and  piezoelement  of  the  output 
transducer;  another  delay,  6,  is  introduced  as  waves  are 
reflected  back  into  the  wearplate  from  the  same 
interface;  and  a  third  delay,  y  ,  is  introduced  as  waves 
propagate  through  the  backing  layer  and  are  reflected 
from  the  interface  between  the  backing  layer  and  the 
potting  material.  All  of  these  time  delays  are  in 
addition  to  the  time  delays  (transit  times)  assiciated 
with  the  transducer  elements  themselves. 

(2)  The  values  for  the  time  delay  6  were  found  to  be  0.08 
nsec  for  those  transducers  with  titanium  carbide 
wearplates,  and  C.05  nsec  for  those  with  aluminum 
wear plates. 

(3)  The  values  for  the  time  delay  6  were  found  to  be  0.20 
nsec  for  those  transducers  with  titanium,  carbide 
wearplates,  and  0.125  nsec  for  those  with  aluminum 
wearplates. 

(4)  Tlie  values  for  the  time  delay  Y  were  found  to  be  0.75 
nsec  for  the  AE-101  transducer,  with  0.635  cm  (0.250  in) 
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backing  thickness,  and  1.50  usee  for  the  AE-102,  with 
1.270  cm  (0.500  in)  backing  thickness.  The  PC-500  has 
a  backing  thickness  of  0.635  cm  (0.25C  in) ,  and  showed 
a  time  delay  of  1.60  usee. 

(5)  No  frequency-dependence  was  observed  for  any  of  the  time 
delays  over  the  frequency  range  0. 4-3.0  .MHz. 

(6)  It  was  found  that  the  decay  of  voltage  across  the  output 
piezoelanent  due  to  the  flow  of  current  through  the 
oscilloscope  played  an  important  role  in  determining  the 
shape  of  the  output  signal.  The  value  of  the  time 
constant  RC0  corresponding  to  this  decay  was  found  to  be 
l.O  usee  for  the  standard  AET  FC-500  transducer,  and  5.0 
ysec  for  the  prototype  transducers  with  matched  backing 

layers  (AE-101,  AE-102,  AE-1C3) . 

2 

(7)  The  amplitude  factor  h  Cc/S,  where  h,  CD  and  S  represent 
the  piezoelectric  constant,  capacitance  and  surface  area 
of  the  piezoelement,  respectively,  was  found  to  be 
approximately  constant  between  0.4  and  2.0  MHz,  falling 
between  3.5  and  7.0  x  1015  N/m3  for  each  of  the  trans¬ 
ducers  over  that  range. 

(3)  Excellent  correlation  between  predicted  results  from  the 
wave  propagation  model  and  measured  results  was  found 
for  all  of  the  transducers  except  the  AE-100,  provided 
the  correct  time  delays  and  time  constant  were  known. 

Items  (9)  through  (11)  are  further  observations  and  conclusions 


based  on  these  results: 


-62- 


(9)  Replacing  the  titanium  wearplate  with  an  aluninum 
wearplate  reduced  both  of  the  wearplate-piezoelement 
time  delays  by  exactly  the  same  proportion — 37.5%.  This 
shows  clearly  that  these  time  delays  are  functions  of  the 
wearplate  impedance;  a  decrease  in  wearplate  impedance 
results  in  a  decrease  in  the  values  of  the  time  delays 
at  the  interface  between  wearplate  and  piezoelement.  This 
also  causes  an  increase  in  the  resonant  frequency  of  the 
piezoelement  (see  item  (13),  below). 

(10)  The  ratio  between  the  backing-wave  time  delays  for  the 
AE-101  and  AE-102  transduce  s  is  exactly  the  same  as  the 
ratio  of  the  backing  thicknesses  for  these  transducers. 

This  indicates  that  perhaps  this  time  delay  does  not  rep¬ 
resent  an  interfacial  effect,  but  rather  a  correction  for 
errors  in  the  measured  values  for  wavespeed  and/or 
thickness  of  the  backing  layers. 

(11)  Using  the  measured  values  of  backing  attenuation  in  the 
model  gave  predicted  backing  waves  that  were  too  large  by 
a  factor  of  two.  Either  the  measurements  were  inaccurate, 
or  some  other  source  of  backing-wave  attenuation  exists, 
in  addition  to  wave  propagation  through  the  backing  layer. 

Items  (12)  through  (17)  are  statements  concerning  the  influence 
of  certain  transducer  parameters  and  conditions  on  the  output  of  an 
ultrasonic  test  system. 


(12)  The  presence  of  a  layer  of  ultrasonic  acoustic  couplant 
at  interfaces  within  the  system  has  a  negligible  effect 
on  the  propagation  of  stress  waves  within  the  system 
over  the  frequency  range  from  0.4  to  3.0  MHz.  Neither 
the  amplitude  nor  the  phase  of  the  stress  waves  is 
affected  by  the  presence  of  the  couplant  layer. 

(13)  A  mismatch  in  acoustic  impedance  between  the  wearplate 
and  piezoelement  in  a  transducer  results  in  a  decrease 
in  output  amplitude  and  the  production  of  spurious 
reflected  waves  at  those  interfaces  in  the  system. 

An  impedance-mismatched  wearplate  can  also  affect  the 
resonant  frequency  of  the  piezoelement  by  altering  the 
time  delays  at  the  piezoelement-wearplate  interfaces. 

(14)  Impedance  mismatch  between  wearplate  and  specimen  also 
reduces  the  amplitude  of  the  output  signal  and  creates 
spurious  reflections  in  the  output  signal. 

(15)  Impedance  mismatch  at  the  rear  surface  of  the  piezo¬ 
element  in  an  ultrasonic  transducer  allows  the  presence 
of  a  large-amplitude  resonant  peak  in  the  frequency 
response  of  the  transducer. 

(16)  Increasing  the  thickness  of  the  backing  layer  in  a  trans¬ 
ducer  reduces  the  sensitivity  of  the  transducer. 

(17)  Increasing  the  thickness  of  the  piezoelement  in  a  trans¬ 
ducer  increases  the  amplitude  of  the  output  signal,  but 
reduces  the  bandwidth  of  the  signal. 
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5 . 2  RECOMMENDATIONS 

Based  on  the  conclusions  of  the  preceding  section,  several 
sets  of  recoirmendations  can  be  made.  These  are  sumnarized  in  the 
paragraphs  which  follow. 

5.2.1  Recoirmendations  for  Future  Work 

This  study  raised  several  questions  which  deserve  further 
experimentation  and  analysis: 

(1)  Further  investigation  is  required  to  experimentally 
establish  the  existence  of  time  delays  and  amplitude 
reductions  at  the  surfaces  of  the  piezoelement  in  an 
ultrasonic  transducer.  A  specimen  for  this  purpose  could 
be  made  by  bonding  a  piezoelement  between  two  layers  of  a 
material  with  low  attenuation,  such  as  titanium  carbide. 
The  entire  specimen  could  then  be  tested  between  two 
transducers,  and  the  output  signal  analyzed  for  evidence 
of  time  delays  and  amplitude  reduction.  At  the  same 
time,  the  effect  of  changing  the  electrical  load  across 
the  piezoelement  on  these  parameters  could  be  determined 
as  well. 

(2)  A  statistical  study  should  be  performed  on  a  number  of 
allegedly  "identical"  transducers  to  determine  the 
correspondence  between  the  values  of  the  time  delays, 
capacitances,  and  piezoelectric  constants  for  these 
transducers.  This  would  help  to  establish  whether  the 
values  of  these  parameters  could  be  predicted  in  advance 
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for  a  particular  transducer,  or  whether  the  manufacturing 
process  introduces  so  many  variables  that  each  transducer 
must  be  evaluated  individually, 

(3)  An  analysis  similar  to  that  employed  here  should  be 

applied  to  resonant  transducers,  such  as  the  AE-1CC  used 
in  this  study.  The  increased  precision  required  in  the 
resonant  case  would  give  additional  insight  into  the 
behavior  of  ultrasonic  transducers  in  general. 

5.2.2  Recommendations  for  Transducer  Design 

Based  on  the  results  of  this  study,  the  following  guidelines 
are  suggested  for  the  design  of  broadband  piezoelectric  ultrasonic 
transducers: 

(1)  A  broadband  ultrasonic  transducer  should  be  designed  with 
a  backing  layer  that  has  the  same  acoustic  impedance  as 
the  piezoelement.  If  possible,  the  round-trip  transit 
time  of  the  backing  layer  should  be  larger  than  the 
expected  duration  of  the  output  signal.  This  eliminates 
the  possibility  of  interference  between  the  backing  wave 
and  the  primary  wave  in  the  output  signal. 

(2)  The  thickness  of  the  piezoelement  should  be  chosen  so 
that  the  round-trip  transit  time  of  the  piezoelement  is 
equal  to  the  desired  time-resolution  vf  the  output  signal. 
This  gives  the  best  compromise  between  amplitude  and 
bandwidth  in  the  output  signal. 

(3)  The  impedance  of  the  wearplate  should  also  be  chosen  to 
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match  either  the  piezoelement  or  the  specimen.  This 
gives  the  largest-amplitude  output  signal,  and  also 
reduces  the  production  of  spurious  reflected  waves  in  the 
system.  In  addition,  the  wearplate  should  be  designed  as 
thin  as  is  practical,  so  that  its  thickness-mode  resonant 
frequency  is  well  above  the  upper  limit  of  the  bandwidth 
of  the  transducer . 

5.2.3  Recommendations  for  Reporting  of  Ultrasonic  Test  Results 

One  of  the  implications  of  this  study  is  that  there  are  many 
factors  which  influence  the  output  of  a  particular  ultrasonic  test 
system.  When  results  of  ultrasonic  tests  are  reported,  the 
information  given  is  often  incomplete  and  vague,  so  that  it  is 
difficult  for  the  reader  to  determine  exactly  what  has  been  done, 
and  even  more  difficult  to  reproduce  the  work  himself.  To  avoid 
this  problem,  the  following  information  should  be  included  when 
reporting  the  results  of  ultrasonic  tests: 

(1)  The  acoustic  impedance,  wavespeed,  and  dimensions  of  the 
test  specimen  should  be  reported  with  the  test  results. 

(2)  A  description  of  the  transducer  used  should  be  included, 
along  with  either  a  frequency  response  curve  or  a  plot  of 
the  impulse  response  of  the  transducer. 

(3)  The  method  of  attaching  the  transducer  to  the  specimen 
should  be  described,  particularly  if  an  adhesive  is  used. 

(4)  A  description  of  the  device  used  to  measure  the  output 
voltage  should  be  given. 
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Fig.  5  Effect  of  a  finite  time  constant  RC0  on  the  voltage 

produced  by  a  step  function  of  stress  in  the  piezoelement 
of  an  ultrasonic  transducer. 


Discretization  of  the  piezoelement  of  an  ultrasonic  transducer 
into  n  segments,  each  having  transit  time  A  . 


Flowchart  of  the  procedure  for  characterizing  an 
ultrasonic  test  system  using  the  wave  propagation 
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Fig.  8  Predicted  broadband  pulse  response  of  two 

FC-500  transducers  clamped  face-to-face  as  a 
function  of  the  time  delay  B. 
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Fig.  9  Measured  broadband  pulse  response  of  (a)  FC-500, 

(b)  AE-101 ,  (c)  AE-102 ,  and  (d)  AE-103  transducers 
clamped  face-to-face  (expanded  time  scale). 
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Fig.  10  Measured  broadband  pulse  response  of  (a)  FC-500, 

(b)  AE-101 ,  (c)  AE-102,  and  (d)  AE-103  transducers 
clamped  face-to-face. 
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11  Comparison  of  measured  backing  attenuation  with 
values  of  attenuation  obtained  from  the 
characterization. 
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Table  1  Values  of  the  modeling  parameters  corresponding  to  each 
predicted  output  signal. 


the  testing  and  characterization  of  the  AK-101  transducer 


Results  of  the  testing  and  characterization  of  the  AE-102  transducer. 
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APPENDIX  A 


BEHAVIOR  OF  ONE- DIMENSIONAL  PLANE  COMPRESSIONAL  WAVES 
IN  PIEZOELECTRIC  AND  NON PIEZOELECTRIC  MATERIALS 


This  analysis  is  patterned  after  that  of  Redwood  [19] .  One¬ 
dimensional  plane-wave  propagation  in  both  piezoelectric  and  non¬ 
piezoelectric  materials  is  investigated;  the  equations  relating 
stress  and  voltage  in  piezoelectric  materials  are  also  developed. 


A.l  WAVE  PROPAGATION  IN  NONPIEZOELECTRIC  MATERIALS 

In  this  section,  the  one-dimensional,  plane  compressional 
wave  equation  is  derived  for  nonpiezoelectric  materials. 
Consider  a  uniform  elastic  bar  extending  along  the  x-axis  and 
propagating  plane  compressional  waves  as  shown  in  Fig.  Al.  The 
momentum  equation  for  a  differential  element  of  the  bar  may  be 
written  as  follows. 


(T  *  £ d*)  -  t  }  ■  j=s 


where  S  is  the  cross-sectional  area  of  the  bar  in  nr ,  p  is  its 

3  o 

mass  density  in  kg/m  ,  and  T  and  u  are  the  stress  in  N/m"  and 

displacement  in  m,  respectively,  in  the  x-direction.  This  equation 
simplifies  to 


In  addition,  the  relationship  between  stress  and  displacement  in 
the  bar  may  be  written  by  Hooke's  law: 

PREVIOUS  PAGE 
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T 


(A3) 

where  E  is  Voung's  modulus.  Substituting  eqn.  (A3)  into  eqn.  (A2) 
then  yields  the  one-dimensional  wave  equation  for  plane 
compression?.':  wave  propagation  in  a  nonpiezoelectric  material: 


(A4) 


Note  that  the  coefficient  E/o  corresponds  to  the  square  of  the 
wavespeed,  c,  in  the  bar: 


2 

c 


(A5) 


A. 2  WAVE  PROPAGATION  IN  PIEZOELECTRIC  MATERIALS 

In  this  section,  the  equation  for  wave  propagation  in 
piezoelectric  materials  is  derived.  The  momentum  equation  (A2) 
is  unaffected  by  the  piezoelectric  effect,  but  the  Hooke's  law 
relationship  requires  an  additional  term: 


T  =  E  -  hD 

ox 


(A<5) 


This  equation  shows  that  the  presence  of  an  electric  flux  D  in 
the  x-direction  produces  a  so-called  "piezoelectric  stress"  hD, 
where  h  is  the  piezoelectric  constant  of  the  material  in  v/m. 
Substituting  eqn.  (A6)  into  the  momentun  equation  (A2)  yields 


o  u 
.2 


E  p"  U 
2 


h  5D 
D  3x 


of  ~  ox”  "  (A7) 

This  equation  can  be  simplified  by  obtaining  an  alternative 


expression  for  D.  If  there  is  no  free  charge  in  the  material, 


Gauss*  law  may  be  written  as  follows: 


3D  3D  3D 
_JS  +  _I  +  _5.  .  o 

3x  3y  3z 

where  Dx,  Dy,  and  Dz  are  the  electric  flux  densities  in  the  x-, 
y-,  and  z-directions,  respectively.  However,  the  plane-wave 
assunption  (that  plane  waves  remain  plane)  requires  that  the  x- 
and  y-  dimensions  of  the  specimen  be  much  larger  than  the  z- 
dimension,  so  that  all  derivatives  with  respect  to  x  and  y  must 
vanish.  Eqn.  (AS)  then  becomes 


3D 

?x 


0 


(A9) 


and  eqn.  (A7)  reduces  to 

e  u  m  E  3  u 

,  2  "  c  .  2 

at  dx  (A1C) 

which  is  the  one-dimensional  plane  ccnpressional  wave  equation  for 
piezoelectric  materials.  Note  that  this  equation  is  identical  to 
eqn.  (A4),  which  holds  for  nonpiezoelectric  materials. 


A. 3  SOLUTION  TO  THE  WAVE  EQUATION 

In  this  section,  expressions  are  obtained  for  stress  and 
displacement  which  satisfy  the  wave  equation  (A10) .  operational 
calculus  is  used  to  simplify  the  mathematics.  Recall  that  the 
Laplace  transform  is  defined  by 

r  » 

{ g ( t ) }  ■  g'(p)  =  j  e  Pt  g(t)  dt 

>o  (All) 

where  g'(p)  represents  the  Laplace  transform  of  the  time-domain 
function  g(t),  and  p  is  the  Laplace  transform  variable.  For  the 
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remainder  of  this  analysis,  Laplace  transforms  of  quantities  will 
be  referred  to  by  the  names  of  the  quantities  thenselves;  fot 
example,  u'  will  be  referred  to  as  displacement,  and  so  on.  It 
will  be  understood  from  the  superscript  *  that  it  is  actually  the 
Laplace  transform  of  the  variable  that  is  being  considered,  rather 
than  the  variable  itself. 

If  the  value  of  the  function  f(t)  is  0  at  t=P,  it  may  be  easily 
shown  [40]  that 


{T1] 


pg' (p) 


(A12) 


and  consequently  that 


r 

d2g(c) 

9 

dt" 

v.  y 


*  P2g’(p) 


(A13) 


Applying  this  operation  to  eqn.  (A10)  yields  the  transformed 
one-dimensional  plane-wave  equation: 


2  F  /  - 

3x  fc/-  (A14) 

where  u*  is  the  Laplace  transform  of  u(t). 

Hie  solution  to  eqn.  (A14)  can  be  written  as  the  sun  of  two 
exponentials: 


u'  =  A  exp  (-^)+  B  exp 


(A15) 


where  c  is  the  wavespeed  from  eqn.  (A5) ,  and  p  is  the  Laplace 
transform  variable.  The  first  term  of  u'  represents  a  wave 
traveling  in  the  positive  x-direction,  and  the  second  a  wave  in 
the  negative  x-direction;  A  and  B  are  amplitude  factors  that  are 


<3 


evaluated  fron  the  boundary  conditions.  Substituting  eqn.  (A15) 
into  eqn.  (A3)  gives  the  corresponding  solution  in  terms  of 
stress  ./avet . 


T’ 


+  B  exp  (^r)J 


(Alfi) 


2 

where  T'  is  stress  in  the  x-direction. .  Substituting  pc  =E  gives 
the  following: 


5c  p  <  -A  exp 


\  c  , 


+  B  exp 


(f) 


(A17) 


The  coefficient  pc  is  seen  to  be  the  ratio  of  stress  to  particle 
velocity  in  the  material,  and  is  caAied  the  "characteristic 
impedance",  Z: 


Z  =  PC 

Eqn.  (A17)  may  now  be  written  in  terms  of  Z: 


(A1B) 


This  is  the  equation  for  stress  in  the  x-direction  of  a 
nonpiezoelectric  material  which  satisfies  the  one-dimensional, 
plane  compress ional  wave  equation. 


A. A  kAVE  INTERACTION  AT  IMPEDANCE-, MISMATCHED  INTERFACES 

In  this  section,  the  amplitudes  of  the  waves  which  are 
produced  when  a  stress  wave  strikes  an  interface  between  two 
materials  of  different  characteristic  impedance  are  found. 
Consider  the  interface  t  .-tween  two  layers  of  nonpiezoelectric 
material  with  impedances  Z^  and  Z^,  as  shown  in  Fig.  A2.  The 
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right-traveling  displacement  wave  with  amplitude  A^  strikes  the 
interface,  and  produces  a  reflected  wave,  B  ,  and  a  transmitted 
wave,  *  The  boundary  conditions  require  that  displacement  and 
.stress  be  continuous  across  the  boundary: 


u 


1 


X 


=  0 


0 


T  =  T  '  I 

*1  '  l2  1 

x  =  0  X  =  0  (A20) 

where  the  suDScripts  1  anti  2  refer  to  layers  1  and  2,  respectively. 
Substituting  eqn.  (A15)  into  the  displacement  boundary  condition 
gives 

Ai  +  3i  =  A2  +  B: 


(A21) 


and  substituting  ecn.  ( A 1 9 )  into  the  stress  condition  gives 


pZi("Ai  +  E:)  =  pZ  2  (-A2  +  B->) 


(A22) 


If  the  analysis  is  restricted  to  values  of  time  less  than  that 
required  for  wave  A ,  to  reach  the  far  end  of  layer  2,  wave  will 
never  show  up  and  it  may  be  struck  from  eqns.  (A21)  and  (A22) . 
Then,  from  eqn.  (A22) : 


l 

'  U  ('Al  +  Bl) 


and  from.  eqn.  (A21)  : 


Bi  =  A2  -  A1 


(A23) 


(A24) 


Substituting  eqn.  (A24)  into  eqn.  (A23)  gives 
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A 


2 


which  simplifies  to 


22, 


(A25) 


(A26) 


This  equation  gives  the  amplitude  of  the  reflected  displacement 
wave  in  terms  of  the  amplitude  of  the  incident,  wave  and  the 
impedances  of  the  two  layers.  The  amplitude  of  the  transmitted 
wave,  3^  ,  can  also  be  foind: 


3 

A 


1 

1 


(A27) 


Eqns.  (A26)  and  (A27)  may  be  written  in  terms  of  stress  by 
comparing  eqns.  (A15)  and  (A19) .  These  show  that  the  relationship 
between  stress  T*  and  displacement  u'  is  expressed  by 


T'  =  pZu ' 

for  a  right-traveling  wave,  and 


(A2B) 


T'  =  -pZu ' 


(A29) 


for  a  left-traveling  wave.  If  T| ,  T^  and  T^  represent  the 
amplitudes  of  the  incident,  reflected,  and  transmitted  stress  waves, 
respectively,  then  these  amplitudes  are  related  to  the  corresponding 
displacement  amplitudes  ,  and  A,,  by: 


T' 

i 


PVl 


t;  =  -pZli 


T' 

tr 


pZ/2 


(A3J) 


— - : - — — - 
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Making  the  appropriate  substitutions  of  eqns.  (A30)  into  eqns. 
(A26)  and  (A27)  gives  the  following: 


T' 

r 

T1 


Z2~Z1 
Z1  +  22 


TI 


2Z, 


Zl  +  Z2 


(A31) 


These  equations  give  the  amplitudes  of  the  reflected  (T^)  and 
transmitted  (T£r)  stress  waves  which  are  produced  when  a  stress 
wave  of  amplitude  traveling  in  a  median  with  impedance  Z^ 
strikes  another  median  with  impedance  Z2.  Although  the  equations 
are  given  in  terms  of  the  Laplace  transforms  of  the  amplitudes, 
they  hold  for  the  time-domain  amplitudes  as  well.  Since  the 
waves  differ  only  by  an  amplitude  factor,  their  time-dependent 
parts  will  be  identical,  and  will  cancel  out  of  the  numerator  and 
denominator  of  eqns.  (A31) . 


A. 5  RELATIONSHIP  BETWEEN  ELECTRICAL  AND  MECHANICAL  QUANTITIES 
IN  A  PIEZOELECTRIC  PLATE 

In  this  section,  the  fundamental  equation  relating 
displacement  and  voltage  in  piezoelectric  materials  is  derived. 
The  piezoelectric  plate  under  consideration,  henceforth  referred 
to  as  the  piezoelement,  is  depicted  in  Fig.  A3. 

To  obtain  a  relation  between  stress  and  electric  flux  in  the 
piezoelement,  substitute  eqn.  (A15)  into  the  Laplace  transform  of 
eqn.  (A6): 

T’  +  hD*  -  E  *  | -A  exp  (=**)  +  B  exp  (^)J 


(A32) 


The  flux  density  D*  in  the  transducer  is  a  constant  function  of 
x#  as  shown  in  eqn.  (A9) ,  and  is  related  to  the  charge  o'  by 
Gauss'  law: 


D' 


(A33) 


where  S  represents  the  surface  area  of  the  plate.  Substituting 
eqn.  (A33)  into  eqn.  (A32)  gives  the  following  general  equation 
for  stress  in  the  piezoelement: 


r  +  iF  *  p2{“a  ePx  {^r) + ® exp  (?)} 

where,  as  in  eqn.  (A18) *  Z»  c. 

To  determine  the  voltage  across  the  piezoelement,  the 
following  fundamental  equation  must  be  used: 


(A34) 


3x 


h  ,  D* 
-n  -z —  +  — 
3x  c 


(A35) 


where  <f>'(x,p)  is  the  local  potential  in  the  piezoelement,  and  c  is 
the  absolute  permittivity  of  the  piezoelectric  material  in  farads/m 
To  solve  for  V* ,  the  total  voltage  across  the  piezoelement,  eqn. 
(A35)  must  be  integrated  along  the  x-axis  from  0  to  L  (see  Fig.  A3) 


v* 


L  3$'<x,p) 

IT 


dx  *  -hu' 


D’L 


o  x  -  0  e  (A36) 

The  second  term  in  eqn.  (A36)  can  be  simplified  by  realizing  that 
Gauss'  law  can  again  be  used  to  determine  the  capacitance  CQ  of  the 
unstressed  piezoelement: 


.  es 
c  -  — 

o  L 


(A37) 


Substituting  eqns.  (A37)  and  (A33)  into  eqn.  (A3fi)  gives  the 


transducer  voltage  in  tents  of  charge 


x  •  L  n, 

V’  -  -hu'l  + 

x  -  0  Co  (A38) 

This  is  the  fundamental  equation  relating  the  quantities  of 
displacement,  voltage  and  charge  in  a  piezoelectric  transducer. 


A. 6  ELECTRICAL  RESPONSE  OF  A  PIEZOELECTRIC  TRANSDUCER 
TO  AN  INCIDENT  STRESS  WAVE 

In  this  section,  a  step  input  of  stress  is  applied  to  a 
piezoelectric,  transducer  and  the  resulting  voltage  is  fowd. 
Consider  a  piezoelement  and  wear pi ate  bonded  together  as  in  Fig. 
A4.  A^  is  the  displacement  associated  with  the  incident  wave  in 
the  wearplate,  and  B ^  and  A  are  the  amplitudes  of  the  reflected 
and  transmitted  displacement  waves,  respectively,  vhich  are 
produced  when  wave  A^  strikes  the  interface  between  the  layers. 
Eqn.  (A20)  gives  the  boundary  conditions  at  x*0: 


x  -  0 


t:  •  T’l 

A  x  ■  0  z  x  -  0 


(A39) 


where  quantities  without  subscripts  refer  to  the  transducer,  and 
those  subscripted  1  refer  to  the  wearplate.  Substituting  eqns. 

(A15)  and  (A34),  the  boundary  conditions  (A39)  become 
Ax  +  Bx  -  A  *  3 

pZl(”Al  +  Bl}  *  PZ<‘A  ♦  B)  -  ^  (A40) 

If  the  analysis  is  restricted  to  values  of  time  less  than  that 
required  for  wove  A  to  reach  the  far  side  of  the  transducer,  wove 
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B  may  be  ignored.  Eqns.  (A40)  then  become 


A1  +  Bi  *  A 


p2l^-Al  +  “  “P2A 


(Ml) 


(A42) 


and  solving  eqn.  (A41)  for  B1  and  substituting  into  eqn.  (A42) 


then  gives 


-2A^pZj  +  A  2^  ■  -A  Z  - 


(A43) 


This  equation  relates  mechanical  quantities  and  A  to  the 
charge  Q1  on  the  piezoelement.  If  a  resistance  R  is  connected 
across  the  piezoelement,  the  voltage  V'  may  be  expressed  as 


Kut  • *'*  ■ 


(AM) 


where  I'  represents  the  current  in  amperes  through  the  resistance  R, 
O'  is  the  charge  on  the  transducer  in  coulombs,  and  p  is  the 
Laplace  transform  variable.  Substituting  eqn.  (A44)  into  ecpi.  (A43) 
and  collecting  terms  yields: 


Ap(Zx  +  Z) 


or,  solving  for  A: 


out 

-j*r  +  V2! 


(A45) 


c  r  .1 

hV'  1 

'  pRS  +  2pZlAl  M  P(Zi  +  Z>  ‘ 

w  V-  / 


(A46) 


This  equation  relates  the  amplitude  of  the  displacement  wave  A 
and  the  voltage  produced  in  the  piezoelement  to  the  amplitude 
of  the  incident  di splacement  wave,  A^. 
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Eqn.  (A38) ,  reproduced  here,  gives  the  voltage  across  the 
plate  in  terms  of  the  displacements  at  its  faces: 

x  -  L 


Vout  -  -hu' 


(A47) 


Substituting  eqns.  (A1S)  and  (A39) ,  the  voltage  becomes 


'out  •  -hf  •**  (=?)  -  1  }  -  55 


out 

pRC 


o  (A48) 

Note  that  iyc  is  the  "transit  time"  of  the  transducer,  the  tine 
required  for  a  wave  to  pass  through  it.  Denoting  this  transit 
time  by  X,  eqn.  (A48)  becomes 


v;ut--h*{«pW-i}-?5 


out 

pRC 


(A49) 


and  substituting  eqn.  (A46)  gives 


V’  „ 
out 


-h 


p(z1  +  Z) 


hV'  _ 
out 


pRS 


+  2pZ1A1  -  jexp  (-p>.^  -1 


V' 

out 

pRC 


(A5C) 


Solving  for  v'  gives  the  final  result: 


out 


(  _  'I 

r  2  l 

,  1  . 1 

1  P  {l  -  exp(-pX)} 

Z1 +  2 1 

1  _e _ h2{.l_-  ex£(-p_>) 

V  J 

LP  RC  RS(Z.  +  Z)  J 

(A51) 


V(Jut  here  is  the  voltage  generated  across  the  piezoelement  when 
a  displacement  wave  of  amplitude  A  enters  it  from  an  adjacent 
wearplate  of  impedance  Z^.  R  is  the  resistance  across  the  piezo¬ 
element,  and  h,  Z,  C  and  S  are  the  piezoelectric  constant, 

o 

characteristic  impedance,  clamped  capacitance  and  surface  area  of 
the  transducer,  respectively. 


Eqn.  (A51)  may  be  cast  in  terms  of  stress,  if  desired,  if 

represents  the  amplitude  of  the  incident  step  function  of 

*•-  •  •  * 

stress  in  the  wearplate,  the  Laplace  transform  of  the  step  [40]  is 
T^/p  ,  so  that  the  stress  in  the  wearplate  is  expressed  by 


7 '  "Mi4i  - 


Solving  for  A  in  terms  of  Ti  gives 


(A52) 


1 

Substituting  this  value  for  A^  into  eqn.  (A51)  results  in 


(A53) 


v  '  .  - di- 

out  { z1  +  Z) 


(1  -  exp(-pl) } 


hZ(l  -  exp(-pA)? 

rs(zl  +  Z) 


(A54) 


which  gives  in  terms  of  the  amplitude  T*  of  a  step  function 
of  stress  in  the  wearplate. 

Examination  of  eqn.  (A54)  shows  that  the  factor  exp(-pA) 
appears  twice.  This  factor  represents  a  time  delay  [40]  equal  to 
A,  so  that  the  terms  vhich  it  multiplies  do  not  show  up  at  the 
sane  time  as  the  rest  of  the  terms  in  eqn.  (A54) .  Since  the 
analysis  was  restricted  to  values  of  time  less  than  the  transit 
time  of  the  piezoelement,  these  time  delay  terms  must  be  struck 
from  the  equation.  Taking  only  those  terms  which  show  up 


inrnedisfely: 


Vout  •  TzTV' Z)  < p2  +  RC~  - 

A  Q 


Jl _ \ 

!i  *  «J 


(A55) 


The  denominator  of  eqn.  (A55)  may  be  written  as  (pfl/t^) (pfl/T  ) * 


The  inverse  Laplace  transform  is  then  easily  found  in  terms  of  the 
time  constants  t1  and- t2 : 


(A57) 


This  equation  gives  the  time-dependent  output  voltage  which 
is  produced  across  the  piezoelement  by  a  step  input  T±  of  stress 
f ran  the  adjacent  wearplate.  The  impedances  of  the  wearplate  and 
piezoelement  are  and  Z,  respectively,  and  h  is  the 
piezoelectric  constant. 


A.7  MECHANICAL  RESPONSE  OF  A  PIEZOELECTRIC  TRANSDUCER 

TO  AN  INPUT  VOLTAGE 

In  this  section,  the  equations  are  developed  which  give  the 
time-dependent  amplitudes  of  the  stress  waves  produced  by  a  step 
input  of  voltage  to  a  piezoelectric  transducer.  Consider  a 
piezoelement  bonded  to  a  wearplate  on  one  face  and  a  backing  layer 
on  the  other  face  as  seen  in  Fig.  A5.  As  shown  in  the  drawing, 
wave  propagates  into  the  wearplate,  wave  propagates  into  the 
backing,  and  waves  A  and  B  are  the  right-  and  left-traveling  waves 
generated  inside  the  piezoelement  itself.  The  impedances  of  the 
wearplate,  backing,  and  piezoelement  are  designated  Z  ,  Z2#  and 
Z,  respectively. 


The  boundary  conditions  at  the  front  face  of  the  piezoelement 
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are  given  in  eqn.  (A41) : 


Bx  -  A  +  B 


pZjB  -  p2(-A  +  B)  -  ^ 


(A58) 


The  boundary  conditions  at  the  back  face  are  obtained  by 


evaluating  eqns.  (A15)  and  (A34)  at  x«L  and  equating  stresses  and 
displacements  in  the  two  layers: 

*  A  exp  (-pX)  +  B  exp  (pX) 


pzj-A  exp(-pX)  +  B  exp  (pX)  j  -  ■  -pZ. 


(A59) 


Here,  as  in  section  A. 6,  X  represents  the  transit  time  of  the 
piezoelement,  L/c.  Combining  eqns.  (A58)  gives 


^  -  -p^Zj^  +  Z)  +  B(ZX  -  Z)| 


(A60) 


and  eqns.  (A59)  combine  to 


^S~  *  p|a(-Z1  +  Z2)  exp  (-pX)  +  B(Z  +  z2)  exp  (pX)j 


(A61) 


In  order  to  reduce  the  amount  of  algebra  in  the  derivation,  the 
following  substitutions  will  be  made: 


m  *  exp(-pX)  ;  n  *  exp(pX) 


(A62J 


so  that  eqn.  (A61)  becomes 

^g“  ■  P  |a  m  (-Z1  +  Z2>  +  Bn(Z  +  Z2>| 

Eliminating  hQ*  from  eqns.  (Afie)  and  (A63)  then  yields 

A  jm(Z  -  Z2)  -  (Z:  +  Z)J  -  B jn(Z  +  Z£)  +  ( Z x  -  Z)J 


(A53) 


(A64) 


‘Ails  gives  the  relationship  between  the  amplitudes  of  the  two 
displacement  waves  produced  inside  the  transducer. 


The  voltage  across  the  transducer  is  introduced  by  recalling 
eqn.  (A38) : 


x  "  L  o’ 

V’  .  -hu’l  + 

x  -  0  o 


(A65) 


Evaluating  eqn,  (A15)  at  x-L  and  substituting  gives 
—  -  -hlAm  +  Bn  -  (A  +  B) ]  + 

P  1  J  Co  (A66) 
Here,  the  input  voltage  has  been  taken  to  be  a  step  of  amplitude 
V  ,  so  that  V**Vo/p.  Multiplying  through  by  hCQ/S  and  collecting 
terns  gives 


h  C  V 

_ 

Sp 


h2  C 


s  -  |\(m  -  1)  +  B(n  -  1)  |  + 


Eliminating  hQ'  from  eqns.  (A6C)  and  (A67)  then  gives 


hC  V 
o  o 

Sp 


r  2 
h  C 


Av 


■g—  (1  -  m)  -  p(Z1  +  Z) 


r  2 

n  C 


+ 


(1  -  n)  +  p(Z  -  Z.) 


Solving  eqn.  (A64)  for  A  and  inserting  into  eqn.  (A68): 

B 


hC 

o  ° 


Sp 


{m(Z  -  Z2)  -  (Zx  +  Z)} 


|n(Z  +  Z2)'+  (Z2  -  Z)}  {-j-2  (1  *  »)  “  P<z!  +  Z)| 
{m(Z  -  Z2)  -  (Z:  +  Z)}{ — g-2-  (1  “  n)  +  p(Z  -  Z^J 


Combining  terms  then  results  in 


(A<V7) 


(A68) 


(A69) 


'VfwiiifSIjf! 
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•i  r. 


hC  V 
o  o 


m(2  -  Z2)  -  (2  +  Z) 


/A,  . 

|  jz^Cn  “  °)  +  2Z(m  +  n  -  2)  +■  Z2(n  -  o)| 

+  pjm(Z  -  Z2)(Z  -  Z2)  -  n(Z  +  Z^^  +  Z)J  . 

> 


(A70) 


or,  solving  for  B: 


,T,  {°«  -  Z2)  -  *>} 


(A71) 


where 


(CD)  -  P{(Z  +  Z2)(Z1  +  Z)  -  (Z  -  Z2)(Z  -  Z1)  exp(-2pX)j 

"  +  Z2)  exp(-2pX)  -  1  -  2z{exp(-2pX)  -  2exp(-pX)  +  l}J  fAj2) 

Eqn.  (A71)  gives  the  amplitude  of  the  left-traveling  displacement 
wave  produced  in  the  transducer  by  a  step  input  of  voltage,  V  , 
Similarly,  solving  eqp.  (A64)  for  B  and  substituting  into  eqn. 
(AS8)  gives  the  amplitude  of  the  right-traveling  wave: 


,v,  {n(Z  .  22)  +  -  »} 


(A73) 


Now,  since  uf(0,p)«A+B: 


hC0v0  [z(m  +  n  -  2)  +  Z2(n  -  m)j 


x  -  0 


(A74) 


Eqn.  (A74)  gives  the  displacement  at  the  front  face  of  the 


transducer  due  to  a  step  input  of  voltage.  It  may  be  expanded 

by  noting  from  eqns.  (A62)  that  mn»l  and  that  £  -exp(-2p  X ) . 

n 

Multiplying  by  m/m  and  making  these  substitutions: 


[•] 


jz{exp(-2pX)  -2exp(-pX)  -  l]  +  Z^il  -  exp(-2p\)}J 


Eqn.  (A75)  is  similar  to  eqn.  (A54)  in  that  several  of  the  terms 
are  delayed  in  time,  some  by  X  and  some  by  2X.  These  delays  arise 
from  the  propagation  of  waves  A  and  B  through  the  piezoelement. 

After  a  delay  of  X,  wave  B  has  reached  the  front  face  of  the 
piezoelement  and  will  affect  the  displacement  there.  After  a 
delay  of  2X,  wave  A. will  have  been  partially  reflected  from  the  rear 
face  and  will  have  returned  to  the  front  face,  where  it  too  will 
affect  the  displacement.  At  this  point,  however,  only  the  waves 
actually  produced  at  the  surfaces  of  the  transducer  by  the  applied 
voltage  are  of  interest.  The  effects  of  waves  propagating  through 
the  system  will  be  considered  later.  Therefore,  taking  only  the 
part  of  eqn.  <A75)  which  show®  no  delay,  and  substituting  for  {CD} 

from  eqn.  (A72)  gives: 

hC  v  f  -  ,  x 

i'|  -  — f-2  (Z  +  z2>  Jp  |(z1  +  Z)(Z  +  z2)j 

x  "  0  I  .2  ^  -1 

Ph  c  ,  i 

•  "T2,  1Z1  +  22  +  Z2)  (A76) 


+  2Z  +  Z„ 


(A76) 


Dividing  by  (Z  +Z) (Z+Z  )  above  and  below: 


„  ph2C  f(Z,  +  Z)  +  (Z  +  Z_) 


^ 


> 


.  » 

'  i 

s 


l 


7 


•  l 


$ 


;  « 


i 

I 
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another,  the  front  face  then  becomes  the  back  face  and  vice-versa. 


Therefore,  interchanging  Z^  and  Z 2  in  eqn.  (A77): 


hC  V 
o  0 

2  ^o! 

'(Z1  +  Z)  +  (Z  +  Z2)  ’ 

i] 

S(Z  +  z2) 

P  -  s  - 

(Z  +  Z)(Z  +  Z  ) 

L  J 

J 

-1 


(A78) 


This  equation  then  gives  the  displacement  at  the  rear  face  of  the 
transducer  caused  by  a  step  input  of  voltage.  Note  that  a  sign 
change  occurs  when  the  surfaces  of  the  transducer  are  reversed, 
because  the  direction  of  +x  is  effectively  reversed  as  well. 

A  time  constant,  x,  may  be  extracted  from  the  denominator 
of  eqns.  (A77)  and  (A79): 


1 

T 


f 


(Z1  +  Z)  +  (Z  + 
"(Z1  +  Z)(Z  +  z2 


so  that  eqns.  (A77)  and  (A78)  become,  respectively: 


hC  V 
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s(z:  +  Z) 


“  P/" 
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S(Z  +  z2) 

~2  ~~ 
P  -  p/x 


(A8B) 


These  equations  show  that  applying  a  voltage  across 
the  piezoelement  produces  a  displacement  at  each  of  its  faces. 
These  displacements  then  cause  four  stress  waves  to  be  propagated, 
one  in  each  direction  from  each  face.  The  amplitudes  of  these 
waves  are  obtained  from  the  stress-displacement  relationships  of 
ecps.  (A28)  and  (A29),  If  T^L’  and  T^R*  represent  the  amplitudes 
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of  the  stress  waves  produced  In  the  wearplate  and  backing, 
respectively,  and  if  TO'  and  TL'  represent  the  right-  and  left¬ 
traveling  waves  in  the  transducer,  then: 


TL^  -  -pZ^  u'  | 


hC  V  Z, 
oo  I  1 


x  -  0 


TR'  -  pZ  a’ 


5  2!  +  Mp2-p/T 


hC  v  . 
o  o  I  Z 


x  ■  0 


TL’  -  -pZ  u’ 


s  Z1  +  Z  p2-p/t 


hC  V  _ 

oo  Z 


S  Z+Z2  P2-P/T 


1^2  *  PZ2U '  i 


hC  V  Z. 
o  o  I  2 


x  -  L 


S  )Z 


+  Z2 J  (p2  -  p/T 


Finally,  taking  the  inverse  transform  of  each: 
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These  equations  give  the  time-dependent  anplitudes  of  the  four 
stress  waves  produced  by  a  step  input  of  voltage  to  a  piezoelectric 
transducer.  The  value  of  the  time  constant  x  is  found  from  eqn. 
(A79) . 
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The  production  of  reflected  and  transmitted  waves  at 
the  interface  between  two  materials  of  different 
characteristic  impedance. 


-116- 


i 

j 


Wearplate 

Piezoelement  Backing 

A1 

A 

-  Bl  — 

B 

X 

\  *  / 
\ — -vW - / 

“  V(T)+ 

-  0  X 

• 

•  L 

Fig.  A4 


i 


Multilayer  piezoelectric  transducer  used  to  detect 
ultrasonic  stress  waves 


< 


» 


-117- 


Wearplate 

Piezoelement  | 

Backing 

t 

A 

' 

.  B1 

4  ■■ 

- - 

A - - /, 


V(T) 

X  *  0  X  -  L 


Fig.  A5  Multilayer  piezoelectric  transducer  used  to 
generate  ultrasonic  stress  waves. 
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ATTACHMENT  B1 


LISTING  OF  THE  BASIC  MODELING  PROGRAM 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 
39* 
40: 


"DIM  B1R(85 0) ,B1L(850) ,B2R{850) ,B2L(850) 

DIM  P1R(275) ,P1L(275) ,P2R(275) ,P2L(275) 

DIM  W1R(75)  ,W1L(75)  ,W2R(75)  ,W2L(75) 

DIM  SIR (250) ,S1L(250) ,S2R(125) ,S2L(125) 

DIM  S3R(125) ,S3L(125) ,VIN(4000) ,VOUT<4000) 

'  INPUT  DATA 

CP-274000 

ZP-1800000 

CB-173000 

ZU-160000 

'  INTERACTIVE  INPUT 

1  PRINT  "ENTER  INPUT  DATA  (PRESS  RETURN  FOR  DEFAULT  VALUE) 1 
PRINT  *************************************************** 
PRINT  "  " 

PRINT  “DATA  FOR  BACKING  LAVERS:" 

PRINT  " - " 


",LB1 

",LB2 


INPUT"  BACKING  THICKNESS  OF  INPUT  TRANSDUCER  IN  CM  (DEF-0) : 

INPUT" BACKING  THICKNESS  OF  OUTPUT  TRANSDUCER  IN  CM  (DEF-0) : 

INPUT" IMPEDANCE  OF  BACKING  LAVERS  IN  GM/CM 2/SEC  (DEF-1.8E06) :  ",ZB 
IF  ZB-0  THEN  ZB-1.8E06 

INPUT" ATTENUATION  OF  BACKING  LAVERS,  ALPHA  (DEF-1.6):  ", ALPHA 
IF  ALPHA-0  THEN  ALPHA-1.6 

INPUT" FREQUENCY  FILTERING  PARAMETER,  DELTA  (DEF-0) :  ", DELTA 
PRINT  "  * 

PRINT  "DATA  FOR  WEARPLATES:" 

PRINT  " - " 

INPUT"WEARPLATE  THICKNESS  OF  TRANSDUCERS  IN  CM  (DEF-0. 035):  ",LW 
IF  IN-0  THEN  LW-0.035 

INPUT"WEARPLATE  WAVES PEED  OF  INPUT  IN  CM/SEC  (DEF-8.64E05) :  ",CW1 
IF  041-0  THEN  O41-8.64E05 

INPUT*WEARPLATE  WAVES PEED  OF  OUTPUT  IN  CM/SEC  (DEF-8.64E05) :  ",CW2 
IF  CW2-0  THEN  CW2-864000 

INPUT" WEARPLATE  IMPEDANCE  OF  INPUT  IN  GM/O! 2/SEC  (DEF-4.5E06) :  ",ZW1 
IF  ZW1-0  THEN  ZW1-4500000 

INPUT-WEARPLATE  IMPEDANCE  OF  OUTPUT  IN  GM/CM2/SEC  (DEF-4.5E06) :  ",ZW2 
IF  2W2-0  THEN  ZW2-4500000 
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41: 

42: 

43: 

44: 

4S: 

46: 

47: 

48: 

49: 

50: 

51: 

52: 

53: 

54: 

55: 

56: 

57: 

58: 

59: 

60: 

61: 

62: 

63: 

64: 

65: 

66: 

67: 

68: 

69: 

70: 

71: 

72: 

73: 

74: 

75: 

76: 

77: 

78: 

79: 

80: 

81: 

82: 

83: 

84: 

85: 

86: 

87: 

88: 

89: 

90: 


PRINT  "  " 

PRINT  "DATA  FOR  PIEZOELEMENTS:" 

PRINT  " - " 

INPUT" THICKNESS  OF  PIEZOELEMENTS  IN  CM  (DEF-0. 073):  ",LP 
IF  LP-0  THEN  LP-0. 073 

INPUT"TIJ«  CONSTANT  OF  INPUT  TRANSDUCER,  TAU  (DEF-LARGE) :  ",TAU 
IF  TAU-0  THEN  TAU-1E10 

INPUT" TIME  CONSTANT  OF  OUTPUT  TRANSDUCER,  RCO  (DEF-1E-06) :  ",RCO 
IF  RQ>0  THEN  RCO-1E-06 
PRINT  "  " 

PRINT  "DATA  FOR  SPECIMEN (S) :" 

PRINT  " - " 

INPUT  "ARE  THERE  ANY  SPECIMENS  PRESENT?  [Y/N] :  ", SPECS 

IF  SPBC$»"N"  OR  SPEC$*"n"  THEN  10 

INPUT" HOW  MANY  SPECIMEN  LAYERS?  [1,2,3):  ",NSPEC% 

INPUT* THICKNESS  OF  LAYER  1  IN  CM  (DEF-0.635):  ",LS1 
IF  LSI-0  THEN  LSI-0.635 

INPUT" WAVES PEED  OF  LAYER  1  IN  CM/SEC  (DEF-4.81E05) :  ",CS1 
IF  CS1-0  THEN  CS1-4.81E06 

INPUT" IMPEDANCE  OF  LAYER  1  IN  GM/CM2/SEC  (DEF-1.30E06) :  ",ZS1 
IF  ZS1-0  THEN  ZS1-1.303E06 

INPUT" ATTENUATION  OF  LAYER  1  IN  NEPERS/CM  (DEF-0) :  " , ALPHA1 
IF  NSPBC%-1  THEN  10 

INPUT" THICKNESS  OF  LAYER  2  IN  CM  (DEF-0) :  ",LS2 
INPUT" WAVESPEED  OF  LAYER  2  IN  CM/SEC  (DEF-0):  ",CS2 
INPUT" IMPEDANCE  OF  LAYER  2  IN  GM/CM2/SEC  (DEF-0):  ",ZS2 
INPUT  ATTENUATION  OF  LAYER  2  IN  NEPERS/CM  (DEF-0):  ",ALPHA2 
IF  NSPEC4-2  THEN  10 

INPUT* THICKNESS  OF  LAYER  3  IN  CM  (DEF-0):  ",LS3 
INPUT" WAVESPEED  OF  LAYER  3  IN  OV'SBC  (DEF-0) :  "  ,CS3 
INPUT" IMPEDANCE  OF  LAYER  3  IN  GM/CM2/SEC  (DEF-0):  ",ZS3 
INPUT" ATTENUATI CN  OF  LAYER  3  IN  NEPERS/CM  (DEF-0):  " ,ALPHA3 
10  PRINT  "  " 

PRINT  "DATA  FOR  INPUT  WAVEFORM:" 

PRINT  " - " 

INPUT" INPUT  AUTOMATICALLY  OR  FROM  FILE?  [A/F] ;  ",INP$ 

IF  INP$-"A*  OR  INP$-"a"  THEN  20 

INPUT" NAME  OF  INPUT  DATA  FILE:  ",NIMP$  :  GOTO  30 

20  INPUT" INPUT  WAVEFORM  (S,G,P,I)  (NO  DEFAULT):  ",ANS$ 

IF  ANS$«"I"  OR  ANS$-"i"  THEN  30 
IF  ANS$*"P"  OR  ANS$-"p"  THEN  FREQ-0.50  :  GOTO  30 
INPUT* FREQUENCY  OF  INPUT  IN  MHZ  (DEF-1.0):  ",FREQ 
IF  FREQ-0  THEN  FREQ-1.0 
30  PRINT  "  " 

PRINT  "NAMES  OF  OUTPUT  FILES:" 

PRINT  " - " 

INPUT" OUTPUT  FILE  NAME  FOR  INPUT  WAVEFORM:  ",IN$ 

INPUT" OUTPUT  FILE  NAME  FOR  OUTPUT  WAVEFORM:  ",OUT$ 

PRINT  "  " 

PRINT  "NUMERICAL  DATA  FOR  SIMULATION:" 


91:  PRINT  " - " 

92:  40  INPUT- SIZE  OF  TIME  STEP  IN  SBC  (DEF-1.333E-08) :  ", DELTAT 

93:  IF  DELTAT-0  THEN  DELTAT*1.33333E-08 

94:  INPUT" FINAL  TIME  IN  MICROSEC  (DEF-50) :  ",TMAX 

95:  IF  TMAX=0  THEN  TMAX*50 

96:  NEND%»CINT(TMAX/DELTAT/1000000) 

97:  IF  NEND%>4000  THEN  NEND»4000 
98: 

99:  '  PRELIMINARY  CALCULATIONS 
100: 

101:  NP%«CINT(LP/(DELTAT*CP))-1 

102:  IF  LP*0  THEN  »>%«1 

103 :  NW1%CINT  (LW/  ( DELTAT* CW1 ) )  -1 

104:  IF  LW*0  THEN  NWlt-1 

105 :  NW2%*CINT ( LW/ (DELTAT*CW2 ) ) -1 

106:  IF  LW»0  THEN  NW2%»1 

107:  NB1%*CINT(LB1/ (DELTAT*Cb) ) -1 

108:  IF  LB1-0  THEN  NB1«>0 

109 :  NB2%-CINT ( LB2/ ( DELTAT*  CB) ) -1 

110:  IF  LB2-0  THEN  NB2%«0 

111:  IF  NSPEC%-0  THEN  45 

112:  NS1 %*CINT { LSI/ ( DELTAT*CS1 ) ) -1 

113:  IF  NSPEC%-1  THEN  45 

114:  NS2%CINT{LS2/(DELTAT*CS2))-1 

115:  IF  NSP£C%*2  THEN  45 

116:  NS3%*CINT ( LS3/ { DELTAT*CS3) ) -1 

117:  45  IF  NP%>275  THEN  PRINT  "NP  EXCEEDS  MAXIMUM  ARRAY  SIZE!"  :  GOTO  50 
118:  IF  NW1%>75  THEN  PRINT  "NW1  EXCEEDS  l-AXIMUM  ARRAY  SIZE!"  :  GOTO  50 
119:  IF  NW2%>75  THEN  PRINT  "NW2  EXCEEDS  MAXIMUM,  ARRAY  SIZE!"  :  GOTO  50 
120:  IF  NB1%>850  THEN  PRINT  "NB1  EXCEEDS  l-AXIMUM  ARRAY  SIZE!"  :  GOTO  50 

121:  IF  NB2%>850  THEN  PRINT  "NB2  EXCEEDS  MAXIMUM  ARRAY  SIZE!"  :  GOTO  50 

122:  IF  NS1%>250  THEN  PRINT  "NS1  EXCEEDS  NAXIMUM  ARRAY  SIZE!"  :  GOTO  50 

123:  IF  NS2%>125  THEN  PRINT  "NS2  EXCEEDS  MAXIMUM  ARRAY  SIZE!"  :  GOTO  50 

124:  IF  NS3%>125  THEN  PRINT  "NS3  EXCEEDS  MAXIMUM  ARRAY  SIZE!"  :  GOTO  50 

125:  GOTO  60 

126  :  50  PRINT  "DELTAT  IS  TOO  SMALL.  TRY  A  LARGER  VALUE.  .  ."  :  GOTO  40 
127:  60  NMAX%*NP% 

128:  IF  NW1%>NMAX%  THEN  NMAX%«NW1% 

129:  IF  NW2%>NMAX%  THEN  NMAX%»NW2% 

130:  IF  NB1%>NMAX%  THEN  NMAX%-NB1% 

131:  IF  NB2%>NMAX%  THEN  NMAX%-NB2% 

132:  IF  NS1%>NMAX%  THEN  J*IAX%«NS1% 

133:  IF  NS2%>N!AX%  THEN  NMAX%»J62% 

134:  IF  NS3%>NMAX%  THEN  NMAX%»NS3% 

135:  PRINT  "  " 

136:  PRINT  «********************************************" 

137:  PRIOT  "  " 

138:  PRINT  "NUttEN  OF  ITERATIONS  REQUIRED  (NEND)  -  ",NEND% 

139:  PRIOT  "NUfCER  OF  SEGMENTS  IN  PIEZXLEMENTS  (NP)  -  ",NPt 

140:  PRINT  "NUTCER  OF  SEGMENTS  IN  WEARPLATES  (NW1,NW2)  -  ",NW1%;NW2% 


141:  PRINT  "SEGMENTS  IN  BACKING  LAYERS  (NB1,NB2)  •  ",NB1%;NB2% 

142:  PRINT  "SEGMENTS  IN  SPECIMEN  LAYERS  (NSlfNS2,NS3)  ■  " ,1€1%;NS2%;NS3% 

143:  TIMET-  (NMAX%+1) "DELTAT 

144:  IF  LB1-0  THEN  ERRB1-0  :  GOTO  70 

145:  ERRB1* (ABS ( { (NB1%+1) *DELTAT) - (LB1/CB) ) /TIMET) *100 

146:  70  IF  LB2-0  THEN  ERRB2-0  :  GOTO  80 

147:  ERRB2-  (ABS(  ( (NB2%+1)*DELTAT)-(LB2/CB) ) /TIMES')  *100 

148:  80  IF  IS1-0  THEN  ERRS 1-0  :  GOTO  90 

149 :  ERRS1* (ABS ( ( (NS1%+1 ) *DELTAT) - (LS1/CS1) ) /TIMET) *100 

150:  90  IF  LS2-0  THEN  ERRS2-0  :  GOTO  100 

151 :  ERRS2- (ABS ( ( (NS2%+1 ) *DELTAT) - ( LS2/CS2) ) /TIMET) *100 

152:  100  IF  LS3-0  THEN  ERRS3-0  :  GOTO  110 

153 :  ERRS3- (ABS ( ( (NS3%+1) *DELTAT) - (LS3/CS3) ) /TIMET) *100 

154:  110  ERRW1-(ABS(((NW1%+1)*DELTAT)-(LW/CW1))/TIMET)*100 

155:  ERRW2- (ABS ( ( (NW2%+1) *DELTAT) -(LW/CW2) ) /TIMET) *100 

156:  ERRP»(ABS( ( (NP%+1) *DELTAT)-(LP/CP) ) /TIMET) *100 

157:  PRINT  "  " 

158:  PRINT  ******************************************************** 

159:  PRINT  "  " 

160:  PRINT  "ERRORS  IN  TRANSIT  TIMES  OF  ELEMENTS:  " 

161:  PRINT  " - — " 

162:  PRINT  "BACKING  OF  TRANSMITTER:  ";ERRB1;"%" 

163:  PRINT  "BACKING  OF  RECEIVER:  ";ERRB2;"t" 

164:  PRINT  "WEARPLATE  OF  TRANSMITTER:  ",-ERFWlj "%" 

165:  PRINT  "WEARPLATE  OF  RECEIVER:  ";ERPW2;"%" 

166:  PRINT  "PIEZOELEMENT  OF  EACH:  ";ERRP;"I" 

167:  PRINT  "SPECIMEN  LAYER  1:  ";ERRS1;*%" 

168:  PRINT  "SPECIMEN  LAYER  2:  ";ERRS2;*%" 

169:  PRINT  "SPECIMEN  LAYER  3:  ";ERRS3;"%" 

170:  PRINT  "  " 

171:  INPUT" ARE  THESE  ERRORS  ACCEPTABLE?  [Y/N] :  ",OK$ 

172:  IF  OK$-"Y"  OR  OK$-"y"  THEN  120 

173:  PRINT  "TOY  ANOTHER  TIME  STEP.  .  .  "  :  GOTO  40 

174:  120  PRINT  ■  " 

175:  PRINT  ■****************************************" 

176:  RHOl-EXP ( -ALPHA* LB 1) 

177:  RHO 2-EXP (-ALPHA* LB2) 

178:  SRH01-EXP(-ALPHA1*LS1) 

179:  SRH02«EXP(-ALPHA2*LS2) 

180:  SRHQ3*EXP (-ALPHA3*  LS3) 

181: 

182:  ’  INPUT  WAVEFORM  AS  A  FUNCTION  CF  TIME 
183: 

184:  PRINT  "  " 

185:  IF  INP$»"A"  OR  INP$»"a"  THEN  130 
186:  PRINT  "  *  READING  INPUT  DATA  FILE" 

187:  OPEN  *I",1.NIMP$  s  INPUT|1,NI% 

188:  FOR  I%-1  TO  NI%  :  INPUT#l,VtN (It)  :  NEXT  It  :  CLOSE  #1 
189:  GOTO  150 
190:  130  PRINT  " 


*  CREATING  INPUT  WAVEFORM" 


191:  VZN(1)"100 

192:  IF  ANS$»"I"  OR  ANSS-"!"  THEN  150 
193:  FOR  Ifl  TO  NEH5% 

194:  IF  (ANS$*"P*  OR  ANS$»"p")  AND  I%>-150  THEN  150 
195:  IF  (ANS$*"G"  OR  ANS$-"g")  A*®  I%>-750  THEN  150 
196:  VIN (1%)— 1000*006 (I%*DELTAT*FREQ*6283200) 

197:  IF  (ANS$»"P"  OR  ANS$-"p")  THEN  VIN(I%)«500+VIN(I%)/2 
190:  140  NEXT  1% 

199: 

200:  1  CONVERSION  OF  INPUT  VOLTAGE  TO  STRESS 
201: 

202:  150  PRINT  "  *  CONVERTING  INPUT  VOLTAGE  TO  STRESS" 

203:  VIN (0)*0 

204:  FOR  I%»1  TO  NENDt 

205:  VIN (NEND%+1-I%}*VIN (NEND%+1-I%)-VIN (NEND%-I%) 

206:  NEXT  1% 

207:  FOR  I%«1  TO  N£ND% 

208:  VIN(I%)*VIN (I1)+(EXP (DELTAT/TAU) *VIN (I%-1) ) 

209:  NEXT  1% 

210* 

211:  '  ITERATION  TO  CALCULATE  TIME  RESPONSE 
212: 

213:  PRINT  "  *  WAVE  PROPAGATION  t  COMPLETED" 

214:  FOR  IT%«1  TO  NEND% 

215:  PDQNE%«INT ( IT%/NEND%*100 ) 

216:  LOCATE  23,24  :  PRINT  POONE% 

217: 

218:  '  REFLECTIONS  AT  IMPEDANCE-MISMATCHED  INTERFACES 
219: 

220:  IF  NB1%«0  THEN  160 

221:  B1RO0* ( (ZU-ZB) / (2U+ZB) ) *B1 LI0 

222 :  B1P1RO* ( ( ( 2*ZP) "B1P1RI ) + ( ( ZB-ZP) *B1P1LI ) + (ZP*VIN ( IT% ) ) ) / ( ZB+ZP) 

223:  B1P1L0(((ZP-ZB)*B1P1RI)+<(2*ZB)*B1P1LI)-(ZB*VIN(IT%)))/(ZB+ZP) 

224:  GOTO  170 

225:  160  P1RO0*(( (ZU-ZP)*P1LI0)+(ZP*VIN(IT%)))/(ZU+ZP) 

226:  170  P1W1R0»(((2*ZW1)*P1W1RI)+<(ZP-ZW1)*P1W1LI)-(ZW1*VIN(IT%)))/(ZP+ZW1> 
227 :  P1W1L0-  ( ( ( ZW1-ZP) *P1W1RI ) + ( ( 2*ZP) "P1W1LI )  +  ( ZP*VIN ( IT%) ) ) / ( ZP+ZW1 ) 

228:  IF  NS1%«0  THEN  200 

229 :  W1S1R0- { ( 2*  ZS1 ) *W1S1RI+ { ZW1-ZS1 ) * W1S1 LI ) / ( ZW1+ZS1 ) 

230:  W1S1U>((ZS1-ZW1)*W1S1RI+(2*ZW1)*W1S1LI)/(ZW1+ZS1) 

231:  IF  NS2%>0  THEN  180 

232:  S1W2R0* ( ( 2*ZW2) *S1W2RI+ ( ZS1-ZW2) *S1W2LI ) / ( ZS1+ZW2) 

233 :  S1W2L0- ( ( ZW2-ZS1 ) *S1W2RI+ ( 2*ZS1 ) *S1W2LI )  /( ZS1+ZW2) 

234:  GOTO  210 

235:  180  S1S2R0* ( ( 2*ZS2) *S1S2RI+ ( ZS1-ZS2) "S1S2LI ) / ( ZS1+ZS2) 

236 :  S1S2L0* {(ZS2-ZS1 ) *S1S2RI-M2*ZS1 ) *S1S2LI) /( ZS1+ZS2) 

237:  IF  NS3%>0  THEN  190 

238 :  S2W2RO* ( (2*ZW2) "S2W2RI+ (ZS2-ZW2) *S2W2LI ) / ( ZS2+ZW2) 

239:  S2W2L0> ( (ZW2-ZS2) "S2W2RI+ <2*ZS2) *S2W2LI)/( ZS2+ZW2) 

240:  GOTO  210 
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241:  190  S2S3RO«<(2*ZS3)*S2S3RI+(ZS2-ZS3)*S2S3LI)/(ZS2+ZS3) 
242:  S2S3LO((ZS3-ZS2)*S2S3RI+(2*ZS2)*S2S3LI)/(ZS2+ZS3) 

243:  S3W2RO { (2*ZW2) *S3W2RI+(ZS3-ZW2) *S3W2LI)/{ZS3+ZW2) 

244:  S3W2LO < (ZW2-ZS3) *S3W2RI+ (2*ZS3) *S3W2LI)/(ZS3+ZW2) 

245:  GOTO  210 

246:  200  W1W2R0-((2*ZW2)*W1W2RI+(ZW1-ZW2)*W1W2LI)/(ZW1+ZW2) 
247 :  W1W2LO  { ( ZW2-ZW1 ) *W1W2RI+  <2*ZW1 )  *W1W2LI ) /( ZW1+ZW2) 

248:  210  W2P2RO  ( (2*ZP)  *W2P2RI+  ( ZW2-ZP)  *W2P2LI )  /  ( ZW2+ZP) 

249 :  W2P2L08  ( (ZP-ZW2) *W2P2RI+  <2*ZW2)  *W2P2LI)  /( ZW2+ZP) 

250:  IF  NB2%»0  THEN  220 

251 :  P2B2RO-  ( ( 2*ZB) *P2B2RI+ ( ZP-ZB) *P2B2LI ) / ( ZPfZB) 

252:  P2B2LCN  ( (ZB-ZP) ‘P2B2RI+ (2*ZP) *P2B2LI) /(ZP+ZB) 

253:  B2LO0* ( ( Z0-Z8) / (ZB+ZU))*B2RI0 
254:  GOTO  230 

255:  220  P2LO0*((ZU-ZP)/(ZPfZU))*P2RI0 
256: 

257:  '  DELAY  LINE  EQUATIONS 
258: 

259:  230  FOR  I%*1  TO  NMAX%-1 
260:  IF  I%>-NS1%  THEN  240 
261:  SIR (NS1%+1-I%) *S1R (NS1%-I% ) 

262:  S1L(I%)«S1L(I*+1) 

263:  240  IF  I%>«NS2%  THEN  250 
264:  S2R(NS2%+1-I%)«S2R(NS2%-I%) 

265:  S2L(I%)»S2L(I%+1) 

266:  250  IF  I%>-NS3%  THEN  260 
267:  S3R(NS3%+1-I%)*S3R(NS3%-I%) 

268:  S3L(I%)-S3L(I%fl) 

269:  260  IF  I%>»NP%  THEN  270 
270:  P1R(NP%+1— I%)=P1R(NP%— 1%) 

271:  P1L(I%)*P1L(I%+1) 

272:  P2R ( NP%+1-I% ) “P2R (NP%— 1%) 

273:  P2L(I%)*P2L(I%+1) 

274:  270  IF  I%>»NW1%  THEN  280 
275:  W1R(NW1%+1-I%)-W1R(NW1%-I%) 

276:  W1L(I%)-W1L(I%+1) 

277:  280  IF  I%>*NW2%  THEN  290 
278:  W2R (NW2%+1-I%) *W2R (NW2%-I%) 

279:  W2L(I%)«W2L(I%+1) 

280:  290  IF  I4>*NB1%  THEN  300 
281:  B1R(NB1%+1-I%)«B1R(NB1%-I%) 

282:  B1L(I%)»B1L(I%+1) 

283:  300  IF  I%>«NB2%  THEN  310 
284:  B2R (NB2 %+l-I% ) =B2R {NB2%-I% ) 

285:  B2L(I%)»B2L{I%+1) 

286:  310  NEXT  1% 

287: 

288:  *  BOUNDARIES  OF  DELAY  LINES 
289: 

290:  II  NBlt-0  THEN  320 


O 
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291:  B1R(1)*B1RO0 

292 :  B1LI0® (1-DELTA) * ( B1L ( 1 ) *RH01 ) +DELTA* B1LI 0 

293  :  B1P1RI* ( 1-DELTA) *(B1R(NB1%) *RH01 ) +DELTA* B1 PI RI 

294:  B1 L ( NBl % ) =B1 Pi LO 

295:  P1R(1)=B1P1R0 

296:  B1P1LI=P1L(1) 

297:  GOTO  330 

298:  320  P1R(1)=P1RO0 

299:  P1LI0=P1L(1) 

300:  330  P1W1RI=P1R(NP%) 

301:  P1L(NP%)=P1W1LC 
302:  W1R(1)*P1W1R0 
303:  P1W1LI=W1L(1) 

304:  IF  NS1%=0  THEN  360 
305:  W1S1RI=W1R(NW1%) 

306:  W1L(NW1%)=W1S1L0 
307:  SlR(l)aWlSlRO 
308:  W1S1LI»S1L(1)*SRH01 
309:  IF  NS2%>0  THEN  340 
310:  S1W2RI=S1R(NS1%)*SRH01 
311:  S1L (NS1%) =*S1W2L0 
312:  W2R(1)»S1W2R0 
313:  S1W2LI*W2L(1) 

314:  GOTO  370 

315:  340  S1S2RI»S1R(NS1%)*SRH01 

316:  S1L(NS1%)=S1S2L0 

317:  S2R(1)*=S1S2R0 

318:  S1S2LI=S2L ( 1 ) *SRIC2 

319:  IF  NS3%>0  THEN  350 

320:  S2W2RI=S2R(NS2%)*SRHD2 

321:  S2L(NS2%)*S2W2LO 

322:  W2R(l)*S2W2RO 

323:  S2*'2LI=W2L(1) 

324:  GOTO  370 

325:  350  S2S3RI=S2R(NS2%)*SRH02 
326:  S2L (NS2%) *S2S3LO 
327:  S3R(l)*S2S3RO 
328:  S2S3LI*S3L { 1 ) *SRH03 
329:  S3W2RI=S3R(NS3%) *SRH03 
330:  S3L(NS3%)*S3W2LO 
331:  W2R(l)*S3W2RO 
332:  S3W2LI*W2L(1) 

333:  GOTO  370 

334:  360  W1W2RI=W1R(NW1%) 

335:  W1L(NW1%)*W1W2L0 
336:  W2R ( 1 ) *W1W2R0 
337:  W1W2LI»W2L(1) 

338:  370  W2P2RI»W2R(NW2%) 

339:  W2L(NW2%)  M42P2L0 
340:  P2R(1)«VV2P2R0 


t 
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341:  W'2F2LI-P2L(1) 

342:  IF  NB2%*0  THEN  380 
343:  P2B?RI=P2R(NP%) 

344:  P2L (NP%) *P2B2LO 
345:  B2R { 1 ) =P2B2R0 

346:  P2B2LI= (1-DELTA) * (B2L(1) *RH02) +DELTA*P2B2LI 
347:  B^Ri.3-  (1-DELTA)  *(B2R(NB2%)*RHO2)+DELTA*B2RI0 
348:  B2L(NB2%)-B2LO0 

349 :  VOUT ( 1T%) = (W2P2LI-W2P2RO) + (P2B2RI-P2B2LO) 

350:  GOTO  390 

351:  380  P2RI0=P2R{NP%) 

352:  P2L(NP%)^P2LO0 

353 :  VC'JT  ( IT% , • (w2P2LI-W2P2RO)  +  (P2RI0-P2LO0) 

354:  390  NEXT  IT% 

355: 

356:  '  CONVERSION  OF  STRESS  TO  OUTPUT  VOLTAGE 
357: 

358:  PRINT  "  *  CONVERTING  STRESS  TO  OUTPUT  VOLTAGE" 

359:  FOR  I%=1  TO  NEND% 

360 :  VOUT ( 1% ) = ( 2 *DELTAT* VOUT ( I % ) / ZP)  +  (VOUT( I%-1 ) *EXP (-DELTAT/RCO) ) 
361:  NEXT  J.% 

362: 

363:  '  WRITING  OF  OUTPUT  DATA  FILES 
364: 

365:  PRINT  "  *  WRITING  OUTPUT  DATA  FILES" 

366:  OPEN  "0*,1,IN$  :  WRITE# 1,NEND% 

367:  FOR  1%=1  TO  NEND%  :  WRITE#1 ,VIN<I%)  :  NEXT  1%  :  CLOSE  #1 
368:  OPEN  "0",1,0UT$  :  WRITE#1,NEND% 

369:  FOR  I%=1  TO  NEND%  :  WKITE#l,VOUT(I%)  :  NEXT  1%  :  CLOSE  #1 
370:  PRINT  "  " 

371:  PRINT  "****************************************" 

372:  PRINT  "  " 

373:  INPUT"WANT  ANOTHER  SIMULATION?  [Y/N]  :  ", REDOS 
374:  IF  REDO$="Y"  OR  R£DO$="y"  TH£N  1 
375:  END 


ATTACHMENT  B2 


SAMPLE  SIMULATION  USING  THE  MODELING  PROGRAM 


ENTER  INPUT  DATA  (PRESS  RETURN  FOR  DEFAULT  VALUE) 

************************************************* 

DATA  FOR  BACKING  LAYERS: 


BACKING  THICKNESS  OF  INPUT  TRANSDUCER  IN  CM  (DEEM)) :  .635 
BACKING  THICKNESS  OF  OUTPUT  TRANSDUCER  IN  CM  (DEF=0) :  .635 
IMPEDANCE  OF  BACKING  LAYERS  IN  GM/CM2/SEC  (DEF-1.8E06) : 
ATTENUATION  OF  BACKING  LAYERS,  ALPHA  (DEF*1.6) :  1.8 
FREQUENCY  FILTERING  PARAMETER,  DELTA  (DEF=0) : 

DATA  FOR  WEARPLATES: 


WEARPLATE  THICKNESS  OF  TRANSDUCERS  IN  CM  (DEF=0.C35) : 
WEARPLATE  WAVESPEED  OF  INPUT  IN  OVSEC  (DEF=8.64E05) : 
WEARPLATE  WAVESPEED  OF  OUTPUT  IN  CM/SEC  (DEF-8.64E05) : 
WEARPLATE  IMPEDANCE  OF  INPUT  IN  GM/CM2/SEC  (DEF=4.5E06) : 
WEARPLATE  IMPEDANCE  OF  OUTPUT  IN  GM/CM2/SEC  (DEF=4.5E06) : 

DATA  FOR  PIEZOELEMENTS: 


THICKNESS  OF  PIEZOELEMENTS  IN  CM  (DEF=0.073): 

TIME  CONSTANT  OF  INPUT  TRANSDUCER,  TAU  (DEF= LARGE) : 
TIME  CONSTANT  OF  OUTPUT  TRANSDUCER,  ROD  (DEF=lE-06) : 

DATA  FOR  SPECIMEN (S): 

ARE  THERE  ANY  SPECIMENS  PRESENT?  [Y/N] :  N 
DATA  FOR  INPUT  WAVEFORM: 


INPUT  AUTOMATICALLY  OR  FROM  FILE?  [A/F] :  F 
NAME  OF  INPUT  DATA  FILE:  IN. DAT 

NAMES  OF  OUTPUT  FILES: 


OUTPUT  FILE  NAME  FOR  INPUT  WAVEFORM:  IN .OUT 
OUTPUT  FILE  NAME  FOR  OUTPUT  WAVEFORM:  OUT. OUT 


I 
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NUMERICAL  DATA  FOR  SIMULATION: 


SIZE  OF  TIME  STEP  IN  SEC  (DEF*1.333E-08) : 

FINAL  TIME  IN  MICROSEC  (DEF*50) :  10 

******************************************** 

NUMBER  OF  ITERATIONS  REQUIRED  (NEND)  *  750 

NUMBER  OF  SEGMENTS  IN  PIEZOELEMENTS  (NP)  * 

NUMBER  OF  SEGMENTS  IN  WEARPLATES  (NW1.NW2)  » 

SEGMENTS  IN  BACKING  LAYERS  (NB1 ,NB2)  »  274  274 

SEGMENTS  IN  SPECIMEN  LAYERS  (NS1,NS2,NS3)  * 

****************************************************** 


ERRORS  IN  TRANSIT  TIMES  OF  ELEMENTS: 


j 


* 

■ 


BACKING  OF  TRANSMITTER: 
BACKING  OF  RECEIVER: 
WEAKPLATE  OF  TRANSMITTER: 
WEARPLATE  OF  RECEIVER: 
PIEZOELEMENT  OF  EACH: 
SPECIMEN  LAYER  1: 

SPECIMEN  LAYER  2: 

SPECIMEN  LAYER  3: 


.1053446  % 
.1053446  % 
1.389166E-02  % 
1.389166E-02  % 
6.618141E-03  % 
0  % 

0  % 

0  % 


ARE  THESE  ERRORS  ACCEPTABLE?  [Y/N] :  Y 

a*************************************** 


*  READING  INPUT  DATA  FILE 

*  CONVERTING  INPUT  VOLTAGE  TO  STRESS 

*  WWE  PROPAGATION  100  %  COMPLETED 

*  CONVERTING  STRESS  TO  OUTPUT  VOLTAGE 

*  WRITING  OUTPUT  DATA  FILES 

**************************************** 


19 

2  2 
0  0 


WANT  ANOTHER  SIMULATION  [Y/N]  :  N 


13. 


Wave 

|  Propagation 

Algorithm 


(see  Fig.  B2) 


Fig.  B1  Flowchart  of  BASIC  modeling  program 
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Fig.  Cl  Sketch  of  Acoustic  Emission  Technology  (AET) 
FC-500  broadband  ultrasonic  transducer. 
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Transducer 


Dimensions  and  properties  of  transducer  elements  for  FC-500 
AE-100,  AE-101 ,  AE-102,  and  AE-103  transducers. 


Ultrasonic  test  system  set  up 
broadband  pulse  input  signals 
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Fig.  C5  Sketch  of  specimens  used  to  study  effects  of 
couplant  and  adhesive-bonding  layers  in  an 
ultrasonic  test  system. 
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t  APPENDIX  D 

RESULTS  OF  TRANSDUCER  EXPERIMENTS 


t 


Time  (usee) 


Fig.  D2  Measured  ouput  signals  corresponding  to  (a) 
broadband  pulse,  (b)  0.4  MHz  tone  burst,  (c) 
0.6  MHz  tone  burst,  and  (d)  1.0  MHz  tone  burst 
input  signals  for  two  AE-100  transducers 
clamped  face-to-face. 


Time  (usee) 


Fig.  D3  Measured  output  signals  corresponding  to  (a) 
broadband  pulse,  (b)  0.4  MHz  tone  burst,  (c) 
0.6  MHz  tone  burst,  and  (d)  1.0  MHz  tone  burst 
input  signals  for  two  FC-500  transducers 
clamped  face-to-face. 


Fig.  DA  Measured  output  signals  corresponding  to  (a) 


broadband  pulse,  (b)  0.4  MHz  tone  burst,  (c) 
0.6  MHz  tone  burst,  and  (d)  1.0  MHz  tone  burst 


Time  (usee) 

Fig.  D5  Measured  output  signals  corresponding  to  (a) 
broadband  pulse,  (b)  0.4  MHz  tone  burst,  (c) 
0.6  MHz  tone  burst,  and  (d)  1.0  MHz  tone  burst 
input  signals  for  two  AE-102  transducers 
clamped  face-to-face. 


o  (d)  10  20  30  40 

lime  (usee) 

Fig.  D6  Measured  output  signals  corresponding  to  (a) 
broadband  pulse,  (b)  0.4  MHz  tone  burst,  (c) 
0.6  MHz  tone  burst,  and  (d)  1.0  MHz  tone  burst 
input  signals  for  two  AE-103  transducers 
clamped  face-to-face. 


18  Measured  output  signals  corresponding  to  (a)  0.6  MHz  and  (b)  1.0  MHz  tone 
burst  input  signals  for  1.270  cm  (0.500  in)  aluminum  plate  clamped 
between  two  AE-103  transducers. 


Fig.  D9  Measured  output  signals  corresponding  to  (a)  broadband  pulse  and  (b)  0. 

tone  burst  input  signals  for  1.270  cm  (0.500  in)  aluminum  plate  clamped 


Fig.  i/lO  Measured  output  signals  corresponding  to  la)  0.6  MHz  and  (b>  i.O  MHz 

tone  burst  input  signals  for  1.270  cm  (0.500  in)  aluminum  plate  clamped 
between  two  FC-500  transducers. 


Measured  output  signals  corresponding  to  (a)  broadband  pulse  and  (b)  0.4  MHz 
tone  burst  input  signals  for  1.270  era  (0.500  in)  aluminum  plate  clamped 
between  two  AK-101  transducers. 


Measured  output  signals  corresponding  to  (a)  0.6  MHz  and  (b)  i.O  MHz  tone 
burst  input  signals  for  1.270  cm  (0.500  in)  aluminum  plate  clamped 


between  two  AE-101 


Measured  output  signals  corresponding  to  (a)  0.6  MHz  and  (b)  1.0  MHz  tone 
burst  Input  signals  for  two  0.635  cm  (0.250  In)  aluminum  plates  with 
intermediate  couplant  layer  clamped  between  two  AE-101  transducers. 


Fig.  D15  Measured  output  signals  corresponding  to  (a)  broadband  pulse  and  (b)  0.4  MHz  tone 
burst  input  signals  for  two  0.635  cm  (0.250  in)  aluminum  plates  with  intermediate 
epoxy  layer  clamped  between  two  AE-10)  transducers. 


Measured  output  signals  corresponding  to  (a)  0.6  MHz  and  (b)  1.0  MHz  tone 
burst  input  signals  for  two  0.635  cm  (0.250  in)  aluminum  plates  with  intermediate 
epoxy  layer  clamped  between  two  AE-101  transducers. 
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APPENDIX  E 

CHARACTERIZATION  OF  AN  ULTRASONIC  TEST  SYSTEM  USING 
THE  WAVE  PROPAGATION  MCDEL 

In  this  appendix,  the  wave  propagation  model  developed  in 
Chapter  2  is  used  to  characterize  a  particular  ultrasonic 
through- transmission  test  system,  consisting  of  two  AET  PC-500 
transducers  clamped  face-to-face.  The  characterization  is 
accomplished  by  using  the  wave  propagation  model  to  predict  the 
output  signal  for  a  particular  input  signal,  and  then  comparing 
that  predicted  output  signal  to  the  measured  output.  The  values 
of  various  parameters  which  serve  as  inputs  to  the  model  are  then 
changed  one  at  a  time  to  obtain  the  best  agreement  between 
predicted  and  measured  output  signals.  In  the  following  sections, 
this  process  will  be  referred  to  as  "correlation." 

Fig.  El  shows  the  measured  output  of  the  system  for  a 
broadband  pulse  input.  Observe  that  the  output  consists  of  two 
distinct  wave  packets,  called  the  primary  wave  and  backing  wave, 
respectively.  The  primary  wave  corresponds  to  wave  interaction 
in  the  piezoelements  and  wearplates  of  the  system,  and  the  backing 
wave  is  a  result  of  waves  traveling  through  the  backing  layers 
[24].  Since  the  two  wave  packets  represent  the  influence  of 
different  transducer  components,  they  can  be  correlated  separately. 

The  correlation  of  each  wave  packet  takes  place  in  four 
steps: 

(1)  Experimental  observations  are  made  concerning  the 
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particular  wave  under  consideration. 

(2)  The  wave  propagation  model  is  used  to  generate  a 
preliminary  prediction,  based  on  the  current  best 
estimates  for  all  modeling  parameters. 

(3)  The  values  of  various  parameters  are  selectively  changed 
to  improve  the  correlation  between  predicted  and  measured 
results. 

(4)  A  "best"  prediction  is  made,  which  represents  the  best 
correlation  which  can  be  obtained  between  measured  and 
predicted  results  without  adding  additional  parameters 
to  the  model. 

These  four  steps  make  up  the  subsections  of  the  following 
primary-wave  and  backing-wave  studies  for  the  characterization  of 
two  rc-500  transducers  face-to-face. 

E.l  CORRELATION  OF  PRIMARY  WWE  USING  BROADBAND  PULSE  INPUT 
The  primary  wave  is  correlated  using  the  broadband  pulse 
as  an  input  signal.  The  short  duration  of  the  pulse  makes  it 
the  best  signal  to  use  for  studying  the  primary  wave,  and  the 
fact  that  it  covers  a  relatively  broad  range  of  frequencies 
does  not  matter,  since  the  properties  of  the  wearpiate  and 
piezoelement  are  not  frequency-dependent. 

E.1.1  Experimental  Observations  Concerning  the  primary  Wave 
Fig.  El  shows  the  broadband  input  pulse  and  the  resulting 
output  signal  for  the  system.  The  first  wave  packet  of  the  output 
is  of  interest  here,  and  is  shown  on  an  expanded  time  scale  in 
Fig.  E2.  It  is  evident  from  these  plots  that  the  conversion  of 
input  to  output  in  the  system  is  not  a  trivial  process.  A  logical 
first  step  toward  accurately  predicting  this  output  signal  is  the 
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developnent  of  an  intuitive  feel  for  the  input-output  behavior  of 
the  system. 

Fig.  E3  is  a  schematic  of  the  propagation  of  stress  waves 
through  the  system.  The  solid  lines  represent  the  paths  of  the 
various  wavefronts  in  the  system  as  they  are  generated  at  the 
surfaces  of  the  input  piezoelement,  propagated  through  the  various 
transducer  layers,  reflected  from  impedance-mismatched  interfaces, 
and  finally  detected  at  the  surfaces  of  the  output  piezoelement. 
Since  this  section  deals  only  with  the  primary  wave,  waves 
entering  the  backing  layer  are  ignored  and  are  shown  as  dotted 
arrows  in  the  drawing. 

In  Chapter  2  (section  2.1)  it  is  shown  that  applying  a 
voltage  impulse  to  the  input  transducer  at  time  t=0  causes  four 
stress  impulses  to  be  generated  at  the  surfaces  of  the  piezo¬ 
element  at  that  instant.  Observe  from  Fig.  E3  that  the  stress 
impulse  produced  at  the  front  face  of  the  input  piezoelement 
(P  W  )  which  strikes  the  front  face  of  the  output  piezoelement 
(w^)  produces  a  positive  voltage  A(ll) ,  the  first  impulse  of 
the  output  signal.  This  wave  then  propagates  through  the  output 
piezoelement  and  strikes  its  rear  surface  (P^),  producing  a 
negative  voltage  A(21)  at  precisely  the  same  instant  as  the 
wave  from  the  rear  of  the  input  piezoelement  reaches  the  front 
surfe.ce  of  the  output  piezoelement.  This  last  wave  then 
continues  on  into  the  output  piezoelement,  producing  a  positive 
voltage  A(31)  as  it  strikes  the  rear  surface  (P^).  Fig.  E4 
shows  schematically  how  the  primary  wave  may  be  constructed  by 
simply  superposing  the  input  pulse  weighted  by  A(ll) ,  A(21), 


and  A(3i)  at  the  appropriate  arrival  times. 

Note  also  the  presence  of  multiple  reflections  produced  at 
the  wear plate- piezoelement  interfaces,  which  are  called  A(12) , 

A(13) ,  and  so  on.  Analysis  shows  that,  for  the  PC-500,  only  the 
first  two  of  these  reflections  make  a  significant  contribution 
(greater  than  1%  of  the  amplitude  of  the  first  impulse)  to  the 
output  signal.  It  is,  therefore,  the  nine  wavefronts  shown  in 
Fig.  E3  which  combine  to  produce  the  output  signal  seen  in  Fig.  E2. 

E.1.2  Preliminary  predictions  of  the  primary  Wave 

Fig.  E5  shows  the  broadband  pulse  response  of  the  system  as 
predicted  by  the  model,  using  the  measured  values  given  in  Fig.  C2 
as  inputs  to  the  model.  The  values  of  all  modeling  parameters 
used  in  each  simulation  are  given  in  Table  1;  in  general,  only 
the  particular  parameters  of  interest  will  be  mentioned  in  the 
discussion. 

Fig.  E5  shows  fair  agreement  between  predicted  and  measured 
output  signals.  The  predicted  initial  time  delay  is  only  0.06 
usee,  as  compared  to  the  measured  delay  of  0.22  usee.  This 
discrepancy  of  0.16  ysec  could  be  the  result  of  an  error  in  the 
calculated  wearplate  transit  time  (wavespeed  or  thickness),  since 
the  initial  time  delay  is  simply  twice  the  wearplate  transit  time, 
as  seen  in  Fig.  E3.  An  alternative  explanation  would  involve  a 
time  delay  introduced  at  either  the  couplant  layer  or  at  the 
wearplate-piezoelement  interfaces.  It  is  also  possible  that  the 
piezoelement  takes  some  time  to  respond  to  stress  and  produce  a 
voltage. 


Other  than  the  time  delay,  the  most  obvious  discrepancy  in 
Fig.  £5  occurs  after  2  ysec,  with  the  appearance  of  disturbances 


in  the  measured  signal  that  do  not  show  up  in  the  predicted 
signal.  Because  of  the  similarity  between  the  first  2  usee  of 
the  two  signals,  however,  it  seems  that  the  three  main  peaks 
(corresponding  to  A(ll),  A(21),  and  A(31)  in  Fig.  £3)  fall  at 
about  the  right  times  and  with  reasonably  correct  amplitudes  in 
the  predicted  output.  The  only  source  for  disturbances  after  2 
sec  would  be  the  multiple  reflections  A(12) ,  A(13) ,  A( 22 ) , 

A(23),  A (32)  and  A(33)  which  are  generated  at  the  wearplate- 
piezoelement  interfaces.  These  do  not  appear  in  the  predicted 
output  of  Fig.  E5  because  of  the  extremely  short  transit  time 
of  the  wearplates  (0.04  usee),  so  that  the  predicted  multiple 
reflections  are  only  0.16  usee  apart,  whereas  those  in  the 
measured  output  are  almost  0.6  ysec  apart.  It  seems  likely 
that  a  time  delay  is  introduced  as  the  waves  are  reflected  from 
the  wearplate-piezoelement  interfaces. 

£.1.3  Effect  of  Wearplate-Piezoelement  Time  Delays 
on  the  Primary  Wave 

Fig.  E6  illustrates  the  effect  of  including  a  time  delay, 

S,  in  the  simulation,  showing  predicted  output  signals  for  values 
of  8  from  0.08  to  0.32  usee.  These  curves  do  indeed  show 
disturbances  beyond  2  usee,  as  expected.  In  particular,  the 
signal  for  which  8“C*24  ysec  bears  a  marked  resemblance  to  the 
measured  output  signal  of  Fig.  £5. 

It  is  interesting  to  note  the  effects  of  introducing  these 


time  delays.  l£  only  the  reflected  waves  are  delayed,  as  is  the 
case  in  Fig.  E6,  the  arrival  times  of  wavefronts  A(ll) ,  A(21) ,  and 
A (31)  remain  unchanged,  because  none  of  these  waves  undergoes  any 
reflection  (see  Fig.  E3) .  The  time  interval  between  A( 11)  and 
A(12)  is  affected,  however,  and  so  is  the  interval  between  A(12) 
and  A(13) ,  A(21)  and  A(22) ,  and  so  on. 

Now  consider  delaying  the  transmitted  waves  at  the  wearplate- 
piezoelement  interfaces.  This  has  the  effect  of  delaying  every 
wavefront  by  precisely  the  same  amount,  since  each  wave  is 
transmitted  across  each  interface  once  and  only  once.  A  time 
delay,  6,  of  0.08  psec  for  the  transmitted  waves  therefore 
resolves  the  discrepancy  in  the  initial  time  delay  noted  earlier. 
It  is  concluded,  then,  that  waves  reflected  from  the  wearplate- 
piezoelement  interfaces  are  delayed  by  0.24  psec,  and  that  waves 
transmitted  across  these  interfaces  are  delayed  by  0.08  psec. 

E.1.4  Effect  of  the  Time  Constant  RC0  on  the  Primary  Wave 

The  predicted  output  signal  of  Fig.  E5  is  based  on  the 
assumption  that  the  output  piezoelement  produces  a  voltage 
proportional  to  the  integral  of  the  stress  inside  it.  This  is  a 
good  approximation  if  the  resistance  across  the  piezoelement  is 
large,  because  electrically  the  piezoelement  behaves  like  a 
capacitor,  and  will  continue  to  accumulate  voltage  if  no  path  is 
provided  for  current  to  flow  from  one  plate  of  the  capacitor  to 
the  other.  If,  however,  there  is  a  relatively  small  resistance 
across  the  piezoelement,  the  voltage  will  decay  with  time 
according  to  the  factor  exp(-t/RC0) ,  where  C0  is  the  capacitance 


of  the  piezoelement,  and  R  is  the  resistance  across  it.  The 
output  voltage  is  then  no  longer  proportional  to  the  integral  of 
the  stress  inside  the  piezoelement,  but  to  the  convolution  of 
stress  with  an  exponentially-decaying  function  of  time  (see 
Chapter  2,  section  2.5). 

If  the  value  of  RC0  in  the  simulation  is  decreased,  the  shape 
of  the  output  signal  changes  drastically,  as  seen  in  Fig.  E7. 

A  time  constant  which  is  of  the  same  order  of  magnitude  as  other 
time  parameters  in  the  system  allows  the  output  voltage  to  decay 
in  "real  time,"  so  that  it  does  not  accumulate  as  fast  as  it  would 
if  RC  were  larger.  If  the  time  constant  is  made  snail  enough, 
the  voltage  never  accumulates  at  all,  and  the  output  voltage  is 
proportional  to  the  stress  itself  rather  than  to  the  integral  of 
stress,  as  shown  in  Fig.  E7(a).  The  other  three  curves  of  Fig. 

E7  show  that  increasing  the  time  constant  has  two  effects:  it 
increases  the  amplitude  of  the  peak  at  1.1  usee,  and  decreases  the 
amplitude  of  the  peak  at  1.8  usee.  Comparing  these  predicted 
output  signals  to  the  measured  output  of  Fig.  E5  indicates  that 
the  value  of  the  time  constant  for  the  system  under  consideration 
is  about  1.0  usee. 

E.1.5  Best  Prediction  of  the  Primary  Wave 

The  best  predicted  output  signal  corresponding  to  a  broadband 
pulse  input  is  presented  for  comparison  with  the  measured  output 
signal  in  Fig.  E8.  The  important  parameters  used  in  the 
simulation  are  the  time  delays  6*0.24  usee  and  6-0.08  usee,  and 
the  time  constant  RCO*1.0  usee.  Qualitatively,  the  measured  and 
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predicted  curves  appear  to  correlate  very  well.  Quantitatively, 

2 

they  indicate  that  the  amplitude  parameter  h  cys  has  the  value 
5.47  x  1015  N/to3  for  this  case. 

E.2  CORRELATION  OF  THE  BACKING  V&VE  USING  TONE  BURST  INPUT 
It  is  easier  to  correlate  the  backing  wave  using  a  tone 
burst  input  signal  for  several  reasons.  First,  the  attenuation 
of  the  backing  material  used  in  the  transducers  is  frequency- 
dependent,  and  the  narrow  bandwidth  of  the  tone  burst  signals 
enables  the  determination  of  attenuation  at  several  different 
frequencies.  Second,  it  is  easier  to  determine  the  exact  time 
delay  associated  with  the  backing  wave  when  a  tone  burst  is  used, 
since  the  time  delay  directly  affects  the  phase  angle  between  the 
primary  wave  and  the  backing  wave.  Although  the  primary  wave  has 
been  correlated  with  a  broadband  input  pulse,  the  optimum  values 
for  6,  3,  and  RC0  have  been  found  to  be  the  same  for  tone  bursts 
across  the  entire  frequency  range  (0. 4-2.0  MHz)  as  for  the  pulse. 

E.2.1  Experimental  Observations  Concerning  the  Backing  Wave 
Fig.  E9  shows  the  measured  0.4  MHz  tone  burst  input  and 
output  signals  for  two  PC-500  transducers  face-to-face.  Observe 
that  the  input  signal  takes  three  full  cycles  to  settle  down  to  a 
steady  amplitude,  and  that  this  transient  shows  up  in  the  output 
signal  as  well.  Note  also,  however,  that  the  output  shows  other 
disturbances  at  about  10  ysec  after  the  input  is  first  applied. 
Recall  from  Fig.  El  that  10  ysec  is  the  observed  time  delay 
associated  with  the  backing  layer,  so  that  both  the  disturbance  at 


j 


10  usee  and  ne  reduced-amplitude  oscillation  from  25  to  35  usee 
are  due  to  waves  traveling  through  the  backing  layer. 


E.2.2  Preliminary  Predictions  Including  the  Backing  Wave 

Fig.  E10  shows  the  measured  output  signal  and  the  predicted 
output,  the  latter  obtained  by  using  the  measured  values  for 
thickness  and  attenuation  of  the  backing  material  at  0.4  MHz. 
Observe  that  the  first  three  cycles  of  the  predicted  output 
compare  very  favorably  with  the  measured  output,  but  that  when  the 
backing  wave  arrives  at  about  10  usee,  the  amplitude  of  the 
predicted  output  is  greatly  reduced.  One  would  suspect  that  the 
value  of  backing  attenuation  is  too  high — that  the  predicted 
amplitude  of  the  backing  wave  is  too  small.  The  plot  shows, 
however,  that  when  everything  except  the  backing  wave  has  died  out 
(30  usee) ,  the  predicted  backing  wave  actually  has  a  larger 
amplitude  that  the  measured  one.  Observe  also  that  the  last  full- 
amplitude  peak  of  the  measured  output  (before  the  amplitude  begins 
to  taper  off)  occurs  at  about  31  usee,  but  only  at  29.5  usee  for 
the  predicted  output.  These  observations  indicate  that  the 
predicted  backing  wave  arrives  about  1.5  usee  too  early, 
necessitating  the  introduction  of  a  time  delay  similar  to  those 
used  in  the  previous  section. 

E.2.3  Effect  of  Time  Delays  on  the  Backing  Wave 

Fig.  Ell  shows  the  predicted  output  signals  for  several 
values  of  time  delay,  y ,  introduced  as  waves  propagate  through 
the  backing  layer  and  reflect  from  the  backing-potting  interface. 
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These  time  delays  were  simulated  by  increasing  the  value  of 
backing  thickness  in  the  input  data  for  the  model,  in  order  to 
also  illustrate  the  dependence  of  the  amplitude  of  the  tacking 
wave  on  the  thickness  of  the  backing  layer  itself.  The 
difference  in  amplitude  over  the  range  of  backing  thickness  shown 
in  Fig.  Ell  is  quite  small,  but  one  can  see  that  the  backing  wave 
decreases  slightly  in  amplitude  for  increasing  backing  thickness. 

The  more  important  effect  seen  in  Fig.  Ell  is  the  rather 
strong  dependence  of  the  steady-state  output  amplitude  on  the 
time  delay.  Note  that  the  variation  in  steady-state  amplitude  as 
a  result  of  changes  in  the  time  delay  is  much  greater  than  the 
variation  in  the  amplitude  of  the  backing  wave  itself.  This  is 
due  to  the  effect  of  the  time  delay,  y  ,  on  the  phase  angle  between 
the  primary  wave  and  the  backing  wave.  A  change  of  0.1  ysec  in 
the  time  delay  is  equivalent  to  a  15  degree  shift  in  phase  at  fl.4 
MHz,  so  that  over  the  range  of  time  delays  in  Fig.  Ell,  the  phase 
angle  between  the  primary  and  backing  waves  varies  by  45  degrees. 
This  can  make  the  difference  between  the  backing  wave  having  a 
strong  effect  on  the  steady-state  amplitude,  as  in  plot  (a),  and, 
as  in  plot  (d) ,  the  steady-state  signal  being  virtually 
independent  of  the  backing  wave,  depending  only  on  the  primary- 
wave  amplitude. 

Fig.  E10  shows  quite  graphically  the  extreme  case,  that  of 
nearly  totally  destructive  interference  between  the  backing  and 
primary  waves.  Note  that  the  output  shows  a  drastic  decrease  in 
amplitude  due  to  the  arrival  of  the  backing  wave  at  10  ysec,  and 
that  it  exhibits  a  sort  of  phase  reversal  at  about  25  ysec,  when 
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the  primary  wave  dies  out  and  the  backing  wave  takes  over. 
Comparing  Fig.  E10  and  Fig.  Ell (a)  shows  the  sort  of  effect  that 
relatively  small  changes  in  the  time  delay  can  have  on  the  output 
of  the  system.  This  effect  is,  of  course,  frequency-dependent  as 
well,  because  the  phase  angle  corresponding  to  a  particular  time 
delay  is  a  function  of  the  frequency. 

It  is  concluded,  then,  that  there  is  a  time  delay  introduced 
as  waves  are  propagated  through  the  backing  layer  and  reflected 
from  the  backing-potting  interface.  This  time  delay  is 
particularly  significant  in  that  it  affects  the  phase  angle 
between  the  primary  wave  and  the  backing  wave,  and  consequently 
can  affect  the  steady-state  output  amplitude  to  a  great  extent. 
More  discussion  of  the  steady-state  output  amplitude  is  given  in 
Appendix  G. 

E.2.4  Effect  of  Backing  Attenuation  on  the  Backing  Wave 
Fig.  E12  shows  the  result  of  varying  the  attenuation 
parameter,  cxD,  of  the  backing  material.  As  expected*  the  only 
effect  on  the  output  is  a  change  in  the  amplitude  of  the  backing 
wave  and,  to  a  lesser  degree,  the  steady-state  anplitude.  From 
these  two  studies  it  is  concluded  that  the  best  correlation 
between  theory  and  experiment  at  0.4  MHz  is  obtained  for  the 

values  V*1.40  usee  and  a  =1.8  nepers/ cm. 

B 

E.2.5  Best  Prediction  Including  the  Backing  Wave 

Fig.  E13  shows  the  best  predicted  output  signal  for  a  P.4 
MHz  tone  burst  input,  along  with  the  measured  output  for 
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comparison.  Once  again,  the  correlation  is  impressive,  and  the 
value  of  h2Co/S  is  found  to  be  3.80  x  1015  N/m3,  only  slightly 
smaller  than  that  found  from  the  broadband  pulse  results. 

To  complete  the  characterization,  the  backing  wave  can  be 
correlated  at  as  many  frequencies  as  desired,  using  the  same 
procedure  as  described  in  this  section.  For  example,  Fig.  E14 
shows  the  measured  input  and  output  signals  for  a  0.6  MHz  tone 
burst  input,  and  Fig.  E15  shows  the  best  prediction,  obtained 

from  the  parameter  values  7*1.40  visec  and  afi=2.5  nepers/am.  The 

2  i  s  3 

amplitude  factor  h  Cc/S=4.72  x  10  N/m  for  this  case. 

The  values  of  the  modeling  parameters  for  the  PC-500,  as 
well  as  for  all  of  the  other  transducers,  are  presented  in  Tables 
2  through  6  and  are  discussed  in  Chapter  4. 


(b)  Ti-.e  (-see) 

Fig.  E2  Measured  broadband  pulse  (a)  input  and  (b)  outp 
signals  for  two  FC-500  transducers  clamped 
face-to-face  (expanded  time  scale). 


Fig.  E3  Schematic  of  (a)  the  propagation  of  stress  waves  and  (b)  the  production 
of  output  voltage  In  an  ultrasonic  test  system  consisting  of  two 
transducers  clamped  face-to-face. 


(b)  Time  (usee) 

Fig.  E5  Comparison  of  (a)  measured  and  (b)  predicted  output 
signals  corresponding  to  broadband  pulse  input  for 
two  FC-500  transducers  clamped  face-to-face 
(5-0,  6-0,  RCo-0). 


Time  (usee) 


Fig.  E6  Predicted  broadband  pulse  response  of  two 

FC-500  transducers  clamped  face-to-face  as  a 
function  of  the  time  delay  8. 


Time  (psec) 

Fig.  E8  Comparison  of  (a)  measured  and  (b)  predicted  output 
signals  corresponding  to  broadband  pulse  input  for 
two  FC-500  transducers  clamped  face-to-face 
(6*0.08  psec, 6*0. 24  psec,  RC  *1.0  psec) 


Time  (..sec) 

Fig.  E10  Comparison  of  (a)  measured  and  (b)  predicted 
output  signals  corresponding  to  0.4  MHz  tone 
burst  input  for  two  FC-500  transducers  clamped 
face-to-face  (A-0,  a  *1.6  nepers/cm). 


3 

2 

1 


Tine  (usee) 


Fig.  E13  Comparison  of  (a)  measured  and  (b)  predicted  output 
signals  corresponding  to  0.4  MHz  tone  burst  input 
for  two  FC-500  transducers  clamped  face-to-face 
(A-1.40  psec,  a  *1.8  nepers/cm). 
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APPENDIX  F 


DISCUSSION  OF  MEASURED  STEADY- STATE  OUTPUT  AMPLITUDE 


All  of  the  transducers  were  clamped  face-to-face  and  tested 
with  100  V  peak-to-peak  continuous  sinusoidal  input  signals.  The 
raw  data  are  summarized  in  Table  7  ,  and  are  plotted  against 
frequency  in  Figs.  FI,  F2,  and  F3.  These  three  plots  illustrate 
the  effects  of  three  different  transducer  parameters — backing 
impedance,  wearplate  impedance,  and  backing  thickness — on  the 
steady-state  output  amplitude  of  an  ultrasonic  test  system.  Each 
of  these  parameters  is  discussed  in  detail  in  one  of  the  following 
sections. 


F.l  EFFECT  OF  BACKING  IMPEDANCE  ON  STEADY-STATE  OUTPUT  AMPLITUDE 

Fig.  FI  shows  the  steady-state  output  amplitude  vs.  frequency 
for  both  the  unbacked  transducer  (AE-100)  and  the  PC-500,  which 
has  an  impedance-matched  backing  layer.  This  plot  illustrates 
very  clearly  how  the  presence  of  a  matched  backing  layer  affects 
the  frequency  response  of  the  transducer.  The  FC-500  shows  a 
maximan  peak  of  about  10  volts  at  1.9  MHz,  but  the  AE-100  has  a 
peak  output  of  over  50  volts,  at  about  1.6  MHz.  The  very  large 
peak  in  the  response  of  the  unbacked  transducer  corresponds  to  a 
condition  in  the  piezoelement  known  as  standing-wave  resonance. 

Hie  standing-wave  (or  half-wave)  resonance  condition  occurs 
whenever  the  frequency  of  the  oscillations  within  the  piezoelement 
is  equal  to  an  even  multiple  of  the  reciprocal  of  the  piezoelement 
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transit  time,  so  that  waves  reflected  back  into  the  piszoelement 
from  the  surfaces  of  the  piezoelement  interfere  constructively 
(in  phase)  with  the  waves  incident  upon  the  surfaces,  A  matched 
backing  layer  absorbs  the  incident  waves  at  one  surface  of  the 
piezoelement,  so  that  they  are  not  reflected  Lack  into  the  piezo¬ 
element,  and  standing  waves  cannot  develop.  It  is  not  surprising, 
then,  that  a  strong  resonant  peak  exists  in  the  unbacked  trans¬ 
ducer,  nor  that  the  peak  is  effectively  damped  out  by  the  presence 
of  the  backing  layer  in  the  FC-500.  This  matched  backing  layer 
gives  the  transducer  a  much  more  flat  response,  so  that  its  useful 
bandwidth  is  increased  by  a  factor  of  two. 


F.2  EFFECT  OF  WEARPLATE  IMPEDANCE  CN  STEADY-STATE  OUTPUT  AMPLITUDE 

Fig.  F2  shows  the  steady-state  response  for  the  transducer 
with  aluminum  wearplate  (AE-101)  compared  to  that  for  the  standard 
FC-50C.  The  presence  of  the  aluninum  wearplate  appears  to  increase 
the  steady-state  amplitude  by  a  significant  amount,  especially  at 
higher  frequencies  (greater  than  1.5  MHz).  In  fact,  the  increase 
in  amplitude  is  so  great  at  the  high  end  that  the  peak  amplitude 
for  the  AE-101  comes  at  about  2.3  MHz,  compared  to  1.9  MHz  for  the 
FC-500 . 

There  appear  to  be  two  different  factors  involved  here. 

The  difference  between  the  steady-state  amplitudes  for  the  two 
transducers  at  relatively  low  frequencies  (below  1.5  MHz)  is  due 
chiefly  to  the  improved  impedance-match  between  wearplate  and 
piezoelement  in  the  AE-101,  which  increases  the  transmission 
coefficient  for  stress-wave  propagation  across  the  wearplate- 
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piezoelement  interfaces  in  the  system.  Thus,  a  greater  portion 
of  the  stress-wave  amplitude  which  is  generated  at  the  input 
piezoelement  actually  arrives  at  the  output  piezoelement,  with  the 
result  that  the  output  amplitude  is  increased.  This  effect  alone 
would  produce  two  curves  of  steady-state  output  amplitude  that 
differed  by  a  constant  multiplicative  factor  at  all  frequencies. 
This  is  obviously  not  the  complete  picture,  though,  for  it  is  clear 
that  the  alixninum  wearplate  has  some  effect  on  the  resonant 
frequency  of  the  transducer  as  well. 

In  general,  an  increase  in  resonant  frequency  corresponds  to 
a  decrease  in  piezoelement  transit  time.  Since  there  are  large 
time  delays  8  for  reflected  waves  which  approach  the  wearplate- 
piezoelement  interfaces  from  the  outside  of  the  piezoelement 
(see  Chapter  4,  subsection  4.1.2),  it  seems  probable  that  there  are 
time  delays  for  waves  which  strike  the  surfaces  of  the  piezoelement 
from  the  inside.  Although  these  time  delays  could  only  increase 
the  piezoelement  transit  time,  recall  that  the  value  of  6  for  the 
AE-101  was  about  half  of  the  value  for  the  PC-500.  The  difference 
between  the  period  corresponding  to  1.9  MHz  and  that  corresponding 
to  2.3  MHz  is  0.39  usee,  and  the  difference  between  the  wearplate- 
piezoelement  time  delays  fer  the  PC-500  and  the  AE-101  is  0.115 
usee,  larger  by  slightly  more  than  20%.  It  seems  perfectly 
reasonable,  then,  that  the  increase  in  resonant  frequency  for  the 
aluminum  backing  layer  is  due  to  smaller  time  delays  at  the  wear- 
plate-piezoelement  interfaces. 

In  conclusion,  then,  the  AE-101  shows  a  consistently  larger 


steady-state  amplitude  than  the  PC-500,  due  to  the  relatively 
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good  impedance  match  between  the  wearplate  and  piezoelement  in 
the  AE-101.  In  addition,  the  resonant  frequency  of  the  AE-101  is 
higher  than  that  for  the  PC-500,  due  to  the  smaller  time  delays 
introduced  at  the  wearplate-piezoelement  interfaces  in  the  AE-101. 

F.3  EFFECT  OF  BACKING  THICKNESS  ON  STEADI-STATE  OUTPUT  AMPLITUDE 

Fig.  F3  shows  a  plot  of  steady-state  output  amplitude  vs. 
frequency  for  the  AE-102,  AE-103,  and  PC-500  transducers.  These 
transducers  are  allegedly  identical,  except  for  the  thickness  of 
their  backing  layers.  Notice  that  the  PC-500  and  the  AE-102  have 
slightly  greater  amplitude  than  the  AE-103  transducer  over  almost 
the  entire  frequency  range  of  the  experiments  (0. 4-3.0  MHz). 

Observe  also  that  the  peak  amplitude  of  the  prototype  transducers 
(AE-102  and  AE-103)  occurs  at  about  1.6  MHz,  while  the  peak  in  the 
response  of  the  standard  PC-500  comes  near  1.9  MHz.  In  addition, 
the  shapes  of  the  response  curves  of  the  AE-102  and  AE-103  are 
qualitatively  similar  to  one  another,  yet  the  curve  for  the  PC-500 
is  significantly  different.  These  observations  indicate  some 
variations  in  construction  between  the  transducers. 

Since  the  behavior  of  the  PC-500  is  so  different  from  the 
AE-102  and  AE-103,  the  PC-500  cannot  really  be  compared  with 
the  two  other  transducers  here.  It  is  still  possible,  though,  to 
draw  some  conclusions.  Recall,  first,  that  in  Chapter  4 
(subsection  4.2.3)  it  was  noted  that  the  amplitude  of  the  pulse 
response  of  a  transducer  appears  to  decrease  for  increasing  backing 
thickness.  This  is  also  the  case  here,  where  the  AE-103  (1.905  cm 
(0.750  in)  backing)  shows  a  characteristically  lower  steady-state 


amplitude  than  the  AE-102  (1.270  an  (0.500  in)  backing).  Thus,  the 
correlation  of  backing  thickness  with  output  amplitude  is  substan¬ 
tiated.  The  relatively  low  resonance  frequency  of  the  prototype 
transducers,  on  the  other  hand,  is  much  lower  than  would  be  expected 
based  on  the  results  for  the  AE-101  transducer  discussed  in  the 
previous  section. 

In  conclusion,  the  steady-state  amplitude  is  observed  to 
decrease  for  increasing  backing  thickness.  No  effect  of  backing 
on  the  resonance  frequency  is  observed;  the  important  variation  of 
resonance  frequency  occurs  between  the  PC-503  and  the  others,  not 
between  identical  transducers  differing  only  in  backing  thickness. 
Apparently,  the  piezoelement  used  in  the  FC-500  is  different  from 
that  in  the  prototype  transducers. 
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steady-state  output  amplitude  vs.  frequency  for  AE-101  and 
ransducers  corresponding  to  continuous-wave  input. 


of  steady-state  output  amplitude  vs.  frequency  for  AE-103,  AE-102 
FC-500  transducers  corresponding  to  continuous-wave  input. 


APPENDIX  G 


SIMULATED  EFFECTS  OF  WEAR  PLATE  IMPEDANCE,  BACKING  IMPEDANCE, 
AND  PIEZOELEMENT  THICKNESS  ON  TRANSDUCER  RESPONSE 


In  this  appendix,  the  effects  of  three  different  parameters 
on  the  broadband  pulse  response  of  a  system  are  determined.  The 
system  consists  basically  of  two  PC-500  transducers  face-to-face, 
except  for  variations  in  the  particular  parameters  being  studied. 


G.l  EFFECT  OF  WEARPLATE  IMPEDANCE 


Fig.  G1  consists  of  four  predicted  output  signals  which  show 
the  result  of  varying  the  wearplate  impedance  in  the  system.  The 
standard  FC-500  corresponds  to  case  (d),  vrfiere  zyz?=2.5.  Case 
(b)  represents  transducer  AE-101,  which  has  an  aluminun  wearplate 
rather  than  the  standard  titanium  carbide  wearplate  of  the  PC-500. 
Case  (c)  corresponds  to  a  system  in  which  the  wearplates  are 
impedance-matched  to  the  piezoelements,  and  in  case  (a),  the  ratio 
of  wearplate  impedance  to  piezoelement  impedance  is  the  reciprocal 
of  the  standard  ratio  for  the  FC-500. 

Notice  that  curves  (a)  and  (d)  appear  to  be  identical.  This 
is  at  first  surprising,  but  recall  that  waves  traversing  the  two 
wearplates  must  cross  two  wearpiate-piezoelement  interfaces,  once 
from  piezoelement  to  wearplate,  and  once  from  wearplate  to 
piezoelement.  Note  that  the  transmission  coefficient  for  a 
wearplate-to-piezoelement  interface  vrfiere  Zw  is  half  of  Zp  is 
identical  to  the  transmission  coefficient  for  a  piezoelement-to- 
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wearplate  interface  where  Zt,  is  twice  Z  .  This  means  that  the 
product  of  the  transmission  coefficients  for  the  two  interfaces 
should  be  the  same,  regardless  of  the  order  in  which  the  inter¬ 
faces  are  traversed. 

The  other  interesting  observation  to  be  made  is  that  the 
"efficiency"  of  the  system  drops  significantly  for  large  impedance 
mismatch  between  piezoelement  and  wearplate.  An  impedance  ratio 
of  0.75,  as  in  case  (b) ,  decreases  the  efficiency  only  slightly, 
but  a  decrease  of  almost  30%  is  observed  for  Z^Zp«0.4. 

In  conclusion,  then,  the  simplest  and  largest-amplitude 
output  signal  is  obtained  by  impedance-matching  the  wearplate  and 
pioezoelement  in  an  ultrasonic  transducer. 

G.2  TOE  EFFECT  OF  BACKING  IMPEDANCE 

Fig.  G2  shows  four  curves  which  predict  the  result  of  varying 

the  backing  impedance  of  the  system.  The  standard  case,  of  course, 

corresponds  to  Z  /Z  =1.0,  where  the  backing  is  impedance-matched 
‘  B  P 

to  the  piezoelement. 

The  most  obvious  feature  of  Fig.  G2  is  the  resonance 
which  occurs  when  the  backing  is  mismatched.  It  is  interesting, 
however,  that  the  resonant  frequencies  are  not  all  the  same.  Case 
(a)  seems  to  have  a  fairly  high  frequency  of  oscillation,  while 
cases  (c)  and  (d)  show  a  much  lower  frequency. 

In  case  (a) ,  the  backing  has  a  lower  impedance  than  the 
piezoelement.  This  means  that  a  wave  in  the  piezoelement  which 
strikes  the  piezoelement-backing  interface  will  be  partly  reflected 
and  partly  transmitted.  More  importantly,  however,  the  reflected 


wave  will  be  inverted  as  compared  to  the  incident  wave — i ts 
amplitude  will  be  a  negative  fraction  of  the  incident  amplitude 
(see  eqns.  (2-15)).  When  the  reflected  wave  strikes  the 
piezoelement-wearplate  interface,  it  will  then  produce  a  voltage 
pulse  which  is  opposite  in  sign  to  the  pulse  produced  at  the 
piezoelement-backing  interface  one  piezoelement  transit  time 
earlier.  Thus,  the  predominant  frequency  in  plot  (a)  is  deter¬ 
mined  by  the  piezoelement  transit  time  (0.27  ysec).  Note  the 
large  peak  at  about  2  ysec  and  its  negative  counterpart  at  about 
2.5  ysec  as  evidence  of  this. 

Figs.  G2(c)  and  (d)  are  fairly  similar  to  one  another  except 
for  the  rate  of  decay  of  the  resonance.  They  both  correspond  to 
cases  in  which  the  backing  has  a  larger  impedance  than  the 
piezoelement,  so  that  the  wave  reflected  from  the  piezoelement- 
backing  interface  has  the  same  sign  as  the  incident  wave.  This 
means  that,  rather  than  a  negative  peak  appearing  in  the  output 
shortly  after  each  positive  peak,  the  positive  peaks  are  repeatea 
at  slightly  reduced  amplitude.  Then,  since  the  wearplate 
impedance  is  also  much  larger  than  the  piezoelement  impedance,  the 
wave  is  reflected  from  the  piezoelement-wearplate  interface  as 
well,  with  a  positive,  further- reduced  amplitude.  This  positive 
reinforcement  accounts  for  the  broadening  of  the  peaks  in  these 
two  output  signals. 

The  0.5  MHz  resonance  frequency  in  these  last  two  cases  of 
Fig.  G2  corresponds  to  repeated  reflection  back  and  forth 
between  the  piezoelement-backing  interfaces  of  the  two  transducers. 
The  round-trip  transit  time  is  almost  2  usee,  which  is  of  course 
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the  period  associated  with  0.5  MHz  resonance.  The  reason  that 
the  resonance  takes  longer  to  decay  in  case  (d)  than  in  case  (c) 
is  that  the  reflection  coefficient  at  the  piezoelement-backing 
interface  is  higher  for  larger  backing  impedance.  This  means 
that  more  of  the  energy  will  be  reflected  back  into  the 
piezoelement  by  a  very  stiff  backing,  so  that  it  will  take  more 
cycles  between  the  piezoelement-backing  interfaces  for  the 
amplitude  of  the  reflected  wave  to  fall  to  a  negligible  level. 

It  is  therefore  clear  that  the  only  backing  layer  that  is 
effective  in  damping  transducer  resonance  is  the  impedance-matched 
backing. 

G.3  THE  EFFECT  OF  PIEZOELEMENT  THICKNESS 

Fig.  G3  shows  four  predicted  pulse  output  signals  for  four 
systems  with  piezoelements  of  varying  thickness.  Case  (c) 
corresponds  to  the  standard  PC-500  piezoelement. 

The  most  striking  feature  of  Fig.  G3  is  the  variation  in 
"smoothness"  of  the  output  signals.  Plot  (a)  is  very  rough  and 
jagged,  while  plot  (d)  shows  almost  no  detail  whatsoever.  This  is 
a  consequence  of  the  fact  that  the  output  voltage  is  approximately 
proportional  to  the  integral  of  the  stress  inside  the  piezoelement 
of  the  output  transducer.  A  thicker  piezoelement  integrates  over 
a  longer  time  interval,  so  that  the  small,  abrupt  variations  in 
the  stress  waveform  make  a  lesser  contribution  to  the  output 
voltage. 

Another  consequence  of  integrating  over  a  longer  time 
interval  is  that  the  amplitude  of  the  integral  increases  as  the 
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time  interval  increases.  Fig.  G3  shows  chat  increasing  the 
piezoelement  thickness  eightfold  results  in  a  factor  of  twenty 
increase  in  output  amplitude.  This  increase  is  highly  desirable, 
but  the  tradeoff  is  that  the  bandwidth  of  the  output  signal  is 
reduced  by  lengthening  the  piezoelement.  Similarly,  a  very  thin 
piezoelement  produces  a  voltage  signal  that  is  very  close  to  the 
actual  stress  waveform,  but  it  is  of  very  low  amplitude.  The 
PC-500  seems  to  be  well  suited  to  the  particular  broadband  input 
pulse  used  in  this  study,  giving  enough  detail  in  the  output 
without  creating  a  jagged  signal  that  is  difficult  to  interpret 
or  quantify. 

In  conclusion,  then,  the  proper  piezoelement  thickness  will 
provide  a  useful  amount  of  detail  in  the  output  signal,  while 
also  maintaining  a  reasonable  amplitude  level. 
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Fig.  G3  Predicted  broadband  pulse  response  of  two  FC-500 
transducers  clamped  face-to-face  as  a  function  of 
piezoelement  thickness. 


appendix  h 


EXPERIMENTAL  STUDY  OF  VWE  INTERACTION  AT  INTERFACES 

This  appendix  is  a  discussion  of  experimental  results  which 
deal  with  the  behavior  of  stress  waves  at  the  interfaces  between 
the  various  lasers  in  an  ultrasonic  test  system.  The  effect  of 
a  thin  layer  of  ultrasonic  couplant  on  the  output  of  the  system  is 
investigated  for  both  impedance-matched  and  impedance-mismatched 
interfaces;  the  effect  of  a  thin  epoxy-bonding  layer  is  also 
discussed. 

H.l  BEHAVIOR  AT  IMPEDANCE-MISMATCHED  COUPLANT  LAYER 

Fig.  HI  shows  two  measured  broadband  pulse  output  signals. 

A  1.270  cm  (0.500  in)  thick  aluminum  plate  serves  as  the  specimen 
in  each  system,  and  two  different  types  of  transducers  are  used — 
AE-103  transducers  for  system  A,  and  AE-101  transducers  for 
system  B.  The  significant  difference  between  the  transducers 
is  that  the  AE-103  has  a  titanium  carbide  wearplate,  and  the 
AE-101  has  an  alum i nun  wearplate.  Thus  one  wearplate  is 
impedance-matched  to  the  specimen,  and  the  other  is  mismatched. 

The  two  output  signals  are  quite  similar  except  for  two 
important  differences.  First,  observe  the  interval  between  7  and 
9  ysec.  System  B  shows  only  a  slight  perturbation  at  7.5  ysec, 
followed  quickly  by  another  plateau,  but  system  A  shows  a  distinct 
reflection  arriving  at  that  time,  if  this  reflection  is  caused 
by  the  largest  peak  in  the  signal  (which  starts  at  about  3.5  usee) , 
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there  appears  to  be  a  time  delay  of  approximately  4  usee  between 
the  arrival  of  the  first  peak  and  the  arrival  of  the  reflection. 
The  theoretical  transit  time  of  a  1.270  cm  (0.500  in)  thick 
a 1  uninun  plate  is  2.04  usee,  so  that  a  measured  time  delay  of 
4  usee  is  within  2t  of  the  predicted  value.  Therefore,  the 
reflection  at  7.5  usee  in  System  A  is  clearly  caused  by  the 
impedance  mismatch  at  the  wearplate-specimen  interface. 

The  second  difference  between  the  plots  is  the  difference  in 
overall  amplitude.  System  B  has  a  peak  amplitude  of  over  0.S 
volts,  while  the  peak  value  for  system  A  is  barely  0.4  volts. 

This  may  also  be  attributed  to  the  reflection  at  the  impedance- 
mismatched  interfaces — a  portion  of  the  input  energy  in  system  A 
is  reflected  back  into  the  input  wearplate  by  the  specimen,  and 
another  portion  is  reflected  back  into  the  specimen  by  the  output 
wearplate.  Thus,  only  a  fraction  of  the  energy  radiated  by  the 
input  piezoelement  is  actually  transmitted  to  the  output 
piezoelement,  resulting  in  a  reduced  measurement  of  output 
amplitude  as  compared  to  system  B.  This  occurs  to  all  waves, 
however,  so  that  the  shape  of  the  output  signal  is  not  really 
affected,  except  for  the  reflection. 

Using  the  impedance  values  from  Fig.  C2,  eqns.  (2-17) 
predict  that  the  amplitude  should  be  reduced  by  a  factor  of  0.696. 
The  observed  reduction  is  0.67,  accurate  to  within  4%  of  the 
predicted  value. 

It  seems  clear,  then,  that  the  presence  of  a  couplent  layer 
does  not  affect  the  validity  of  the  impedance  model  for  wove 
interaction  at  surfaces.  In  addition,  no  time  delay  appears  to  be 


introduced  when  waves  ere  reflected  from  impedance-mismatched 
interfaces  where  a  couplant  layer  is  present.  Further  experiments 
and  some  wave  propagation  simulations  of  the  system  should  be 
performed  to  confirm  these  conclusions. 

t 

H. 2  BEHAVIOR  AT  IMPEDANCE-MATCHED  COUPLANT  LAYER 

In  order  to  isolate  the  effect  of  the  couplant  layer,  two 
specimens  have  been  prepared  {see  Fig.  C5) .  One  is  simply  the 

I. 270  cm  (0.S00  in)  thick  aluminum  plate  of  the  preceding 
experiment,  and  the  other  consists  of  two  0.635  cm  (0.250  in) 
thick  aluminum  plates  with  a  layer  of  couplant  between  them.  The 
measured  broadband  pulse  output  signals  corresponding  to  these 
two  systems  are  shown  in  Fig.  H2. 

The  curves  are  strikingly  similar,  qualitatively  as  well  as 
quantitatively.  It  appears  that  there  is  no  ;time  delay  introduced 
by  the  extra  couplant  layer  in  system  B,  as  far  as  can  be  seen 
from  the  resolution  of  the  plot  (about  0.01  psec) .  There  is, 
however,  a  very  slight  reflection  produced  at  the  couplant  layer — 
it  shows  up  at  5.5  usee  in  the  output  signal  -for  system  B.  This 
is  consistent  with  the  calculated  value  of  1.02  usee  for  the 
transit  time  of  a  0.635  cm  (0.250  in)  aluninum  plate. 

The  results  of  this  experiment  further  reinforce  the  earlier 
conclusion  that  the  presence  of  a  couplant  layer  in  an  ultrasonic 
system  does  not  introduce  additional  time  delays  into  the  system 
response.  And,  finally,  although  slight  reflections  are  generated 
at  impedance-matched  Interfaces  where  oouplant  is  used,  the 
quantitative  impedance  model  for  wave  interaction  at  surfaces  is 


essentially  unaffected  by  the  praaanea  of  a  oouplant  layer. 

H.3  THE  EFFECT  OF  ADHESIVE  BONDING  LAYERS 

An  experiment  similar  to  the  previous  one  has  been  performed 
using  an  adhesive  bonding  layer  between  the  two  0.635  an  (0.250 
in)  a 1  uni nan  plates,  instead  of  a  couplant  layer.  The  measured 
broadband  pulse  response  of  this  system  is  shown  in  Fig.  H3.  The 
output  for  a  solid  1.270  an  (0.500  in)  al uninun  plate  is  al9o  shown 
for  comparison. 

The  output  for  the  adhesive-bonded  specimen  shows  some 
distortion  that  has  not  been  observed  in  any  of  the  previous 
experiments.  The  first  peaks  are  slightly  distorted  near  3  ysec, 
and  the  trailing  edge  of  the  largest  peak  is  broadened  as  well. 

The  particular  causes  of  these  distortions  are  difficult  to 
isolate,  but  perhaps  the  epoxy  layer  has  a  transit  time  or  a  time 
delay  that  is  close  enough  to  the  piezoelement  transit  time  to 
interfere  with  the  arrival  of  the  first  three  wavefronts  A(ll), 
A(21) ,  and  A(31)  (see  Fig.  E3) .  Regardless  of  this  distortion, 
however,  there  are  three  clear  observations  that  can  be  made. 

First,  no  additional  time  delay  between  the  application  of 
the  input  voltage  and  the  arrival  of  the  first  output  peak  is 
introduced  by  the  epoxy  layer.  Second,  a  significant  reflection 
shows  up  in  the  measured  output  at  5.5  usee.  Since  the  calculated 
transit  time  of  the  0.635  cm  (0.250  in)  plate  is  2.04  usee,  and 
the  transit  time  of  the  wearplate  is  less  than  0.06  usee,  this 
reflection  is  a  consequence  of  waves  bouncing  back  and  forth 
between  the  epoxy  layer  and  the  piezoelements.  No  apparent  time 
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delay  is  associated  with  these  reflections,  either. 

*  Finally,  there  is  a  difference  in  amplitude  between  the  two 
output  signals  in  Fig.  H3.  The  amplitude  of  the  largest  peak 
(3.8  ysec)  for  system  B  is  only  0.48  of  the  corresponding  value 

*  for  system  A.  This  factor  may  be  used  to  calculate  the  impedance 

4 

of  the  epoxy  from  eqns.  (2-17) ,  giving  a  value  Z  «2.0  x  10 

2 

gm/cm  /s.  This  value,  an  approximate  wavespeed,  and  the  measured 
thickness  of  the  epoxy  layer  (0.005  cm)  could  be  used  to  attempt 
to  simulate  the  output  of  this  system  using  the  computer  model. 

In  conclusion,  the  presence  of  an  adhesive  bond,  however 
thin,  does  have  a  significant  effect  on  the  interaction  of  stress 
waves  at  an  interface.  In  addition  to  generating  multiple 
reflections,  the  epoxy  layer  reduces  the  amplitude  of  the  waves 
transmitted  across  the  interface.  For  the  particular  specimen 
tested  in  this  study,  no  time  delays  were  observed  at  the  epoxy 
layer,  either  for  reflection  or  transmission. 
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Comparison  of  broadband  pulse  response  for  (a)  1.270  cm  (0.S00  In) 
aluminum  plate  and  (b)  two  0.635  cm  (0.250  in)  aluminum  plates 
with  intermediate  couplant  layer  clamped  between  AE-101  transducers. 


Comparison  of  broadband  pulse  response  for  (a)  1.270  cm  (0.500  in) 
aluminum  plate  and  (b)  two  0.635  cm  (0.250  in)  aluminum  plates 
with  intermediate  epoxy  layer  clamped  between  AE-101  transducers. 
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APPENDIX  I 

SIMULATION  OF  FREQUENCY-DEPENDENT  BACKING  ATTENUATION 
USING  A  LOW-PASS  DIGITAL  FILTER 

In  this  appendix,  a  scheme  is  proposed  for  the  efficient 
simulation  of  the  strongly  frequency-dependent  attenuation  of  the 
backing  material  used  in  the  transducers  in  this  study.  This 
enables  the  prediction  of  the  backing  wave  produced  by  a  broadband 
input  signal. 

1.1  PRELIMINARY  PREDICTION  OF  BACKING  WAVE 

Fig.  II  shows  the  complete  measured  and  predicted  broadband 
pulse  response  of  two  EC-500  transducers  clamped  face-to-face. 

The  predicted  output  signal  is  based  on  the  parameters  given  in 
Table  2  for  a  0.4  MHz  tone  burst  input.  This  is  only  a  crude 
approximation,  since  the  input  pulse  has  appreciable  frequency 
components  from  0. 0-1.0  MHz,  and  Fig.  C6  shows  that  the  backing 
attenuation  varies  significantly  with  frequency  at  the  high  end 
of  that  range. 

The  predicted  output  signal  shows,  as  expected,  a  backing 
wave  arriving  at  about  10  usee  after  the  application  of  the  input 
voltage.  It  also  shows  a  second  backing  wave  arriving  at  about 
20  usee,  with  greatly  reduced  amplitude.  The  time  delays  for 
each  of  these  backing  waves  appear  to  be  correct,  but  the 
amplitude  and  shape  are  distinctly  different  from  those  in  the 
measured  output  signal,  observe  that  th  t  first  backing  wave  in 


H*«viou*2fe« 


-218' 


the  predicted  signal  is  merely  a  reduced-amplitude,  inverted 
version  of  the  primary  wave,  while  the  measured  output  shows  a 
very  smooth  backing  wave,  without  the  jagged  peaks  of  the  primary 
wave.  This  behavior  can  be  simulated  by  selective  filtering  of 
the  high-frequency  components  in  the  predicted  backing  wave. 


1.2  DESIGN  CP  TOE  FILTERING  ALGORITHM 

A  one-parameter  digital  low-pass  filter  is  described  by 

T  (t)  »  (1  -  A)T  (t)  +  (A)T  (t  -  At) 
o  i  o  (11) 

where  T  represents  the  stress  wave  entering  the  filter,  T 

1  V 

represents  the  filtered  stress  wave  leaving  the  filter,  At  is  the 
time  interval  of  the  filtering  process  and  A  is  a  dimensionless 
filtering  parameter.  This  filter  resists  arbitrarily  abrupt 
(high-frequency)  changes  in  the  stress  by  weighting  the  previous 
value  of  stress  into  the  current  value,  so  that  the  filtered 
output  stress  is  a  sort  of  time  average  of  the  input  and  output  at 
a  particular  instant.  Setting  A«0  means  that  TQ(t)  *  T^t),  so 
that  no  filtering  is  introduced  at  all.  For  A*l,  on  the  other 

hand,  T  (t+At)*T  (t),  and  the  value  of  T  remains  constant 

O  0  o 

regardless  of  variations  in  T^.  Thus,  increasing  A  from  0  to  1 
varies  the  cutoff  frequency  from  infinity  all  the  way  down  to  0. 
Note  that  the  cutoff  frequency  depends  on  At  as  well  as  on  the 
filtering  parameter,  A. 

Fig.  12  shows  the  predicted  broadband  pulse  response  of  the 
system  for  several  different  values  of  A.  Of  these  four  plots, 
curve  (d) ,  for  which  A-0.975,  appears  to  look  the  most  like  the 


measured  output  signal 


1.3  PREDICTION  OF  THE  BACKING  WAVE  USING  TOE  FILTER 

Fig.  13  stews  the  final  prediction  of  the  complete  output 
signal,  including  the  filtered  backing  wave.  The  results  are 
very  encouraging,  with  the  predicted  and  measured  backing  waves 
showing  excellent  agreement  with  one  another. 

Fig.  14  is  a  plot  of  the  attenuation  of  the  filter  versus 
frequency.  Comparing  this  plot  to  Fig.  C6  (measured  backing 
attenuation  versus  frequency)  stews  that  the  predicted  curve  is 
much  steeper  than  the  measured  curve,  and  that  the  predicted 
cutoff  frequency  is  somewhat  low.  Probably  a  two-  or  three- 
parameter  filter  is  required  to  accurately  duplicate  the  measured 
frequency-dependence  of  the  backing  attenuation. 
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Fig.  12  Predicted  broadband  pulse  response  for 
two  FC-500  transducers  clamped 
face-to-face  as  a  function  of  attenuation 
filtering  parameter. 
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Fig.  14  Plot  of  attenuation  vs.  frequency  for 
the  one-paraaeter  digital  low-pass 
filter  00.975). 
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